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INTRODUCTION 

Background 

Breast  cancer  is  the  leading  cause  of  cancer  deaths  among  women  in  many  parts  of  the  world  and 
currently  no  single  imaging  modality  has  both  high  sensitivity  and  specificity  for  breast  cancer.  Although 
x-ray  mammography,  the  current  gold  standard  for  breast  imaging,  has  high  sensitivity,  false-positive 
mammograms  have  a  substantial  cumulative  risk.  New  imaging  modalities  with  better  specificity  will 
reduce  these  risks. 

Electrical  properties  of  malignant  tissues  are  significantly  different  from  the  properties  of  normal  and 
benign  tissues.  It  has  been  reported  that  the  electrical  impedance  of  malignant  tumors  decreases  by  a 
factor  of  20  to  40  with  respect  to  the  normal  or  benign  tissue  [1].  Therefore,  electrical  conductivity 
information,  which  is  inversely  related  to  impedance,  may  be  used  in  tumor  detection  and 
characterization.  At  the  present  time,  well-established  breast  screening  methods  have  high  sensitivity  but 
suffer  from  poor  specificity.  Specifically,  x-ray  mammography,  which  is  accepted  as  gold  standard  for 
breast  cancer  screening,  provides  high  sensitivity  but  has  a  high  rate  of  false-positives  [2-4],  Similarly, 
although  the  sensitivity  of  breast  MRI  has  been  extremely  high,  the  specificity  for  the  detection  of 
abnormalities  is  variable  [5,  6],  In  addition,  these  techniques  fail  to  detect  the  breast  cancer  in  some  cases, 
such  as  dense  breast  tissue,  which  is  common  among  younger  patients,  and  patients  undergoing  hormone 
replacement  therapy  during  post-menopausal  period  [7],  Accurate  measurements  of  conductivity  can  be 
used  to  detect  tumors  and  achieve  higher  specificity  rates  compared  to  the  currently  used  techniques. 

Electrical  Impedance  Tomography  (EIT)  has  been  developed  in  1980s  to  reconstruct  the  conductivity 
distribution  inside  a  volume  conductor.  In  this  method,  a  current  distribution  inside  the  object  is  generated 
by  injection  or  induction  and  peripheral  voltage  measurements  are  acquired  to  find  the  internal 
conductivity  distribution  [8].  Although  it  is  widely  used  for  dynamic  imaging,  it  has  been  reported  that 
the  technique  has  potential  application  in  breast  cancer  detection  [9]. 

Electrical  Impedance  Spectroscopy  (EIS)  is  a  similar  technique  developed  to  diagnose  malignant 
breast  tumors  [10-14],  There  is  an  FDA  approved  device  that  employs  this  technique  [15],  However, 
conductivity  images  reconstructed  using  this  technique  are  of  low  spatial  resolution  and  are  insensitive  to 
cases  where  tumor  is  located  3-3. 5cm  away  from  the  surface  or  just  under  the  nipple.  Moreover,  the 
spatial  resolution  is  not  uniform  throughout  the  imaging  region. 

MREIT  is  a  recently  developed  method  for  conductivity  imaging  [16-21],  It  uses  the  magnetic  flux 
density  measurements  acquired  from  MR  phase  images  in  conductivity  reconstruction.  Magnetic  flux 
density  generated  due  to  applied  currents  can  be  measured  with  high  and  space  independent  spatial 
resolution  using  MRI  for  DC  [22],  AC  [23,  24]  and  RF  [25]  currents.  Note  that,  only  the  component  of  the 
magnetic  flux  density  in  the  direction  of  the  main  field  of  the  MRI  system  can  be  measured.  Then,  the 
inverse  problem  of  finding  conductivity  or  current  density  from  magnetic  flux  density  can  be  solved.  Only 
relative  conductivity  values  can  be  reconstructed  if  the  magnetic  flux  density  measurements  are  used 
alone.  In  order  to  find  the  absolute  conductivity  values,  at  least  one  voltage  measurement  from  the 
boundary  is  required.  The  reconstruction  algorithms  can  also  be  divided  into  two  groups  depending  on  the 
data  type  required.  First  group  uses  magnetic  flux  density  directly  whereas  in  second  one,  current  density 
distribution  is  required  in  the  image  reconstruction.  Algorithms  that  uses  current  density  distribution 
requires  the  measurement  of  all  three  components  of  the  magnetic  flux  density,  thus,  rotation  of  the  object 
inside  the  magnet  is  required. 
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BODY 


In  this  multidisciplinary  post-doctoral  training  grant,  five  specific  aims  were  proposed. 

Aim  1.  Training  stage:  In  order  to  be  able  to  carry  out  the  study  proposed,  principal  investigator  will 
receive  training  in  tumor  biology  and  DCE-MRI. 

Aim  2.  Development  of  3D  reconstruction  algorithms  for  MREIT:  In  parallel  with  the  training  stage,  new 
studies  to  improve  the  MREIT  images  will  be  carried  out. 

Aim  3.  Test  of  3D  reconstruction  using  phantoms:  Phantom  experiments  will  be  carried  out  for  the  3D 
case. 

Aim  4.  Performing  animal  experiments  using  both  techniques:  Comparative  animal  experiments  using 
ENU  (N-Ethyl-N-Nitrosourea)  induced  malignant  and  benign  breast  tumors  in  rats  will  be  carried  out. 
Aim  5.  Statistical  analysis:  Efficacy  of  MREIT  and  comparison  with  DCE-MRI  will  be  evaluated  by  ROC 
analysis  of  the  acquired  data  and  results  from  the  pathological  examination. 

Statement  of  Work 

The  tasks  to  achieve  the  specific  aims  listed  above  are  covered  in  the  statement  of  work  as  outlined 
below. 

Task  1.  Acquire  the  necessary  training  in  tumor  biology  and  dynamic  contrast-enhanced  magnetic 
resonance  imaging  (Months  0-24) 

a)  Auditing  courses  and  attending  seminars  on  tumor  biology  (Months  0-12). 

b)  Training  in  DCE-MRI  and  preparation  of  animals  for  the  training  experiments  (Months  13-16). 

c)  Experimental  training  in  DCE-MRI  (Months  17-21). 

d)  Statistical  analysis  of  DCE-MRI  data  (Months  22-24). 

Task  2.  Implement  a  3D-reconstruction  algorithm  for  MREIT  and  test  with  simulations  (Months  0-18). 

a)  3D  Mesh  generation  from  MR  slice  images  (Months  0-2). 

b)  Implementation  of  3D  FEM  algorithm  for  arbitrary  meshes  (Months  3-6). 

c)  Implementation  of  the  reconstruction  algorithms  and  testing  with  simulation  data  (Months  6-10). 

d)  Incorporation  of  a  priori  anatomical  information  in  image  reconstruction  (Months  11-12). 

e)  Determination  of  optimal  electrode  locations  for  arbitrary  meshes  (Months  13-14). 

f)  Making  necessary  updates  in  the  algorithm  (if  necessary)  by  assessing  the  performance  with  the 
preliminary  experimental  results  (Months  15-18). 

Task  3.  Modify  existing  hardware  and  pulse-sequences  and  test  3D  reconstruction  with  phantoms  and 
animals  (Months  13-24). 

a)  Construction  of  3D  conductivity  phantoms  (Month  13). 

b)  Modifying  and  testing  the  pulse  sequence  for  multi-slice  imaging  (Months  14-15). 

c)  T esting  the  system  and  adjusting  the  imaging  parameters  using  phantoms  (16-18). 

d)  Carrying  out  animal  experiments  using  MREIT  to  determine  the  optimal  parameters  and  electrode 
locations  (Months  19-24). 

Task  4.  Perform  animal  experiments  using  both  techniques, 
a)  Carrying  out  experiments  on  50  ENU  induced  rats  using  both  techniques  for  hypothesis  testing. 
(Months  25-32). 

Task  5.  Receiver  Operating  Characteristics  (ROC)  analysis. 

a)  Analysis  of  MREIT  and  DCE-MRI  data  (Months  25-32). 

b)  Comparison  with  pathology  results  and  ROC  analysis  (32-36). 
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Task  1  Training  in  tumor  biology  and  dynamic  contrast-enhanced  magnetic  resonance  imaging 

As  stated  in  the  description  of  the  research  training  section  of  the  proposal,  the  principal  investigator 
attended  a  basic  cancer  biology  course  (UCI,  BioSci  125  Molecular  Biology  of  Cancer)  offered  by  her  co¬ 
mentor  Dr.  Eva  Lee. 

For  training  in  tumor  biology  and  applications  of  different  imaging  modalities  in  breast  cancer 
imaging,  the  principal  investigator  attended  several  seminars.  One  of  the  comments  in  the  peer  review 
panel  summary  was  the  lack  of  expansion  of  the  training  in  tumor  biology  in  other  imaging  modalities. 
This  comment  was  taken  into  account  in  the  selection  of  seminars  attended.  The  detailed  list  of  seminars 
were  given  in  previous  reports. 

Experimentation  and  analysis  methods  in  dynamic  contrast  enhanced-magnetic  resonance  imaging 
(DCE-MRI)  were  also  part  of  the  training  in  the  second  year  of  the  award.  Principal  investigator  acquired 
training  in  data  acquisition  (adjustment  of  dynamic  pulse  sequences,  synchronization  with  contrast  agent 
injection)  and  performed  experiments  herself. 

Task  2.  Implementation  a  3D-reconstruction  algorithm  for  MREIT  and  test  with  simulations 
Forward  Problem  Formulation 

Calculation  of  surface  potential  values  and  inner  magnetic  flux  density  distribution  for  a  known 
conductivity  distribution  and  boundary  conditions  is  called  the  forward  problem  of  MREIT.  The  relation 
between  conductivity  and  potential  distributions  in  the  object  is  given  by  the  Laplace's  equation, 

V(a(x,y,z)V^(x,y,z))=0  (. x,y,z)eV 

where  cr(x,  y,  z )  is  electrical  conductivity  distribution,  </>(x,  y,  z )  is  the  electrical  potential  distribution  and 
V  is  the  domain  representing  the  object  to  be  imaged.  For  the  case  of  current  injected  on  the  surface  via 
two  current  electrodes  Neumann  boundary  condition  holds,  which  is, 

j  positive  electrode 
-  j  negative  electrode 
0  elsewhere 

where  j  is  the  magnitude  of  the  applied  current  and  d(j)  /  dn  is  the  change  in  electrical  potential  in  the 
direction  of  the  outward  normal.  In  order  to  calculate  the  magnetic  flux  density  distribution,  B(x,  y,  z) 
inside  the  object,  current  density  distribution,  J(x,y,z )  must  be  calculated  first.  Since  the  conductivity  is 
assumed  to  be  known  in  the  forward  problem,  electric  field  intensity  distribution,  E{x,y,z),  and  the 
corresponding  j(x,y,z)  are  found  using 


r(x,y,z) 


d<f>(x,y,z) 


dn 


E{x,  y,  z)  =  -V  (f){x,  y,  z)  and  j(x,  y,  z)  =  <r(x,  y,  z)E(x,y,z ) 


Finally,  B\x,y,z)  is  calculated  using  the  Biot-Savart  relation: 


B{x,y,z) 


M>  c  j(x,y)dV  xaR 

4x>  R2 
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where  B{x ,  y,  z)  is  the  magnetic  flux  density  at  measurement  points,  //o  is  the  permeability  of  free  space, 
R  is  the  distance  between  source  and  field  points,  and  aR  is  the  unit  vector  pointing  from  source, 
(x',y,z'),  to  the  field  point,  (x,y,z).  Only  the  z-component  of  the  magnetic  flux  density  is  measured 
using  MRI  and  is  used  in  image  reconstruction  and  it  is  given  by 


Bz(x,y,z) 


Mo  r  J_x (*' » / » - /) - Jy (*' » y > - *0 
4  n*  R3 


Note  that  only  x  and  y  components  of  the  current  will  create  a  magnetic  field  in  the  z-direction  and 
therefore  are  required  in  the  formulation. 

For  most  irregular  geometries  and  complex  conductivity  distributions,  it  is  not  possible  to  solve  above 
equations  analytically.  Finite  element  method  (FEM)  is  used  to  solve  the  boundary  value  problem 
defined.  [26].  First  order  elements  are  used  to  discretize  the  imaging  volume,  V.  Each  element  is  assumed 
to  have  uniform  conductivity,  and  potential  distribution  in  each  element  is  approximated  to  vary  linearly 
within  the  element.  Consequently,  electric  field  and  current  density  inside  each  element  are  constant.  In 
the  calculation  of  magnetic  flux  density  from  current  density  using  the  Biot-Savart  relation,  each  finite 
element  is  taken  to  act  as  a  current  source  located  at  its  center  of  gravity,  and  the  measurement  points  are 
also  selected  as  the  elements  centers.  Although  magnetic  flux  values  are  calculated  at  element  centers, 
practically  magnetic  flux  density  is  measured  on  a  rectangular  grid,  therefore,  an  interpolation  is  required 
before  reconstruction. 

Mesh  Generation 

The  accuracy  of  the  image  reconstruction  algorithms  highly  depends  on  the  quality  of  the  match 
between  the  true  physical  properties  and  the  corresponding  model.  Earlier  attempts  for  using  exact 
boundary  information  in  reconstruction  included  some  transformations  from  simple  geometries  that  were 
still  coarse  approximations.  Another  approach  is  to  force  the  structure  into  a  more  regular  shape,  such  as 
placing  a  rat  in  a  cylinder  or  compress  breast  between  plates,  in  which  the  anatomical  structures  shift 
from  their  true  locations.  In  cases  where  multi-slice  anatomical  images  are  available,  it  is  possible  to 
generate  true  geometry  from  these  images. 

In  MREIT,  although  we  only  use  the  phase  maps  acquired  during  current  injection  for  conductivity 
imaging,  it  is  possible  to  record  a  set  of  high-resolution  anatomical  images  and  use  this  information  for 
geometry  formation.  Another  advantage  is  that  the  subject  is  kept  exactly  in  the  same  position  during  both 
experiments  therefore  movement  errors  are  minimum.  High-resolution  anatomical  images  are  used  not 
only  for  geometry  extraction  but  also  for  the  determination  of  exact  boundary  conditions  using  MR 
visible  markers. 

In  Figure  1,  a  set  of  rat  MR  images  with  markers  is  given.  First,  the  boundary  points  (dotted  red  line) 
are  manually  selected  for  the  first  slice.  The  number  of  points  in  all  slices  is  kept  constant  and  determined 
by  shifting  the  locations  from  the  previous  slice.  After  the  boundary  points  for  each  slice  are  determined, 
they  are  coupled  through  slices  and  using  the  loft  function  of  FEMLAB,  the  3D  geometry  in  Figure  2a  is 
formed.  The  important  step  in  geometry  formation  is  the  selection  of  appropriate  weights  and  coupling  of 
edges  in  each  slice  so  that  a  valid  geometry  can  be  formed.  Once  the  geometry  is  defined,  a  finite  element 
mesh  with  tetrahedral  elements  is  generated  as  shown  in  Figure  2b. 
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Figure  2.  (a)  3D  geometry  formed  using  boundary  points  defined  in  Figure  1.  (b)  Corresponding  finite 
element  mesh  with  24617  tetrahedrons. 
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After  the  formation  of  the  3D  geometry  and  the  corresponding  finite  element  mesh,  next  step  is  the 
solution  of  the  forward  problem  using  the  method  explained  in  the  formulation  section.  In  Figure  3a,  the 
simulated  electrical  potential  distribution  for  a  pair  of  electrodes  placed  oppositely  on  top  left  and  bottom 
right  of  the  geometry  is  given.  As  expected,  a  high  potential  is  observed  under  the  positive  electrode  and  a 
low  potential  is  observed  above  the  negative  electrode.  In  Figure  3b,  gradient  of  potential  distribution  is 
seen.  Note  that  due  to  the  negative  sign  in  equation  (3),  the  current  flow  lines  will  be  in  the  opposite 
direction  and  currents  are  flowing  from  positive  electrode  to  the  negative  electrode. 


Figure  3.  (a)  Electric  field  distribution  for  a  single  pair  of  electrodes,  (b)  Electrical  current  density  distribution. 

Iterated  Senstivity  matrix  reconstruction  using  Tikhonov  regularization 

In  the  reconstruction  process,  a  uniform  conductivity  is  assumed  as  the  initial  distribution  and  the 
forward  relation  between  conductivity  and  magnetic  flux  density  is  linearized  around  this  conductivity 
distribution  as  explained  in  detail  in  Birgul  et.al  [20],  Resulting  matrix  equation  that  approximates  the 
relation  between  conductivity  and  magnetic  flux  density  perturbations  is  given  by: 

Ab  =  SA  a 

where  Ab  is  the  difference  between  measured  magnetic  flux  density  and  the  magnetic  flux  density 
corresponding  to  initial  distribution,  Aa  is  the  change  with  respect  to  initial  conductivity  distribution  and 
S  is  the  sensitivity  matrix  that  gives  the  relation  between  changes  in  magnetic  flux  density  and 
conductivity.  The  sensitivity  matrix  is  calculated  using  semi-analytical  method  [20],  The  sensitivity 
matrix  is  ill-posed  in  general  and  due  to  noise  in  the  data,  some  form  of  regularization  is  required  to  be 
able  to  solve  the  system  in  above  equation.  Including  Tikhonov  regularization  parameter, 

(STS+AI)Aa=  STAb 

where  I  is  the  identity  matrix.  In  the  first  two  years  of  this  study,  to  choose  the  optimum  regularization 
parameter,  the  matrix  equation  was  solved  for  different  values  of  X  using  the  built-in  cgs  function  of 
MATLAB™  with  default  settings  (l(f6  tolerance  level,  20  iterations,  and  no  pre-conditioning)  and  the 
optimum  regularization  value  is  selected  as  the  one  minimizing  the  difference, 

m 

m}nY\Bmeas,i-BcalcAA\ 

X  <=1 
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where  m  is  the  total  number  of  measurement  points,  Bmeas  is  the  measured  magnetic  flux  density,  and  Bca\c 
is  the  flux  density  calculated  using  reconstructed  conductivity.  We  have  observed  that  the  regularization 
parameter  chosen  with  this  approach  was  sometimes  not  the  optimum,  especially  for  the  irregular  outside 
boundary  cases  and  we  started  using  the  following  criteria  for  the  selection  of  X 

A  =  trace(STS)/« 

where  trace  is  the  sum  of  the  diagonal  elements  of  the  matrix  product  and  n  is  the  dimension  of  the  matrix 
product.  For  iterations,  calculated  conductivity  distribution  is  assigned  as  the  initial  value  and  the  steps 
starting  with  sensitivity  matrix  calculation  are  repeated  until  the  change  in  conductivity  two  consecutive 
iterations  are  below  a  defined  threshold.  The  flowchart  of  the  algorithm  is  given  in  Figure  4. 


Figure  4.  Flowchart  of  the  reconstruction  algorithm 

Compared  to  sensitivity  based  reconstruction  algorithms  reported  before,  the  major  differences  of  the 
algorithm  used  in  this  study  are: 

•  using  conjugate  gradient  solver  (CGS)  instead  of  truncated  SVD.  This  eliminated  the 
requirement  of  large  decomposition  matrices,  which  enables  faster  and  memory  efficient 
reconstruction. 

•  introducing  Tikhonov  regularization  instead  of  truncation.  This  eliminates  operator  dependence. 

•  iterating  the  algorithm  to  improve  accuracy  even  for  large  conductivity  perturbations. 

It  is  known  that  only  relative  conductivity  values  can  be  reconstructed  using  magnetic  field 
measurements  alone.  Using  the  non-linear  reconstruction  outlined  above,  the  conductivity  distribution  can 
be  found  with  a  scale  and  peripheral  voltage  measurements  can  be  used  to  calculated  this  scaling  factor 
[20].  In  this  study,  voltages  are  not  measured  and  only  relative  conductivity  distributions  are  found. 
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Effects  of  Limited  Volume  Coverage  on  Accuracy  of  MREIT 

The  goal  of  this  study  was  to  investigate  the  effects  of  limited  volume  coverage  on  the  accuracy  of 
conductivity  reconstruction  in  MREIT.  Since  the  currents  injected  into  an  object  will  be  distributed  in  the 
whole  volume,  limited  volume  coverage  in  data  acquisition  or  reconstruction  will  impact  the  accuracy  of 
resulting  conductivity  maps.  On  the  other  hand,  the  magnitude  of  current  density  will  decrease  rapidly  as 
the  distance  from  the  injecting  electrodes  increases.  Moreover,  if  one  is  interested  in  a  particular  volume 
of  interest  (VOI),  the  magnetic  fields  generated  by  currents  away  from  the  VOI  could  be  negligible. 
Therefore,  we  investigated  the  decrease  in  magnitude  of  current  density  as  the  distance  from  the 
electrodes  increases.  We  also  investigated  the  contribution  of  these  weaker  currents  to  the  magnitude  of 
the  magnetic  fields  generated  within  the  VOI. 

We  performed  several  simulation  studies  to  calculate  the  current  density  and  magnetic  fields  inside  a 
cylindrical  volume  (d  =  4.5cm  h  =  7cm)  with  uniform  conductivity  distribution.  Two  electrodes  were 
placed  at  135°  and  315°  around  the  central  transaxial  slice  (slice  0).  In  3D  FEM  model  13773  nodes  and 
73452  of  tetrahedral  elements  were  used.  Figure  5a  shows  the  magnitude  of  the  transverse  current  density 
(./x(r)2+  Jy(r)2)12  across  15  slices  for  regions  labeled  13-19  in  Figure  5b.  bz(r)  was  calculated  for  five 
cases  on  these  15  transverse  planes  that  were  5mm  apart.  For  each  case,  the  current  density  only  within  a 
transverse  slab  of  various  thicknesses  was  used  for  bz(r)  calculations  as  given  in  Table  1. 


Table  1.  Description  of  cases. 

Case  1  Case  2  Case  3  Case  4  Case5 


Figure  6.  bz(r)  contour  maps  for  case  1  (a),  case  5  (b)  are  illustrated.  The  difference  of  the  two  bz(r) 
maps  is  shown  in  (c).  Current  was  injected  between  electrodes  placed  at  135°  and  315°. 
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Figure  6  shows  the  bz(r)  contour  maps  in  slice  0.  Figure  6a  is  the  resulting  magnetic  field  when  the 
transverse  current  density  in  4.2cm  thick  slab  was  used  (Case  1).  Similarly,  Figure  6b  is  the  bz(r)  for  Case 

5.  Figure  6c  is  the  difference  field  when  6b  is  subtracted  from  6a.  When  Case  1  is  taken  as  the  reference, 
the  field  difference  maps  (Abz(r))  between  Case  1  and  others  give  us  a  measure  of  the  errors  made  in  the 
calculation  of  bz(r)  when  the  current  density  outside  the  selected  slab  is  ignored.  Therefore,  we  defined  a 
region  of  interest  (ROI)  that  encompassed  the  white  contour  band  in  the  upper  right  hand  comer  of  Figure 

6.  The  mean  bz(r)  inside  the  ROI  for  Case  1  is  taken  as  the  reference  and  the  mean  of  Abz(r)  for  each  case 
is  calculated  and  divided  by  this  reference  mean.  Results  are  summarized  in  Table  2. 


Table  2.  Normalized  mean  valued  with  respect  to  Case  1. 

Case  2  Case  3  Case  4 
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The  result  in  Figure  5  shows  that  the  current  density  in  the  plane  that  is  21mm  away  from  the 
electrodes  has  significant  magnitude,  (e.g.  in  region  19,  the  magnitude  of  the  current  is  roughly  40%  of 
that  of  slice  0).  Moreover,  the  magnetic  field  in  the  ROI  in  slice  0  is  underestimated  by  almost  30%  if 
only  currents  within  a  12mm  slab  are  taken  into  account.  This  will  lead  to  underestimation  of  conductivity 
maps  in  this  slice.  Therefore,  one  has  to  acquire  MREIT  data  that  covers  sufficiently  large  volume  to 
obtain  accurate  conductivity  maps.  On  the  other  hand,  reconstruction  of  conductivity  maps  from  large 
data  sets  may  require  18-20  hours  of  computing  time.  Therefore,  optimum  volume  coverage  has  to  be 
found  to  obtain  a  balance  between  acceptable  accuracy  and  computational  efficiency  by  carrying  out 
simulations  and  experiments.  The  results  presented  here  were  obtained  from  a  uniform  conductivity 
phantom.  The  results  from  objects  of  nonuniform  conductivity  will  be  different  and  models  should  be 
developed  accordingly  to  find  the  optimum  volume  coverage  for  studies  such  as  in  vivo  experiments. 

Conductivity  Maps  Reconstructed  with  Partial  B-field  Coverage  in  MREIT 

In  MREIT,  ideally,  the  whole  3D  volume  should  be  imaged  and  current  density  and  resulting  magnetic 
fields  have  to  be  calculated  within  this  whole  volume;  however,  due  to  the  limitations  of  data  acquisition 
and  computationally  expensive  3D  numerical  models  for  reconstruction,  calculations  are  done  within  a 
limited  volume,  especially  for  large  objects.  In  previous  section,  we  investigated  the  effects  of  calculating 
current  density  in  a  limited  volume  on  the  accuracy  of  conductivity  maps.  We  also  investigated  the  impact 
of  reconstruction  of  conductivity  maps  from  magnetic  fields  calculated  within  a  limited  slab.  If  one  is 
interested  in  a  particular  VOI,  the  currents  are  injected  through  electrodes  placed  around  the  VOI.  The 
current  density  and  resulting  magnetic  fields  will  be  most  dominant  inside  the  VOI  and  diminish  further 
away.  Therefore,  one  has  to  find  the  optimum  volume  coverage  and  field  calculations  to  find  the  best 
trade-off  between  accuracy  and  efficiency.  This  will  complement  our  previous  study  where  optimum 
volume  coverage  for  current  density  calculations  were  sought,  because  optimum  volume  coverage  for 
current  density  calculations  and  magnetic  field  calculations  could  be  different  and  both  will  have  impact 
on  accuracy  differently. 

We  performed  several  simulation  studies  to  reconstruct  the  conductivity  maps  inside  a  cylindrical 
volume  (d=4.5cm  h=4cm)  with  two  inserted  objects,  (d=1.5cm  h=1.2cm  each,  one  is  placed  at  z=0,  the 
other  at  z=1.2).  Four  electrodes  were  placed  90°  apart,  around  the  central  transaxial  slice.  Same  3D  FEM 
mesh  defined  in  previous  section  is  used.  Current  density  calculation  was  made  inside  a  transaxial  slab 
that  extended  from  -1.2cm  to  1.2  cm.  For  all  cases  reported  here,  the  current  in  all  7  slices  were  used  to 
calculate  bz(r)  so  that  the  variability  introduced  by  changing  the  current  density  volume  was  eliminated; 
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but  various  subsets  of  the  bz(r)  were  used  in  reconstruction  and  its  impact  on  resulting  conductivity  maps 
was  investigated. 

The  first  row  of  images  in  Figure  7  shows  the  magnitude  of  the  actual  conductivity  distribution  in  6  of 
the  7  slices  within  the  24mm  slab.  Current  was  injected  into  the  third  slice  from  left.  The  second  row 
shows  the  case  where  only  bz(r)  inside  three  slices  was  used  in  conductivity  image  reconstruction; 
similarly,  the  third  row  is  the  conductivity  map  when  the  fields  in  all  7  slices  were  used. 

It  is  observed  from  the  second  row  of  Figure  7  that  if  the  fields  inside  a  limited  volume  is  used,  the 
high  conductivity  object  in  the  current  injection  slice  becomes  dominant  and  replicates  itself  in  all 
surrounding  slices  and  the  second  object  (seen  in  slices  5  and  6  in  the  top  row)  becomes  almost 
completely  invisible.  When  magnetic  fields  inside  a  larger  volume  are  used  in  reconstruction,  both  objects 
reconstruct  reasonably  well.  Although  we  used  the  current  density  and  the  magnetic  field  in  only  a  24mm 
portion  of  the  whole  volume  of  40mm,  both  high  conductivity  objects  are  detected.  This  has  important 
implications  in  applications  such  as  tumor  imaging  with  MREIT  because  one  focuses  in  the  tumor  area  to 
obtain  the  conductivity  distribution  inside  the  tumor  and  surrounding  tissues,  but  possible  metastasis  or 
multiple  tumor  foci  might  be  missed  if  sufficiently  large  volume  is  not  modeled  in  reconstruction. 
Although  substantial  gain  in  reconstruction  speed  will  be  obtained  with  limited  volume  coverage,  critical 
lesions  could  be  missed.  Moreover,  a  lesion  might  appear  larger  than  it  actually  is  (see  Figure  7,  row  2). 
Considering  that  reconstruction  of  conductivity  maps  from  data  sets  of  even  10-12  slices  may  require  18- 
20  hours  of  computing  time,  speed  in  reconstruction  becomes  an  important  issue.  Therefore,  one  has  to 
acquire  MREIT  data  that  covers  sufficiently  large  volume  to  obtain  accurate  conductivity  maps. 
Moreover,  optimum  volume  coverage  used  in  the  numerical  model  of  reconstruction  has  to  be  found  to 
obtain  a  balance  between  acceptable  accuracy  and  computational  efficiency  by  carrying  out  simulations 
and  experiments. 


True 

Conductivity 

Distribution 


Reconstruction 
using 
3  Slices 


Reconstruction 
using 
7  Slices 


Figure  7.  Top  row:  actual  conductivity  distribution.  One  high  conductivity  object  is  at  3  o’clock  position  (slices 
2-4)  and  the  other  at  9  o’clock  position  (slices  5-6).  2nd  and  3rd  rows  show  the  reconstructed  map  using  bz(r) 
only  in  3  and  7  slices,  respectively. 
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Task  3.  Modification  of  existing  hardware  and  pulse-sequences  and  test  3D  reconstruction  with 
phantoms  and  animals 

In  order  to  be  able  to  measure  magnetic  flux  density  at  low  amplitude  currents,  we  implemented 
various  modifications  in  our  hardware  and  pulse  sequence  and  carried  out  various  phantom  experiments 
for  optimization.  In  our  earlier  studies,  we  were  using  sine  wave  current  synchronized  with  the  first  RF 
pulse  in  a  conventional  spin  echo  sequence,  Figure  8a.  In  this  scheme,  the  pulse  sequence  was  triggering 
an  external  function  generator  to  output  a  burst  sine  wave,  making  it  more  susceptible  to  small  timing 
errors.  In  the  new  scheme,  instead  of  using  sine  wave  we  used  a  train  of  pulses,  Figure  8b,  where  pulse 
sequence  directly  controls  polarity  and  timing  of  each  pulse  sequence,  which  provides  more  accurate 
synchronization.  Another  modification  in  the  pulse  sequence  is  that  no  current  is  applied  during  RF 
refocusing  pulses.  In  older  sequence,  RF  pulses  were  applied  during  application  of  the  burst  sine  wave 
and  it  was  suspected  that  the  small  currents  during  application  of  RF  pulses  were  leading  slight  shifts  in 
the  slice  selection  location.  We  also  maintained  slice  selection  gradient  for  refocusing  pulses  throughout 
the  current  application.  This  functions  as  crusher  gradient,  dephasing  unwanted  secondary/stimulated 
echoes  generated  from  multiple  RF  pulses.  Without  larger  crushers,  resulting  images  exhibited  significant 
interference  patterns.  Our  original  sequence  collected  single  slice  images  only,  therefore,  we  extended  the 
base  sequence  for  multi-slice  acquisition  for  3D  data  collection.  We  tested  various  echo  time,  TE,  values 
and  number  of  cycles  of  injected  current  to  determine  optimum  SNR.  Optimum  TE  is  found  at  T2  value  of 
imaged  object,  which  is  consistent  with  findings  in  the  literature  [22], 

In  addition  to  these  modifications  in  the  pulse  sequence,  we  also  constructed  an  upgraded  current 
source/multiplexer  circuit.  In  the  new  current  source,  selection  of  the  injection  profile  and  pulse  polarity 
is  controlled  by  the  pulse  sequence,  allowing  more  precise  timing  and  synchronization. 
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Figure  8.  (a)  Pulse  sequence  using  AC  current  (b)  Pulse  sequence  using  train  of  pulses. 
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3D  Reconstruction  Example 

A  3D  phantom  was  prepared  using  agarose  powder  (lgr/lOOmL)  as 
background  and  a  spherical  shell  filled  with  the  same  agarose  gel.  Due 
to  insulator  shell,  the  whole  sphere  acted  as  a  nonconducting  object. 
Four  electrodes  of  size  0.5cm  x  0.5cm  were  placed  at  the  center  of  the 
cylinder  every  90°  as  shown  in  Figure  9.  Opposite  pairs  were  used, 
giving  two  current  injection  profiles.  A  FEM  mesh  with  2399  nodes  and  <g 
11781  tetrahedral  elements  was  constructed.  The  magnetic  flux  density  o 
is  calculated  at  a  rectangular  grid  of  24x24x5,  where  2300  of  grid  £ 
points  falling  into  the  cylindrical  region.  A  total  of  4600  measurements  ® 
were  used  in  image  reconstruction.  MR  structural  images  for  the  5 
slices  used  are  given  in  Figure  10  (top  row).  Corresponding 
reconstructed  conductivity  images  are  presented  in  bottom  row,  dark 
region  corresponding  to  insulating  object. 


Figure  10.  MR  anatomical  images  (top  row)  and  corresponding  reconstructed  conductivity  images  (bottom  row) 

With  this  study,  we  demonstrated  the  feasibility  of  using  3D  sensitivity  matrix  reconstruction  in 
multi-slice  magnetic  resonance  electrical  impedance  imaging.  Although  sensitivity  matrix  based 
reconstruction  methods  are  computationally  expensive,  they  perform  well  in  the  presence  of  noise. 
Previously,  we  used  sensitivity  matrix  based  reconstruction  using  2D  FEM.  That  approach  ignores  the 
contributions  of  currents  in  off-center  slices  leading  to  various  artifacts  in  the  conductivity  images; 
conductivity  values  in  the  reconstructed  slice  could  be  under  or  overestimated,  spatial  resolution  and 
point  spread  function  could  be  compromised.  In  some  cases  the  ‘ghosts’  of  objects  in  other  slices  may 
even  appear  in  the  reconstructed  slice  depending  on  how  the  conductivity  and  resulting  3D  current  density 
is  distributed  in  the  object. 

4T  versus  7T  Experiments 

The  magnetic  flux  density  is  extracted  from  MRI  phase  images  and  one  of  the  challenges  in  MREIT  is 
acquiring  phase  information  with  sufficient  SNR  while  maintaining  the  injected  currents  at  safe  levels. 
Sadleir  et.  al.  carried  out  noise  analysis  in  magnetic  flux  density  measurements  at  3T  and  11T  and  they 
suggest  that  the  noise  level  in  magnetic  flux  density  is  reduced  as  field  strength  is  increased  by  a  factor 
approximately  proportional  to  the  increase  in  the  field  strength  [27],  We  carried  out  experiments  at  4T 
and  7T  field  strengths  to  compare  the  effect  of  increased  magnetic  flux  density  SNR  in  the  final 
reconstructed  conductivity  images. 
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MREIT  experiments  were  carried  out  using  two  phantoms  with  the  same  distribution  in  4T  and  7T 
systems.  Since  ion  diffusion  alters  the  initial  conductivity  distribution  [28],  we  prepared  the  phantom 
immediately  before  each  experiment  to  minimize  these  effects  in  the  comparison.  The  pulse  sequence 
parameters  for  both  experiments  were  as:  TR  =  500ms,  TE  =  60ms  FOV  =  7cm,  slice  thickness  =  2mm, 
field  image  matrix  =  64x64,  NEX  =  4.  The  sampling  bandwidth  was  set  to  33.3KHz  and  1.25MHz  for  4T 
and  7T  systems,  respectively.  The  amplitude  of  the  injected  current  was  1mA  and  the  duration  of  the 
positive  and  negative  cycles  were  27.3ms  and  27.7ms.  Reconstructed  conductivity  images  and  profiles 
along  p-line  (as  defined  in  Figure  1 1)  at  both  field  strengths  are  given  in  Figure  12.  Although  the  phantom 
is  slightly  rotated  between  two  cases,  this  does  not  affect  the  relative  positions  of  the  objects  with  respect 
to  each  other.  We  calculated  full-width-at-half-maximum  (FWHM)  for  both  objects  through  x  and  y 
cross-sections  and  also  we  looked  at  peak  contrast  in  the  object.  Note  that  the  SMM  reconstructs  relative 
conductivity  values  if  magnetic  flux  density  measurements  are  used  alone  and  the  values  presented  here 
are  normalized  to  the  background  conductivity.  The  peak  reconstructed  contrast  and  FWHM  values  are 
summarized  in  Table  3  for  both  field  strengths  and  objects.  The  objects  were  not  fully  resolved  in  the  p- 
line  profile  for  4T  case  but  more  separable  in  the  7T  case.  For  all  cases,  the  constructed  contrasts  are 
lower  than  the  expected  ones,  especially  for  the  object  located  at  the  center. 


(c)  (d)  Position  (cm) 


Figure  12.  (a)  image  at  4T  (b)  4T  profile  (along  p-line) 
(c)  image  at  7T  (d)  7T  profile  (along  p-line) 

Table  3.  Comparison  of  reconstructed  peak  contrast  and  FWHM 
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When  the  reconstructed  images  are  compared,  it  is  seen  that  given  the  same  amount  of  experimentation 
time,  the  accuracy  and  the  spatial  resolution  of  the  conductivity  images  at  7T  is  better  compared  to  4T. 
Although  the  reconstructed  contrast  values  are  lower  than  the  expected  values  due  to  ion  diffusion  in  all 
cases,  the  comparison  between  different  field  strengths  shows  13%-22%  improvement  at  7T  compared  to 
the  4T  case. 
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Tasks  4  &  5  Animal  experiments  and  ROC  analysis 
Preliminary  Animal  Experiments 

After  testing  our  new  system  with  phantom  studies,  we  performed  experiments  using  tumor  bearing 
rats.  These  tumors  were  either  induced  by  the  carcinogen  ENU  or  R3230  AC  tumor  grafts.  Before  starting 
the  animal  experiments  for  hypothesis  testing,  we  scanned  14  animals  in  order  to  setup  the  system  and 
define  imaging  parameters.  Since  R3230  tumors  take  a  couple  of  weeks  to  grow,  we  preferred  that  line 
line  in  development  stage  instead  of  ENU  line,  which  is  expected  to  grow  in  3-4  months.  For  animal 
imaging,  a  special  animal  holder  was  prepared  from  acrylic  sheets.  In  this  holder  electrodes  were  placed 
on  acrylic  hollow  tubes  filled  with  C11SO4  solution  mark  the  electrode  positions  precisely  in  the  images. 
The  precise  localization  of  electrode  positions  is  critical  for  accurate  assignment  of  boundary  conditions 
and  thus,  effects  the  reconstruction  performance.  Current  carrying  wires  ran  along  these  tubes,  which 
were  in  z-direction.  This  is  essential  to  minimize  interference  from  the  magnetic  fields  generated  by 
current  in  the  wires.  The  animal  was  anesthetized  by  injection  of  ketamine  and  xylazine  and  placed  inside 
the  holder.  The  skin  areas  of  contact  were  shaved  to  guarantee  good  contact.  An  anatomical  image  was 
collected  using  FSE  sequence  prior  to  the  MREIT  images.  A  single  slice  with  6 mm  thickness  was 
collected  from  the  same  anatomical  location  as  the  MREIT  image.  The  data  matrix  was  256  x  256,  FOV  = 
10cm,  Tr  =  4sec,  Te  =  40msec,  and  NEX  =  4  (signal  averages).  MREIT  images  were  collected  using  the 
pulse  sequence  in  Figure  9a  with  Tr  =  500mses,  Te  =  30msec,  NEX  =  8,  64  x  64  data  matrix,  field  of 
view,  FOV  =  10cm,  6mm  slice  thickness,  with  an  AC  current  of  1mA  peak,  100Hz  and  4  cycles. 
Reconstructed  images  from  six  animals  are  shown  in  Figure  13. 


Figure  13.  Anatomical  MR  and  corresponding  conductivity  images  for  six  different  animals. 

In  each  image  pair,  top  image  shows  the  anatomical  MR  image  and  bottom  image  shows  the 
corresponding  conductivity  image.  Circled  areas  show  the  tumor  location.  For  each  animal,  region  of 
interests  were  selected  over  the  tumor  region  and  over  the  rest  of  the  body.  Since  these  images  provide 
only  relative  conductivity  measures,  the  ratio  of  mean  conductivity  in  the  tumor  divided  by  the  mean 
conductivity  in  normal  tissue  was  calculated  for  each  animal.  The  mean  of  these  ratios  across  all  animals 
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was  p  =  2.26,  <7=  0.62.  In  other  words,  the  mean  conductivity  was,  on  the  average,  2.26  times  higher  than 
the  conductivity  in  normal  regions. 

The  major  challenge  we  encountered  in  animal  experiments  was  the  errors  in  the  magnetic  flux  density 
measurements  due  to  motion  artifacts.  The  motion  artifact  in  the  MR  image  is  clearly  seen  in  top  left 
image  pair  in  Figure  13.  For  some  cases,  this  effect  was  worse  and  the  data  had  to  be  discarded. 

ENU  Animal  Preparation  for  Comparative  Studies 

Due  to  motion  artifact  problems  we  have  encountered  in  animal  imaging,  there  has  been  some  delay  in 
tumor  induction  in  animals  for  comparative  studies.  In  order  to  avoid  unnecessary  animal  experiments,  we 
first  injected  ENU  into  15  SD  rats  i.p.  on  August  9,  2006.  Tumors  are  usually  expected  to  show  up  3~4 
months  after  injection.  The  first  tumor  appeared  on  January  5,  2007.  Tumors  grew  up  around  upper  chest 
area.  In  R3230  cases,  the  tumor  location  can  be  controlled  and  all  were  in  the  flank  area  and  therefore 
imaging  was  less  susceptible  to  the  motion  due  to  breathing.  For  the  ENU  cases,  however,  tumor  location 
cannot  be  controlled  and  motion  artifacts  were  elevated.  Moreover,  air  cavities  introduced  some 
additional  problems  and  we  could  not  acquire  any  good  magnetic  flux  density  map  in  the  7  animals  that 
has  grown  tumors  using  any  of  the  spin  echo  sequences  given  in  Figure  8  that  yielded  successful 
conductivity  distribution.  At  that  stage,  since  we  didn’t  have  a  large  set  of  ENU  animals  where  we  can 
choose  a  subset  from  that  would  model  our  ultimate  application  of  breast  cancer  imaging  better,  we 
concentrated  on  improving  the  data  acquisition  scheme  by  searching  faster  imaging  alternatives  that 
would  reduce  the  problems  rising  from  motion  artifacts.  Faster  imaging  would  be  beneficial  in  both  our 
model  and  in  our  real  application  in  not  only  reducing  the  motion  artifacts  but  also  giving  the  option  of 
using  higher  number  of  averaging  that  would  make  imaging  using  lower  currents  possible.  In  the  next 
section,  the  idea  and  the  preliminary  results  for  EPI  based  image  reconstruction  for  MREIT  are  presented. 

EPI  Sequence 

As  an  alternative  to  the  established  SE-based  pulse  sequence,  we  are  now  investigating  the  use  of  a 
modified  single-shot,  spin-echo,  echo-planar-imaging  (SS-SEPI)  pulse  sequence  for  measuring  the 
current-generated  magnetic  flux  density  (Figure  14).  SS-SEPI  acquires  data  for  an  entire  image  in  a  single 
scan,  and  thus  has  the  potential  to  drastically  reduce  the  total  data  acquisition  time,  which  will  resolve  the 
problems  encounterd  in  animal  studies. 


Figure  14.  SS-SEPI  Pulse  Sequence  for  MREIT  Figure  15.  Phantom  Description 

We  tested  the  sequence  using  an  agar  phantom  first.  For  the  test  phantom,  a  hollow  acrylic  disk  with 
an  inner  diameter  of  7cm  and  thickness  of  1cm  was  filled  with  2%  agarose  and  4mM  CuS04.  Within  this 
disk,  a  smaller  cylindrical  region  of  13mm  diameter  was  filled  with  1%  NaCl,  2%  agarose,  and  4mA 
CuSC>4  to  generate  a  high  conductivity  region  (Figure  15).  The  plane  of  the  disk  was  placed  perpendicular 
to  the  main  static  MRI  field.  Two  copper  electrodes  each  6mm  wide  were  placed  opposite  of  each  other 
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along  the  inner  acrylic  wall  and  used  to  inject  currents  into  the  interior  region.  A  4mA  bipolar  current 
pulse  was  injected  into  the  phantom  and  the  resulting  magnetic  flux  density  distribution  measured  using  a 
modified  SS-SEPI  pulse  sequence  (Figure  14).  The  component  of  current-generated  magnetic  flux  density 
parallel  to  the  main  static  MRI  field  (z-component)  introduced  a  phase  shift  in  the  MR  image  given  as: 
cp(r)  =  y(TA+TB)Bz(r).  The  scan  parameters  were:  Ta  =  35ms,  Tb  =  10ms,  Tr  =  3s,  Te  =  60ms,  matrix 
=64X64,  FOV  =  20cm,  and  single  slice  thickness  =  5mm.  Data  was  collected  twice,  each  with  opposite 
polarities  in  the  applied  current  waveform.  The  resulting  phase  maps  were  subtracted  then  divided  by  two, 
so  as  to  cancel  out  any  additional  phase  contributions,  such  as  those  arising  from  small  imperfections  in 
the  hardware  timing. 


Figure  16.  Magnetic  flux  density  and  reconstructed  conductivity  distributions 

Relative  conductivity  maps  reconstructed  using  5  iterations  sensitivity  matrix  reconstruction  are  given 
in  Figure  16.  Comparison  of  the  0mA  and  4mA  cases  shows  that  the  SS-SEPI  pulse  sequence  was  able  to 
detect  the  magnetic  flux  density  generated  by  the  injected  current.  Comparison  with  the  SE  results  shows 
that  the  SS-SEPI  pulse  sequence  was  also  able  to  map  out  the  general  variation  in  the  magnetic  flux 
density  due  to  the  high  conductivity  perturbation.  However,  the  SS-SEPI  images  suffer  from  geometric 
distortions  resulting  from  the  low  sampling  bandwidth  in  the  phase  encode  direction  inherent  in  EPI  based 
pulse  sequences.  As  a  result,  the  reconstructed  conductivity  distribution  also  suffers  from  some  geometric 
distortion.  Never  the  less,  the  high  conductivity  perturbation  could  still  be  reconstructed  from  the  SS- 
SEPI  data.  We  could  further  improve  reconstruction  by  applying  some  sort  of  geometric  correction  to  the 
MRI  data.  An  advantage  of  the  SS-SEPI  pulse  sequence  is  the  significant  reduction  in  data  acquisition 
time.  A  typical  imaging  experiment  takes  about  4  minutes  when  spin  echo  based  pulse  sequence  was  used 
whereas  this  reduces  to  about  12sec  when  EPI  pulse  sequence  is  used.  This  could  allow  for  increased 
signal  averaging  in  an  allotted  study  time,  which  may  be  required  for  improving  the  SNR  when  using  low 
amplitude  injected  currents.  Single  shot  pulse  sequences  are  less  susceptible  to  motion  artifacts,  which 
will  be  of  significance  when  applying  MREIT  to  in  vivo  studies. 

Although  imaging  time  is  considerably  shorter  for  EPI-based  sequences,  they  are  prone  to  geometric 
distortions  caused  by  inhomogeneity  of  the  main  static  MRI  field.  This  is  due  to  the  extended  data 
acquisition  period  and  corresponding  low  sampling  bandwidth  in  the  phase  encode  direction.  Sources  of 
inhomogeneity  include  inherent  imperfections  of  the  main  magnet,  perturbations  caused  by  the  object  to 
be  imaged,  and  differences  in  magnetic  susceptibility  of  various  regions  including  the  object/air  interface. 
Different  correction  methods  such  as  higher  order  shimming,  pixel  shifting,  or  partial  k-space  imaging 
should  be  investigated  from  the  standpoint  of  improving  conductivity  reconstruction. 
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KEY  RESEARCH  ACCOMPLISHMENTS 

•  MREIT  system  for  3D  data  acquisition  was  implemented. 

•  Mesh  generation  and  image  reconstruction  algorithms  in  3D  were  developed. 

•  3D  phantom  studies  were  carried  out  successfully. 

•  Preliminary  animal  experiments  that  demonstrate  conductivity  contrast  for  R3230  tumor  models 
were  carried  out. 

•  Data  acquisition  system  was  generalized  for  imaging  at  different  field  strengths. 

•  Preliminary  images  were  acquired  using  fast  imaging  sequences  that  is  critical  in  in  vivo  imaging. 

•  PI  acquired  training  in  tumor  biology  and  different  imaging  modalities  in  tumor  imaging. 

CONCLUSIONS 

•  We  demonstrated  the  feasibility  of  using  3D  sensitivity  matrix  reconstruction  in  multi-slice 
magnetic  resonance  electrical  impedance  imaging.  Although  sensitivity  matrix  based  reconstruction 
methods  are  computationally  expensive,  they  perform  well  in  the  presence  of  noise. 

•  MREIT  imaging  is  possible  at  current  levels  around  1mA  and  even  lower  currents  can  be  achieved 
using  faster  imaging  sequences  and  averaging. 

•  MREIT  is  capable  of  detecting  objects  as  small  as  3mm  using  low  amplitude  currents. 

•  The  comparison  of  MREIT  experiments  carried  out  at  two  different  field  strengths  suggests  a 
slight  improvement  in  higher  field  strengths. 

•  In  vivo  animal  images  shows  that  there  is  an  increase  in  conductivity  of  tumors  and  further 
investigation  in  improvement  of  existing  hardware  and  pulse  sequences  would  overcome  the  motion 
artifacts. 
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Abstract 

Magnetic  resonance-electrical  impedance  tomography  employs  low  amplitude 
currents  injected  or  induced  inside  an  object.  The  additional  magnetic  field  due 
to  these  currents  results  in  a  phase  in  the  MR  images.  In  this  study,  a  modified 
fast  spin-echo  sequence  was  used  to  measure  this  magnetic  field,  which  is 
obtained  by  scaling  the  MR  phase  image.  A  finite  element  method  with  first 
order  triangular  elements  was  used  for  the  solution  of  the  forward  problem. 
An  iterated  sensitivity  matrix-based  algorithm  was  developed  for  the  inverse 
problem.  The  resulting  ill-conditioned  matrix  equation  was  regularized  using 
the  Tikhonov  method  and  solved  using  a  conjugate  gradient  solver.  The  spatial 
and  contrast  resolution  of  the  technique  was  tested  using  agarose  gel  phantoms. 
A  circular  phantom  with  7  cm  diameter  and  1  cm  thickness  is  used  in  the 
phantom  experiments.  The  amplitude  of  the  injected  current  was  1  mA.  3,  5  and 
8  mm  diameter  insulators  and  high  conductor  objects  are  used  for  the  spatial 
resolution  study  and  an  average  full- width  half-maximum  value  of  4.7  mm 
is  achieved  for  the  3  mm  insulator  case.  For  the  contrast  analysis,  the 
conductivity  of  a  15  mm  object  is  varied  between  44%  and  500%  with  respect 
to  the  background  and  results  are  compared  to  the  ideal  reconstruction. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

MREIT  (magnetic  resonance-electrical  impedance  tomography)  is  a  recently  developed 
imaging  modality  that  reconstructs  conductivity  images  from  the  magnetic  flux  density 
generated  due  to  a  current  distribution  in  a  volume  conductor.  The  main  advantages  of 
MREIT  are  the  non-invasive  measurement  of  the  field  from  inside  the  object  (that  increases 
the  sensitivity  to  inner  regions),  higher  number  of  measurements  without  being  limited  to 
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detection  probes  (that  improves  the  resolution),  and  availability  of  the  MR  image  for  accurate 
modelling  of  the  geometry  and  boundary  conditions. 

Although  it  is  possible  to  detect  locations  of  lesions  in  breast  cancer  using  techniques  such 
as  x-ray  mammography,  due  to  problems  of  specificity,  many  women  undergo  unnecessary 
biopsies.  Therefore,  new  highly  specific  techniques  are  needed  to  decrease  the  false  positive 
results  while  maintaining  the  high  sensitivity  (Elmore  et  al  2002).  Since  the  conductivity 
values  of  malignant,  benign  and  normal  tissues  are  significantly  different,  this  information  can 
be  used  in  tumour  classification  (Surowiec  et  al  1988)  to  improve  specificity. 

MREIT  reconstruction  algorithms  can  be  grouped  into  two  depending  on  the  type  of 
data  they  use.  Earlier  algorithms  generated  current  density  maps  from  the  measurements  as  an 
intermediate  step  towards  conductivity  reconstruction.  The  /-substitution  method  (Khang  et  al 
2002),  current-constrained-voltage-scaled  reconstruction  (CCVSR)  (Birgul  et  al  2003a)  and 
equipotential  projection-based  reconstruction  (Ozdemir  et  al  2004)  are  examples  of  this  type 
of  algorithm.  The  disadvantage  of  these  algorithms  is  the  need  for  two  or  three  components 
of  the  magnetic  flux  density,  which  requires  rotation  of  the  object  inside  the  magnet.  Other 
algorithms  use  magnetic  flux  density  measurements  directly  and  the  measurement  of  only  one 
component  is  sufficient.  Sensitivity -based  reconstruction  (Ider  and  Birgul  1998,  Birgul  et  al 
2003b,  Muftuler  et  al  2004),  the  algebraic  reconstruction  technique  (Ider  and  Onart  2004), 
variational  gradient  Bz  algorithm  (Park  et  al  2004a),  gradient  Bz  decomposition  algorithm 
(Park  et  al  2004b)  and  response  surface  methodology  (RSM)  algorithm  (Gao  et  al  2006)  fall 
into  the  second  group  of  algorithms. 

The  main  advantage  of  the  MREIT  technique  is  the  availability  of  the  volume  data. 
In  conventional  EIT,  the  number  of  peripheral  voltage  measurements  are  physically  limited 
and  therefore  it  is  a  highly  ill-conditioned  problem.  In  contrast,  the  magnetic  flux  density 
can  be  measured  from  within  the  object  with  a  resolution  sufficient  to  match  the  number  of 
unknowns  in  the  discretized  problem  resulting  in  a  better-conditioned  problem  compared  to 
the  conventional  EIT  technique. 

Several  studies  have  been  carried  out  experimentally  using  phantoms  to  evaluate  the 
performance  of  MREIT  (Birgul  et  al  2003b,  Lee  et  al  2003,  Oh  et  al  2004).  The  major 
limitation  in  practice  is  the  amplitude  of  the  current  required  to  achieve  an  acceptable  signal- 
to-noise  ratio.  Lee  et  al  (2003)  reported  relative  L2-e rrors  of  25.5%  and  32.3%  with  denoising 
and  without  denoising  with  an  injected  current  of  28  mA.  Birgul  et  al  (2003a)  carried  out 
experiments  with  25  mA  peak  current  and  reconstructed  conductivity  images  with  13-17% 
error  without  using  any  denoising  techniques  at  0.15  T.  Oh  et  al  (2004)  reported  11-35% 
error  when  the  current  was  24  mA.  Muftuler  et  al  (2004)  presented  results  for  phantoms 
and  an  in  vivo  tumour  bearing  animal  with  current  levels  that  are  as  low  as  2-4  mA  using  a 
3  T  magnet.  In  their  study,  the  ratio  of  the  conductivity  of  a  region  of  interest  to  a  reference 
was  investigated;  therefore,  no  voltage  measurements  were  acquired  resulting  in  relative 
conductivity  images.  In  their  recent  study,  Sadleir  et  al  (2006)  acquired  MREIT  images  using 
an  11  T  MRI  system.  They  have  used  current  levels  changing  between  5  and  20  mA  and 
presented  images  for  tissue  phantoms  of  different  sizes  (60  mm  and  40  mm  both  in  diameter 
and  height).  Recently,  Hamamura  et  al  (2006)  demonstrated  measurement  of  ion  diffusion 
using  the  MREIT  technique  in  a  70  mm  diameter  phantom  with  900  j± A  at  4  T.  It  must  be 
noted  that  the  field  strengths,  phantom  sizes  and  imaging  parameters  are  different  for  all  the 
studies  listed  above.  Therefore,  it  is  not  possible  to  make  a  direct  comparison  in  terms  of 
an  achievable  minimum  applied  current  level.  A  better  measure  of  comparison  would  be  the 
current  density  inside  the  object. 

In  this  study,  we  used  an  iterative  version  of  the  reconstruction  algorithm  outlined 
by  Muftuler  et  al  (2004)  with  Tikhonov  regularization.  This  reduces  noise  sensitivity  of 
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reconstruction  and  enables  one  to  use  lower  current  levels,  which  is  critical  for  human 
applications.  The  pulse  sequence  used  in  data  acquisition  was  also  improved  so  that  the 
noise  in  the  magnetic  flux  density  measurements  for  low  current  levels  were  acceptable. 
Several  agarose  phantoms  were  prepared  to  characterize  the  contrast  and  spatial  resolution 
performance.  The  details  of  data  acquisition,  image  reconstruction  and  reconstructed  images 
are  presented  in  the  following  sections. 

2.  Methods 

Reconstruction  of  conductivity  involves  two  basic  steps.  The  first  step  is  the  measurement 
of  magnetic  flux  density  using  magnetic  resonance  imaging.  This  step  involves  MRI  data 
acquisition  using  a  modified  spin-echo  pulse  sequence  and  generation  of  magnetic  flux  density 
images  from  the  MRI  phase  images  using  scaling.  In  the  second  step,  these  images  are  used 
as  input  data  in  the  inverse  problem  of  finding  conductivity  from  the  magnetic  flux  density 
information.  In  this  section  the  components  of  the  hardware  and  pulse  sequence  used  for  data 
acquisition  are  presented  first.  Then,  the  reconstruction  algorithm  based  on  the  sensitivity 
matrix  is  explained. 

2.7.  Experimental  set-up 

The  acquisition  system  for  the  MREIT  consists  of  a  conventional  whole  body  MRI  system 
and  a  current  source  that  is  synchronized  with  the  RF  pulses  generated  by  the  MRI  system. 

The  MRI  system  has  a  4  T  whole  body  Magnex  magnet  with  a  whole  body  gradient  coil 
set  that  provides  up  to  3  G  cm-1  gradient  fields.  A  13  channel  room-temperature  high-order 
shim  coil  set  was  used  to  minimize  field  inhomogeneities.  This  system  was  interfaced  with 
a  MRRS  console  (magnetic  resonance  research  systems,  Guildford,  UK)  that  has  broadband 
RF  transmit  and  receive  channels.  A  16  leg,  quadrature,  high-pass  birdcage  coil  with  10  cm 
diameter  and  18  cm  length  was  designed  and  built  in-house  for  the  MREIT  experiments. 

The  current  source  uses  pulses  generated  by  the  MRI  console  and  a  voltage -to -current 
converter  that  uses  three  EM741  OPAMPs.  This  current  source  was  triggered  by  a  TTE  pulse 
generated  by  the  scanner  computer.  The  current  cables  were  connected  to  the  object  through 
several  RF  chokes  and  low  pass  filters  to  suppress  RF  coupling.  These  cables  were  fixed  on 
an  acrylic  structure  that  also  supports  the  phantom.  The  orientation  of  the  cables  is  critical 
so  that  only  the  effect  of  current  flowing  through  the  object  creates  a  phase  in  the  MR  image. 
The  rigid  support  guarantees  that  external  currents  flow  parallel  to  the  main  magnetic  field 
(z -direction)  so  that  their  effect  is  zero. 

2.2.  Pulse  sequence 

The  MREIT  data  were  collected  using  the  pulse  sequence  shown  in  figure  1.  In  this  pulse 
sequence,  if  the  RF  pulses  are  not  synchronized  with  the  current  pulses,  the  phases  generated 
in  positive  and  negative  cycles  would  cancel  each  other  out.  By  synchronizing  successive  tt 
pulses  to  half  cycles  of  the  current,  the  phase  shift  accumulates;  however,  it  is  superimposed 
onto  the  phase  component  due  to  main  and  gradient  field  inhomogeneities.  To  extract  the 
current-only  component,  the  experiment  is  repeated  twice  with  opposite  current  polarities. 
When  the  difference  is  taken,  system-dependent  phase  shifts  cancel  out  and  current-dependent 
phase,  O/ ,  is  obtained  as 

<S>i(r)  =  2yNTvuheBz(r),  (1) 

where  y  is  the  gyromagnetic  ratio,  N  is  the  number  of  cycles  of  injected  current,  Bz(r)  is  the 
amplitude  of  the  z -component  of  the  current-generated  magnetic  field  at  point  r  and  rpuise  is  the 
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Figure  1.  Sample  pulse  sequence  used  in  the  MREIT  experiments  (parameters  used  are  given  in 
the  results  section). 


duration  of  each  current  pulse.  Measurement  of  this  phase  shift  allows  for  calculation  of  the 
z-component  of  magnetic  field  distribution.  More  information  about  extraction  of  magnetic 
flux  density  from  MR  images  can  be  found  in  Muftuler  et  al  (2004). 


2.3.  Forward  problem  and  image  reconstruction 


The  forward  problem  for  MREIT  is  the  calculation  of  peripheral  surface  voltage  values 
and  the  magnetic  flux  density  distribution  for  known  conductivity  distribution  and  boundary 
conditions.  The  nonlinear  relation  between  conductivity  and  potential  field  is  given  by  the 
boundary  value  problem  (BVP)  defined  by  Poisson’s  relation  as 

V .  (cr(x,  y)V0(x,  30)  =  0  (x,y)eS  (2) 


for  a  two-dimensional  case,  where  cr(x,y)  and  0(v,  y)  are  the  electrical  conductivity  and 
potential  field  distributions  in  the  imaging  plane,  S.  The  boundary  conditions  must  be  specified 
together  with  the  above  BVP  for  the  solution  of  the  electric  field.  The  electric  current  on  the 
boundary  of  the  imaging  region  is  specified  for  the  MREIT  problem  and  is  imposed  as  the 
Neumann  boundary  condition: 


-ct(x,  y ) 


dHx,  y) 
3  n 


J  on  positive  current  electrode 

—  J  on  negative  current  electrode 
0  elsewhere. 


(3) 


Once  the  potential  field  distribution  is  found,  the  electrical  field  distribution  can  be  calculated 
as 


E(x,  y)  =  y) 


(4) 
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and  the  corresponding  current  density  distribution  is  obtained  using  Ohm’s  relation  in  field 
quantities  as 

J(x,y)  =  cr(x,y)E(x,y).  (5) 

The  magnetic  flux  density  at  the  imaging  slice  generated  due  to  the  current  density  distribution 
found  in  equation  (5)  is  given  by  the  Biot-Savart  relation 

3,  ,  /to  (  J(x,y)dS  x  aR 

S(^)  =  sj — *5 — ■  <6) 

where  /jlq  is  the  permeability  of  the  free  space,  R  is  the  vector  from  source  point  (. x ',  yr)  to  the 
field  point  (x,  y )  and  aR  is  the  unit  vector  in  that  direction.  The  finite  element  method  (FEM) 
with  first  order  triangular  elements  is  used  for  the  discretization  of  the  forward  problem. 

In  the  reconstruction  process,  a  uniform  conductivity  is  assumed  as  the  initial  distribution 
and  the  forward  relation  between  conductivity  and  magnetic  flux  density  is  linearized  around 
this  conductivity  distribution  as  explained  in  detail  in  Birgul  et  al  (2003b).  The  resulting 
matrix  equation  that  approximates  the  relation  between  conductivity  and  magnetic  flux  density 
perturbations  is  given  by 

Ab  =  S  Act,  (7) 

where  Ab  is  the  difference  between  measured  magnetic  flux  density  and  the  magnetic 
flux  density  corresponding  to  initial  distribution,  A  a  is  the  change  with  respect  to  initial 
conductivity  distribution  and  S  is  the  sensitivity  matrix  that  gives  the  relation  between  changes 
in  magnetic  flux  density  and  conductivity.  The  sensitivity  matrix  is  calculated  using  the  semi- 
analytical  method  (Birgul  et  al  2003b).  The  sensitivity  matrix  is  ill-posed  in  general,  and  due 
to  noise  in  the  data,  some  form  of  regularization  is  required  to  be  able  to  solve  the  system  in 
equation  (7).  Including  the  Tikhonov  regularization  parameter,  A,  the  matrix  equation  becomes 

(STS  +  AI)Acr  =  ST  Ab,  (8) 


where  I  is  the  identity  matrix.  The  matrix  equation  is  solved  for  different  values  of  A  using  the 
built-in  cgs  function  of  MATLAB™  with  default  settings  ( \e~ 6  tolerance  level,  20  iterations 
and  no  pre-conditioning).  The  optimum  regularization  value  is  selected  as  the  one  minimizing 
the  difference, 

m 

min  Eh  ^meas,/  BatoAVW  (9) 

i= 1 

where  m  is  the  total  number  of  measurement  points,  7?meas  is  the  measured  magnetic  flux  density 
and  Bcdlc  is  the  flux  density  calculated  using  reconstructed  conductivity.  For  iterations,  the 
calculated  conductivity  distribution  is  assigned  as  the  initial  value  and  the  steps  starting  with 
sensitivity  matrix  calculation  are  repeated  until  the  change  in  conductivity  for  two  consecutive 
iterations  is  below  a  defined  threshold.  The  flowchart  of  the  algorithm  is  given  in  figure  2. 

Compared  to  sensitivity-based  reconstruction  algorithms  reported  before  (Birgul  et  al 
2003b,  Muftuler  et  al  2004),  the  major  differences  of  the  algorithm  used  in  this  study  are  the 
following: 

(i)  Using  conjugate  gradient  solver  (CGS)  instead  of  truncated  SVD.  This  eliminated  the 
requirement  of  large  decomposition  matrices,  which  enables  a  faster  and  memory  efficient 
reconstruction. 

(ii)  Introducing  Tikhonov  regularization  instead  of  truncation.  The  criterion  used  in  selection 
of  the  optimum  regularization  parameter  has  been  explained  before  in  this  section.  This 
eliminates  operator  dependence. 

(iii)  Iterating  the  algorithm  to  improve  accuracy  even  for  large  conductivity  perturbations. 
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Figure  2.  Flowchart  of  the  sensitivity-based  iterative  reconstruction  algorithm. 


Figure  3.  (a)  Home-built  RF  coil,  (b)  Phantom,  (c)  Phantom  description. 


It  is  known  that  only  relative  conductivity  values  can  be  reconstructed  using  magnetic 
field  measurements  alone.  Using  the  nonlinear  reconstruction  outlined  above,  the  conductivity 
distribution  can  be  found  with  a  scale  and  peripheral  voltage  measurements  can  be  used  to 
calculate  this  scaling  factor  (Birgul  et  al  2003b).  In  this  study,  voltages  are  not  measured 
and  only  relative  conductivity  distributions  are  found.  In  the  results  section,  we  scaled  the 
conductivity  values  assuming  known  voltages  to  provide  readers  a  basis  for  comparison. 

3.  Results 

3.1.  Phantom  preparation 

Various  agarose  gel  phantoms  were  used  for  the  experiments.  The  agarose  powder  (OmniPur 
brand)  was  mixed  with  different  concentrations  of  NaCl  to  generate  different  conductivity 
values.  2  g  of  agarose  powder  was  used  per  100  ml  of  water.  Gels  were  poured  into  an  acrylic 
cylinder  with  an  inner  diameter  of  7  cm  and  height  of  1  cm  (figure  3).  Four  electrodes  made 
from  6  mm  wide  copper  tape  with  conducting  glue  were  placed  at  { J }  along  the 
inner  wall,  that  could  provide  up  to  six  different  current  injection  pairs.  It  has  been  shown 
that  at  least  two  independent  current  injection  profiles  are  required  for  the  uniqueness  of  the 
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Table  1.  Measured  conductivity  values  for  different  NaCl  concentrations. 

NaCl  (g/100  mL)  0.25  1.00  1.20  2.00  4.00  7.00  10.00 

Conductivity  (S  m-1)  0.43  1.61  1.88  3.02  5.68  9.08  11.22 


solution  (Ider  et  al  2003)  and  we  selected  two  current  injection  set-ups  that  used  opposite 
electrode  pairs.  Current  carrying  wires  were  placed  precisely  in  the  z -direction  to  eliminate 
any  contribution  from  the  currents  flowing  in  those  wires.  Wires  were  mounted  on  acrylic 
support  beams  to  establish  rigidity. 

For  all  cases,  1  g  NaCl/ 100  g  water  was  used  as  the  background  solution.  Different 
amounts  of  NaCl  were  used  to  create  different  object  values.  Although  NaCl  amounts 
were  chosen  as  multiples  of  the  background  reference,  due  to  nonlinearity  the  conductivity 
values  were  not  linearly  scaled.  For  the  accurate  evaluation  of  the  reconstruction  values,  we 
measured  the  conductivity  of  all  solutions  independently  using  a  four-electrode  conductivity 
measurement  cell.  The  measured  values  are  given  in  table  1. 

3.2.  Pulse  sequence  and  current  parameters 

The  pulse  sequence  in  figure  1  was  used  with  the  following  parameters: 


repetition  time 

Tr 

500  ms 

echo  time 

Te 

58  ms 

field  of  view 

FOV 

100  mm 

number  of  averages 

NEX 

4 

slice  thickness 

Aw 

5  mm 

image  matrix 

64  x  64 

number  of  current  cycles 

N 

5 

current  peak 

h 

1  mA 

current  frequency 

f 

100  Hz 

number  of  profiles 

2. 

The  magnetic  flux  density  Bz  is  measured  on  a  5  mm  slice  centred  in  the  1  cm  agar  slab 
defined  in  figure  1(b).  Note  that  although  NEX  is  listed  as  four,  since  the  pulse  sequence  was 
applied  twice  with  opposite  polarities,  the  effective  number  of  averages  was  equivalent  to  eight. 
The  current  was  applied  at  100  Hz  with  1  mA  peak.  The  current  pulse  can  be  expanded  as  a 
summation  of  infinite  sine  waves  with  frequencies  corresponding  to  odd  harmonics  of  the  main 
frequency.  Due  to  timing  of  the  n  pulses,  only  the  component  at  the  main  harmonic  frequency 
will  accumulate  additional  phase  and  this  scale  must  be  included  in  the  calculations.  For  each 
current  profile,  data  acquisition  took  about  4  min  with  the  above  parameters  resulting  in  8  min 
total  acquisition  time.  For  phantom  studies,  we  did  not  acquire  high  resolution  anatomical 
images  prior  to  the  MREIT  experiment  and  magnitude  images  of  the  MREIT  experiment 
were  used  for  mesh  generation.  For  cases  with  irregular  boundaries,  such  as  small  animal 
imaging,  an  accompanying  high  resolution  anatomical  image  would  increase  the  accuracy  of 
the  forward  solution. 


3.3.  Reconstruction  parameters 

A  finite  element  mesh  with  1089  nodes  and  2048  triangular  elements  was  generated  for  the 
forward  solver.  6  mm  electrode  covers  two  elements  with  three  nodes  at  the  boundary.  Phase 
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Table  2.  Spatial  resolution  analysis. 


case 

FWHMx 

(mm) 

FWHMy 

(mm) 

Peak  object 
(S  m-1) 

Mean  background 
(S  m-1) 

3  mm  ins 

4.2 

5.3 

0.78 

1.58 

3  mm  cond 

8.3 

7.0 

3.01 

1.50 

5  mm  ins 

6.0 

6.0 

0.44 

1.61 

5  mm  cond 

8.7 

7.4 

5.68 

1.42 

8  mm  ins 

7.9 

8.2 

0.27 

1.59 

8  mm  cond 

10.6 

9.5 

6.00 

1.37 

measurements  were  acquired  on  a  rectangular  grid  of  64  x  64  for  the  whole  field  of  view 
and  approximately  1000-1 100  of  these  fell  into  the  object  region  in  the  mask.  Note  that  the 
number  of  measurements,  m ,  depends  on  the  mask  specific  to  each  case  defined  by  the  signal-to- 
noise  ratio.  The  highest  noise  in  magnetic  field  measurements  is  observed  at  the  phantom-air 
interface.  In  iterated  sensitivity  matrix  reconstruction,  it  is  possible  to  reconstruct  conductivity 
values  at  points  even  if  the  magnetic  flux  density  at  that  specific  location  is  not  available.  In 
order  to  avoid  the  errors  that  can  be  introduced  due  to  high  noise,  measurements  near  the 
boundary  are  masked  out  and  are  not  included  in  reconstruction.  The  mask  is  determined  by 
the  noise  level  in  the  MR  images  and  can  be  different  for  each  case.  For  most  cases,  out  of 
1576  pixels  of  the  magnetic  field  images,  1000-1 100  pixels  are  used  in  image  reconstruction. 
The  sensitivity  matrices  for  each  current  profile  were  stacked  to  have  the  combined  sensitivity 
matrix,  which  is  indicated  by  S  in  equation  (8),  of  size  2m  x  2048.  For  the  first  iteration  the 
maximum  singular  value,  crmax,  and  condition  number,  crmax/crmin,  of  STS  are  1.13  x  10  12  and 
5.2  x  1019,  respectively,  when  the  units  of  distance,  conductivity  and  current  amplitude  are 
cm,  S  cm-1  and  mA,  respectively,  with  an  initial  conductivity  distribution  of  0.016  S  cm-1. 
For  a  typical  case,  the  Tikhonov  regularization  parameter,  A.,  is  in  the  order  of  10“ 14  for  the 
first  iteration  which  is  in  the  order  of  1  / 100  of  the  maximum  singular  value.  These  values  were 
not  fixed  for  all  reconstructions  and  the  Tikhonov  regularization  parameter  was  recalculated 
for  each  sensitivity  matrix  using  the  criterion  given  in  equation  (9).  Above  numbers  were  only 
provided  to  give  a  reader  an  idea  about  the  order  of  these  parameters  for  a  typical  experiment 
in  this  study.  The  iterations  are  stopped  when  the  change  in  overall  conductivity  is  below  3% 
of  the  previous  iteration  and  this  corresponds  to  4-6  iterations. 

3.4.  Spatial  resolution 

In  order  to  understand  the  spatial  resolution  for  low  amplitude  currents,  phantoms  that  contain 
objects  with  different  diameters  (3  mm,  5  mm  and  8  mm)  and  extreme  conductivity  cases 
(insulator  (ins)  and  high  conductivity  (cond))  were  reconstructed.  For  the  insulator  cases, 
hollow  cylinders  filled  with  the  same  agarose-NaCl  solution  as  the  background  were  used. 
For  the  high  conductivity  cases,  a  solution  with  10  g  NaCl/100  g  water,  which  gives  a  6.70:1 
contrast  with  respect  to  background,  was  used. 

In  figure  4,  images  reconstructed  for  conductor  cases  are  given  together  with 
corresponding  profiles.  Profile  lines  are  as  defined  in  figure  3.  The  full- width  at  half¬ 
maximum  (FWHM)  values  were  calculated  in  two  orthogonal  profiles  passing  through  the 
centre  of  the  object  (FWHMv  for  the  v-direction  and  FWHMy  for  the  y -direction).  In 
figure  4(d)-(f),  profiles,  baseline,  and  half-maximum  levels  for  the  measurement  are  shown. 
FWHM  values  for  all  cases  are  summarized  in  table  2. 
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Figure  4.  Reconstructed  conductivity  images  for  3,  5  and  8  mm  conductor  objects  and 
corresponding  profiles. 


From  these  values,  it  is  seen  that  for  all  cases,  the  spread  is  more  for  the  conductor 
case  than  the  insulator  case.  Ideally,  the  spread  functions  are  to  be  obtained  using  extreme 
conductivity  cases  using  a  perfect  conductor  (oo  conductivity)  and  a  perfect  insulator 
(zero  conductivity).  The  experimental  insulator  case  matches  this  condition  closest  but  the 
conductor  case  phantom  uses  finite  conductivity  values.  Since  the  NaCl  concentrations  in  the 
object  and  background  are  different,  the  conductivity  of  the  object  decreases  with  time  due 
to  diffusion  (Hamamura  et  al  2006).  Hamamura  et  al  (2006)  carried  out  extensive  analysis 
on  diffusion  effects  in  agarose  phantoms  over  a  12  h  period  by  reconstructing  images  at  11 
different  time  points.  After  1  h,  the  FWHM  value  of  a  12  mm  diameter  object  increases  from 
1 1.98  mm  to  12.44  mm.  This  may  explain  the  difference  between  the  spread  functions  for  ins 
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Table  3.  Contrast  resolution  analysis. 


NaCl 

(g) 

True 

ratio 

Reconstructed 

ratio 

True 

peak 

Reconstructed 

peak 

True 

background 

Reconstructed 

background 

0.25 

0.36 

0.44 

0.43 

0.70 

1.61 

1.60 

1.2 

1.17 

1.32 

1.88 

2.09 

1.61 

1.58 

2 

1.88 

1.70 

3.02 

2.61 

1.61 

1.54 

4 

3.53 

3.45 

5.68 

4.71 

1.61 

1.36 

7 

5.64 

5.03 

9.08 

6.00 

1.61 

1.19 

and  cond  cases.  For  example,  the  8  mm  ins  case  is  equal  to  the  object  size  whereas  for  the 
5  mm  ins  case,  it  is  20%  more  than  the  object  size. 

Another  difference  between  the  conductor  and  insulator  cases  is  the  background  noise 
in  the  reconstructed  images.  For  all  cases,  conductor  images  are  less  noisy  compared  to  the 
insulator  cases  of  the  same  size  object  case.  The  main  reason  for  this  difference  is  the  loss  in 
the  measured  signal  amplitude  in  perfect  insulator  regions.  The  highest  change  in  magnetic 
flux  density  occurs  at  the  object  boundary;  however,  due  to  the  insulator  shell,  this  cannot 
be  accurately  measured  using  MRI.  Especially  for  the  3  mm  case,  where  the  total  volume  of 
the  insulator  shell  becomes  significant  compared  to  the  object  volume,  the  signal  level  drops 
significantly,  resulting  in  low  SNR.  Therefore,  this  limit  may  partially  be  imposed  by  the 
phantom  preparation. 

We  also  looked  at  the  peak  reconstructed  values  (maximum/minimum  conductivity  for 
conductor/insulator  cases)  in  the  object  and  mean  conductivity  in  the  background  and  these 
values  are  also  listed  in  table  2.  Note  that  the  significant  measure  for  spatial  resolution 
analysis  is  the  FWHM  value;  however,  the  numerical  conductivity  values  are  also  listed  to 
provide  readers  an  idea.  For  all  cases,  expected  background  conductivity  is  1.61  S  cm-1.  For 
the  insulator  cases,  the  expected  object  conductivity  is  zero.  For  the  high  conductor  case, 
although  the  true  conductivity  value  is  11.22  S  m-1,  reconstructed  values  are  much  lower, 
which  could  be  expected  due  to  diffusion  and  small  object  size  (PSF  effect).  If  the  peak  values 
for  different  cases  are  compared,  it  is  seen  that  as  the  object  size  increases,  the  value  increases, 
towards  the  true  conductivity  value. 

3.5.  Contrast  resolution 

In  order  to  understand  the  system  performance  at  different  contrast  levels,  15  mm  diameter 
objects  with  five  different  conductivity  values  between  the  two  extreme  cases  used  in  the  spatial 
resolution  analysis  were  selected.  The  true  conductivity  and  ratio  values  and  reconstructed 
peak  object  and  mean  background  values  are  summarized  in  table  3.  In  the  first  case,  the 
object  conductivity  is  set  lower  than  the  background  conductivity  (0.36: 1  contrast  case,  case  I). 
In  the  top  row  of  figure  5,  reconstructed  image  and  one  of  the  profiles  is  plotted.  Due  to  high 
errors  at  the  boundary,  the  object  is  not  clearly  visible  in  the  image;  however,  it  is  well 
resolved  in  the  profile  image.  Next,  we  used  a  low  contrast  conductor  object  where  the  object 
conductivity  was  17%  more  than  the  background  (case  II).  For  this  case,  the  noise  effects  are 
dominant  in  the  image  (figure  5,  second  row)  and  the  object  can  hardly  be  seen  in  the  profile 
images.  Then  we  increased  the  conductivity  further  using  three  contrast  levels,  1.88:1,  3.53:1 
and  5.64: 1  (cases  III,  IV  and  V,  respectively),  and  results  are  presented  in  figure  5  and  table  3. 
The  noise  in  the  background  becomes  less  significant  as  the  object  contrast  increases. 

For  the  conductor  cases  III,  IV  and  V,  the  peak  reconstructed  values  are  lower  than  the 
true  values  which  is  expected  due  to  the  diffusion  effect  explained  in  the  previous  section. 


MREIT  using  low  amplitude  current 


5045 


-i - * - * - * - 

0  20  40  60  80 

Position  (mm) 


5 


0  - ‘ - ‘ - 1 - 

0  20  40  60  80 

Position  (mm) 


Figure  5.  Reconstructed  images  and  corresponding  profiles  for  different  contrast  objects.  From 
top  to  bottom  (I)  0.36:1,  (II)  1.17:1,  (III)  1.88:1,  (IV)  3.53:1,  (V)  5.64:1  contrast  cases.  Values  are 
listed  in  table  3. 
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Figure  6.  Comparison  of  experimental  contrast  ratios  obtained  to  the  ideal  case. 


The  same  effect  results  in  a  higher  conductivity  peak  in  the  object  when  the  object  has  lower 
conductivity  (case  I)  and  the  direction  of  diffusion  is  reversed. 

In  figure  6,  the  graph  for  true  contrast  ratio  versus  reconstructed  conductivity  ratio  is 
presented.  In  this  plot,  the  dashed  line  corresponds  to  ideal  reconstruction.  It  is  seen  that  the 
deviation  from  the  theoretical  case  is  more  for  higher  conductivity  values,  which  could  be  due 
to  the  nonlinear  nature  of  the  governing  differential  equation. 


4.  Discussion 

In  this  study,  contrast  and  spatial  resolution  limits  in  MREIT  are  investigated  using  low 
level  ac  currents.  We  used  1  mA  peak  current  in  our  experiments  and  generated  conductivity 
images  with  better  resolution  and  accuracy  compared  to  other  conductivity  imaging  techniques. 
Current  levels  comparable  to  those  in  our  study  are  used  in  Muftuler’s  (Muftuler  et  al  2004) 
preliminary  study  where  a  couple  of  phantom  results  and  a  small  animal  result  were  presented. 
Since  the  pulse  sequence  and  reconstruction  algorithm  are  improved,  noise  performance  is 
also  improved. 

Small  objects  with  extreme  conductivity  values  were  prepared  for  understanding  the 
spatial  resolution  limits.  3  mm,  5  mm  and  8  mm  insulator  and  high  conductor  objects  were 
used.  Our  results  suggest  that  the  spatial  resolution  is  3  mm  or  better.  Due  to  diffusion  effects, 
for  high  conductivity  cases,  the  object  starts  spreading  in  the  conductivity  map  and  an  accurate 
spread  measurement  is  not  possible.  The  FWHM  of  the  3  mm  ins  object  was  measured  to  be 
4.2  mm.  Therefore,  we  can  state  that  the  resolution  is  around  3  mm  or  better  in  a  7  cm  object 
when  1  mA  current  is  used.  MRI  resolution  for  magnetic  flux  density  measurements  was 
64  x  64  for  a  field  of  view  (FOV)  of  100  mm.  For  this  image  matrix  and  FOV,  the  pixel  size 
in  magnetic  field  images  is  1.56  mm  x  1.56  mm  and  therefore,  this  is  the  theoretical  spatial 
resolution  that  can  be  achieved  with  this  set  of  imaging  parameters.  In  order  to  understand 
the  exact  resolution  limit,  we  tried  using  smaller  insulator  objects,  however,  since  we  cannot 
get  enough  MRI  signal  from  inside  such  a  small  insulator  object,  it  was  not  possible  to  take 
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measurements  from  such  phantoms.  For  in  vivo  applications,  there  may  be  air  spaces  such 
as  the  air  in  the  lung  and  in  the  intestines  with  low  conductivity  regions  which  cause  similar 
problems.  Recently,  Lee  et  al  (2006)  proposed  a  method  for  conductivity  reconstruction  from 
defective  data  in  MREIT  to  handle  such  situations. 

Several  phantoms  with  different  conductivity  values  were  prepared  to  understand  the 
contrast  resolution  limits  of  the  low  current  MREIT.  Conductivity  values  ranging  from 
approximately  one  fourth  of  the  background  to  six  times  the  background  were  used.  For  the  low 
conductivity  case,  although  the  object  was  clearly  resolved,  the  background  error,  especially  at 
the  boundary,  was  very  high.  For  the  low  contrast  conductor  object,  where  conductivity  of  the 
object  was  17%  higher  than  the  background,  a  similar  result  was  obtained  and  the  object  was 
not  resolved.  It  has  been  reported  by  several  groups  that  the  highest  error  in  measurements 
is  observed  at  the  boundary.  In  a  recent  study,  Hamamura  (2005)  suggested  several  masking 
and  correction  algorithms  to  eliminate  these  problems.  For  higher  conductivity  contrast 
cases,  the  boundary  effects  are  not  as  significant.  Diffusion  of  NaCl  from  higher  to  lower 
concentration  also  affects  the  reconstruction  here.  For  low  conductivity  cases,  diffusion  is 
from  the  background  towards  the  object  and  reconstructed  object  conductivity  is  higher.  For 
higher  conductivity  cases,  the  direction  of  diffusion  is  reversed  and  the  reconstructed  values 
are  lower  than  the  expected  values  as  can  be  seen  from  figure  6. 

In  this  study,  in  order  to  reduce  total  data  acquisition  time,  we  only  used  two  current 
injection  profiles,  which  is  the  minimum  required  for  uniqueness  and  was  sufficient.  The 
total  data  acquisition  time  was  8  min  for  this  case.  Depending  on  the  application,  more 
current  injection  profiles  may  be  required.  If  the  data  acquisition  time  is  limited,  trade 
off  between  current  profiles  and  number  of  the  averages  must  be  investigated  for  optimum 
SNR. 

The  first  step  in  the  image  reconstruction  is  the  selection  of  the  initial  conductivity 
distribution.  Theoretically,  as  long  as  a  uniform  conductivity  distribution  is  chosen,  the 
magnetic  flux  density  is  independent  of  its  value.  When  we  did  reconstruction  using  real  data, 
however,  we  observed  that  initial  conductivity  distribution  affects  the  convergence  speed.  For 
the  cases  presented  here,  the  results  converged  in  4-6  iterations  if  the  initial  conductivity 
distribution  was  equal  to  or  lower  than  the  background  and  a  higher  number  of  iterations 
was  required  if  the  initial  conductivity  distribution  was  higher  (approximately  ten-fold)  than 
the  background.  Also  note  that  reconstruction  results  provide  relative  conductivity  values, 
which  is  a  scaled  version  of  the  true  values,  due  to  the  fact  that  only  the  magnetic  flux 
density  measurements  are  used.  For  breast  cancer,  the  relative  conductivity  value  of  tumour 
with  respect  to  background  or  a  reference  organ  would  provide  the  necessary  information. 
However,  if  the  absolute  conductivity  values  are  required,  the  scaling  factor  can  be  determined 
using  a  single  voltage  measurement. 

Recently,  several  researchers  have  studied  the  use  of  MRI  in  breast  cancer  research.  MRI 
is  highly  sensitive  and  can  detect  the  location  and  sizes  of  the  tumours  precisely;  however, 
its  specificity  is  low  and  variable  (Su  et  al  1998).  Our  main  motivation  in  proposing  the 
use  of  MREIT  in  breast  cancer  research  is  its  potential  to  increase  the  specificity  using  the 
conductivity  information.  It  is  reported  that  the  conductivity  of  malignant  tissues  maybe  20- 
40  times  higher  than  benign  tissues  in  some  forms  of  cancer  (Surowiec  et  al  1988).  Therefore, 
we  expect  both  the  spatial  and  the  contrast  resolution  achieved  in  this  study  are  sufficient  for 
applicability  of  the  technique  in  breast  cancer.  The  main  limitation  in  clinical  application  of 
this  technique  is  the  safe  current  amplitude  for  patient  safety.  In  the  future  work,  a  priori 
anatomical  MRI  information,  which  is  available  for  any  MREIT  experiment,  can  be  included 
in  MREIT  reconstruction  for  the  improvement  of  the  accuracy  of  the  conductivity  with  lower 
current  levels  towards  the  aim  of  translation  of  the  technique  to  clinical  applications. 
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5.  Conclusions 

Results  presented  here  are  critical  in  understanding  the  limits  of  MREIT  at  lower  currents 
which  is  required  for  translation  of  the  study  to  animal  and  ultimately  to  human  applications. 
Comparison  of  our  results  with  the  literature  suggests  that  using  multiple  current  pulses 
provides  better  signal  to  noise  in  measurements,  and  sensitivity-based  reconstruction  is  less 
susceptible  to  noise  in  measurements  compared  to  the  gradient-based  algorithms. 
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A  significant  increase  in  electrical  conductivity  of  neoplasticities  compared  to  healthy  tissues 
and  benign  formations  has  been  reported  in  several  studies.  We  previously  reported  prelimi¬ 
nary  results  with  MR  based  Electrical  Impedance  Tomography  (MREIT)  on  several  phantoms 
and  a  single  animal.  In  the  presented  study,  we  applied  the  technique  on  ten  tumor-bearing 
rats  and  collected  MREIT  images  to  investigate  the  potential  of  MREIT  for  characterizing 
malignant  tumors.  Results  show  that  the  tumors  had  significantly  higher  mean  conductivity 
compared  to  the  mean  of  conductivity  in  the  rest  of  the  body.  Although  heterogeneity  of  con¬ 
ductivity  was  observed  in  the  tumor,  the  mean  was  still  higher  than  the  background. 

Keywords:  MREIT;  Electrical  Impedance  Tomography;  MRI;  and  Cancer  imaging. 

Introduction 

Several  in  vitro  studies  have  shown  that  the  electrical  impedance  of  malignant  tis¬ 
sues  is  significantly  higher  than  those  of  normal  and  benign  tissues  (1,2).  There¬ 
fore,  in  vivo  impedance  imaging  of  suspicious  lesions  has  the  potential  to  improve 
the  sensitivity  and  specificity  of  detecting  malignant  tumors.  Imaging  of  conduc¬ 
tivity  distribution  can  be  achieved  by  Electrical  Impedance  Tomography  (EIT)  in 
which  weak  electrical  currents  are  injected  into  or  induced  in  the  body  and  surface 
voltage  measurements  are  made  to  reconstruct  a  conductivity  image  (3).  How¬ 
ever,  in  this  technique  a  limited  number  of  measurement  electrodes  are  placed  on 
the  outer  surface.  This  results  in  poor,  non-uniform  spatial  resolution  due  to  low 
sensitivity  of  the  voltages  to  conductivity  changes  in  the  interior  regions  (4). 

Imaging  conductivity  distribution  for  diagnosing  cancer  has  been  proposed  by 
several  researchers  (5,6).  Cherepenin  et  al.  proposed  a  conventional  EIT  system 
for  breast  cancer  screening,  while  Malich  et  al.  used  a  newly  developed  instru¬ 
ment  called  Electrical  Impedance  Scanning  (EIS).  In  EIS,  electrical  currents  are 
injected  to  the  body  by  a  single  probe  held  by  the  patient  and  an  electrode  array 
placed  over  the  breast  is  used  to  measure  the  distribution  of  currents  on  the  sur¬ 
face  of  the  breast.  This  information  is  used  to  estimate  the  distribution  of  conduc¬ 
tivity  under  the  electode  array.  Although  these  techniques  proved  to  be  useful, 
both  suffer  from  poor  spatial  resolution  and  poor  sensitivity  to  deeper  structures. 

Magnetic  Resonance-Electrical  Impedance  Tomography  (MREIT)  has  been  re¬ 
cently  introduced,  in  which  weak  electrical  currents  are  injected  into  the  tissue 


Abbreviations:  MRI,  Magnetic  Resonance  Imaging;  EIT,  Electrical  Impedance  Tomography; 
MREIT,  Magnetic  Resonance  Electrical  Impedance  Tomography;  FEM,  Finite  Element  Method; 
EIS,  Electrical  Impedance  Scanning;  OPAMP,  Operational  Amplifier;  SVD,  Singular  Value  De¬ 
composition;  NEX,  Number  of  Excitations;  SMM,  Sensitivity  Matrix  Method 
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and  the  resulting  perturbations  in  the  magnetic  flux  density 
are  measured  using  phase  information  in  MR  images.  The 
conductivity  distribution  inside  the  body  is  reconstructed 
from  these  magnetic  flux  density  measurements  (7-9).  Un¬ 
like  conventional  EIT,  the  MREIT  technique  provides  mea¬ 
surements  from  inside  the  object  on  a  uniform,  high  resolu¬ 
tion  grid.  This  provides  uniformly  distributed  high  spatial 
resolution  conductivity  images. 

Several  MREIT  studies  using  biological  tissue  phantoms 
have  been  previously  reported  (10-12).  However,  these  stud¬ 
ies  were  done  on  in  vitro  samples  using  injected  current  lev¬ 
els  of  12  mA  and  above.  For  any  in  vivo  study,  these  current 
levels  are  substantially  higher  than  biologically  safe  current 
limits  and  will  certainly  cause  involuntary  muscle  twitching 
even  if  they  do  not  harm  the  animal.  In  some  of  our  experi¬ 
ments,  we  have  observed  involuntary  muscle  twitching  even 
at  current  levels  as  low  as  2mA,  if  the  electrodes  are  placed 
near  a  major  nerve  bundle  by  coincidence.  In  those  cases  the 
experiment  is  immediately  stopped  and  electrodes  are  relo¬ 
cated  to  avoid  any  discomfort  to  the  animal. 

We  have  recently  reported  our  preliminary  studies  with  phan¬ 
toms  as  well  as  one  in  vivo  experiment  (13).  Here,  we  present 
the  results  of  MREIT  performed  on  ten  animals.  Parameters 
like  variance  and  mean  in  the  tumor  versus  the  rest  of  the 
body  were  investigated.  The  goal  is  to  verify  the  potential  of 
MREIT  to  aid  in  the  diagnosis  of  tumors. 

Methods 

Data  were  collected  using  a  4T  whole  body  MRI  magnet 
(Magnex  Scientific  Inc.,  UK),  which  is  interfaced  with  a 
MRRS  console  (Magnetic  Resonance  Research  Systems, 
Guildford,  UK)  that  has  broadband  RF  transmit  and  receive 
channels.  The  system  is  equipped  with  a  whole  body  gradi¬ 
ent  coil  set  (Tesla  Eng.  UK),  which  provides  up  to  3G/cm 
gradient  fields.  The  clear  bore  of  the  magnet  is  650mm  in 
diameter  with  the  gradient  tube.  The  system  also  includes  a 
13  channel  room-temperature  high-order  shim  system  with 
MXA-13-4  shim  power  supply  (Resonance  Research,  Bil¬ 
lerica,  MA)  to  minimize  field  inhomogeneities.  A  16  leg, 
quadrature,  high-pass  birdcage  coil  was  designed  and  built 
in-house  for  the  MREIT  experiments.  This  coil  has  10  cm 
diameter  and  1 8  cm  length,  which  is  suitable  to  image  large 
size  sprague-dawley  rats  used  in  the  experiments. 

Pulse  Sequence 

The  pulse  sequence  used  for  the  MREIT  experiments  em¬ 
ploys  a  train  of  180°  RF  pulses  following  a  90°  RF  pulse  (Fig. 
1).  No  spatial  encoding  gradients  were  applied  between  the 
180°  RF  pulses  and  the  data  was  collected  with  a  single  read¬ 
out  gradient  only  after  the  last  180°  RF  pulse.  Alternating 


electrical  currents  were  injected  to  the  animals  in  the  form  of 
burst  sine  wave  pulses  in  synchrony  with  the  RF  pulse  train 
where  polarity  of  current  changed  after  each  180°  RF  pulse. 
The  duration  of  each  180°  RF  pulse  was  1.3ms.  A  similar 
sequence  was  proposed  by  Mikac  et  al.  to  obtain  an  image  of 
electric  currents  inside  an  object  (14). 


Figure  1:  The  MREIT  pulse  sequence  for  multi-slice  conductivity  imaging. 

These  currents  flowing  inside  the  animal’s  body  generate  a 
magnetic  flux  density.  The  component  of  the  magnetic  flux 
density  that  is  parallel  to  the  main  static  field  (z-component) 
introduces  a  phase  shift  that  accumulates  over  each  half-cycle 
of  the  sine  wave.  The  signal  equation  for  MREIT  is  given  as: 

s(li,v)  =  ff  M(x,y)  e^xu  +y  v>  eMo  (b(x,y)cos(wt))dt  dxdy 

[1] 

Here,  M(x,y )  is  the  magnetization  of  protons  at  (x,y),  and  the 
Tl,  T2  decays  are  ignored  to  simplify  the  equation  (they  can 
be  considered  as  constant  weighting  factors,  thus  M(x,y )  ac¬ 
counts  for  the  proton  density  as  well  as  the  Tl  and  T2  weight¬ 
ing).  The  constant  phase  shifts  due  to  various  sources  such 
as  static  field  inhomogeneity  and  other  hardware  related 
phase  delays  are  all  summarized  in  the  phase  term  0(JCj);  y  is 
the  gyromagnetic  ratio,  and  u  and  v  are  the  spatial  frequen¬ 
cies  given  by  u  -  y*Gx*TGx  and  v  =  y*Gy*TGy,  where  TGx 
and  TGy  are  the  durations  of  the  gradient  pulses  Gx  and  Gy, 
respectively.  Therefore,  in  the  MREIT  image  the  total  phase 
accumulated  at  a  pixel  at  location  (x,y)  due  to  the  magnetic 
flux  density  generated  by  injected  currents  is: 

cp(x,y)  =  4  •  y  •  N  •  b(x,y)  /  oo  [2] 

where  N  is  the  number  of  cycles  of  injected  current,  oa  is 
the  angular  frequency,  and  b(x,y)  is  the  amplitude  of  the 
magnetic  flux  density  at  point  (x,y)  (z-component  only) 
generated  by  injected  currents.  This  phase  term  in  Equa- 
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tion  [2]  can  be  easily  measured  from  the  MREIT  images. 
Therefore,  the  z-component  of  the  magnetic  flux  density 
distribution  can  be  calculated. 

Since  the  constant  phase  term  0(jcj)  is  unknown,  the  data  has 
to  be  collected  twice  with  the  polarity  of  electric  current  re¬ 
versed.  When  the  phase  of  these  two  MREIT  images  are  sub¬ 
tracted,  static  phase  terms  will  be  eliminated,  leaving  only 
the  phase  term  introduced  by  b(x,y). 

MREIT  Hardware:  Several  hardware  components  were  in¬ 
terfaced  to  the  MRI  system  to  perform  the  experiments.  The 
MRI  console  controls  the  whole  experiment,  which  gener¬ 
ates  the  pulse  sequence,  acquires  the  incoming  data,  and  syn¬ 
chronizes  the  external  units.  The  sine  waves  were  produced 
by  an  HP  ESG-4400B  signal  generator.  These  signals  were 
synchronized  to  the  pulse  sequence  by  a  TTL  pulse  generated 
by  the  scanner  computer.  A  transconductance  amplifier  was 
designed  and  built  using  three  LM741  OPAMP  circuits  to 
convert  the  voltage  from  the  signal  generator  into  a  current 
output.  This  circuit  is  described  in  detail  in  (15),  which  was 
replicated  from  (16).  The  output  current  was  calibrated  to 
lma/2V  and  an  independent  measurement  of  current  was  not 
used.  This  experimental  procedure  was  described  in  more 
detail  in  (13).  After  data  collection,  raw  MRI  data  (k-space) 
was  uploaded  to  another  computer  for  off-line  processing. 

MREIT  Conductivity  Image  Reconstruction 

A  relationship  between  the  conductivity  distribution  inside 
an  object  and  the  measured  magnetic  flux  density  has  to  be 
established  to  be  able  to  reconstruct  images  of  conductiv¬ 
ity.  If  all  three  orthogonal  components  of  the  magnetic  field 
were  known,  Maxwell’s  equations  could  have  been  utilized 
to  calculate  the  current  density  distribution.  However,  this 
is  only  possible  by  rotating  the  object  in  three  orthogonal 
directions  inside  the  MRI  system  which  will  not  be  practi¬ 
cal  with  human  subjects.  Even  with  small  animals,  rotation 
would  introduce  problems  with  registration  of  measure¬ 
ments  taken  with  different  orientations  of  the  body  under 
investigation.  The  geometry  of  the  soft  tissues  would  de¬ 
form  due  to  gravity,  complicating  the  registration  further. 
Moreover,  a  larger  RF  coil  would  have  to  be  used,  which 
would  reduce  the  SNR.  Therefore,  a  method  was  adopted 
that  uses  only  the  z-component  of  the  magnetic  flux  density 
to  calculate  the  conductivity  distribution  (7,  8). 

In  order  to  calculate  the  relationship  between  the  conduc¬ 
tivity  distribution  of  an  object  and  the  magnetic  flux  density 
generated  by  currents  flowing  inside  that  object,  we  first 
start  with  the  calculation  of  the  electric  potential  distribu¬ 
tion  c|)(x,y,z)  in  the  imaging  slice.  This  is  calculated  by 
solving  Poisson’s  equation  with  Neumann  boundary  condi¬ 
tions,  which  is  given  by: 


V  •  (aV(|))(v,y)  =  0  (x,y)  E  D 


a 


d§_ 

dn 


J  on  positive  current  electrode 
-J  on  negative  current  electrode 
0  elsewhere 


[3] 


where  a  is  the  electrical  conductivity,  §  is  the  electric  po¬ 
tential,  and  D  is  the  slice  of  object  to  be  imaged.  The  above 
problem  is  a  nonlinear  equation  since  the  electric  potential 
itself  is  also  conductivity  dependent.  It  is  difficult  to  obtain 
an  analytic  solution  for  Equation  [3];  therefore,  numerical 
methods  should  be  used.  In  this  study,  the  Finite  Element 
Method  (FEM)  was  used  to  calculate  the  distribution  of  the 
electric  potential.  Once  §  is  calculated,  the  electric  field  and 
current  density  distribution  inside  the  imaging  region  may  be 
found  using  the  following  equations: 


£  =  -V( () 
J=oE 


[4] 


Finally,  the  magnetic  flux  density  generated  by  this  ohmic 
current  can  be  calculated  using  the  Biot-Savart  law.  The 
magnetic  field  can  be  written  in  terms  of  the  differential  cur¬ 
rent  element  Idl  as, 


[5] 

where  \i0  is  the  permeability  constant  and  R  is  the  vector  from 
the  source  point  at  {x\  y  ’,  z’)  to  the  field  point  (x,  y,  z).  Using 
Equations  [3]-[5],  one  can  solve  the  forward  problem  where 
the  magnetic  field  is  calculated  for  an  initial  conductivity  dis¬ 
tribution.  In  order  to  solve  the  inverse  problem  of  finding  the 
conductivity  image  from  the  magnetic  flux  density  measure¬ 
ments,  the  sensitivity  matrix  method  (SMM)  was  employed, 
in  which  a  linear  relationship  between  the  conductivity  per¬ 
turbations,  A  a,  and  magnetic  flux  density  perturbations,  AB, 
was  assumed  (8).  This  relationship  is  given  as: 


ABl  ' 
A  B2 

dBl  dB{ 
dCT^  d02 

dZ?2  dZ?2 

3(3-^  802 

dB, 

dB2 

do; 

Aox 

Ao2 

.ABm. 

dBm  dBm 

..  dBm 

Ao3 

ABZ  =  S  •  A  a 

Here,  AB  =  Bmeas  -  Bj  and  Act  =  acalc  -  c^;  where  is  the 
initial  estimate  of  the  conductivity  distribution  and  Bj  is  the 
field  distribution  corresponding  to  For  most  cases,  a  uni¬ 
form  distribution  is  assumed  for  Bmeas  is  the  magnetic 
flux  density  measurements  obtained  from  MRI.  If  there  are 
n  conductivity  elements  in  the  FEM  and  m  measurement 
points  from  MREIT  data,  then  A  a  is  an  nx  1  vector,  AB  is 
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an  mx  1  vector  and  the  sensitivity  matrix  S  is  an  m  x  n  ma¬ 
trix.  The  entry  of  the  sensitivity  matrix  at  ith  row  and  jth 
column  denotes  the  change  in  the  ith  field  measurement  due 
to  a  small  change  in  conductivity  of  the  jth  element.  The  sen¬ 
sitivity  matrix  can  be  calculated  using  either  a  numerical  or 
a  semi-analytical  approach.  In  the  numerical  approach,  each 
column  of  the  sensitivity  matrix  is  calculated  separately  by 
solving  the  forward  problem  by  changing  the  conductivity 
of  a  single  element.  This  approach  requires  the  repetitive 
solution  of  the  forward  problem  and  increases  the  execution 
time.  For  faster  reconstruction,  the  semi-analytical  method 
was  used  for  the  calculation  of  the  sensitivity  matrix  which 
was  described  in  (8).  Once  S  is  calculated,  the  conductivity 
distribution  can  be  approximated  by: 

Ag  =  S1  •  AB  [7] 

Since  SMM  approximates  the  relationship  between  perturba¬ 
tions  in  conductivity  and  perturbations  in  Bz  (z  component 
of  magnetic  flux  density)  with  a  linear  matrix  equation,  this 
reconstruction  method  provides  accurate  results  only  for 
small  conductivity  perturbations.  On  the  other  hand,  several 
researchers  reported  substantial  changes  in  conductivity  be¬ 
tween  normal  and  malignant  tissues  (1,2).  Therefore,  SMM 
will  tend  to  underestimate  the  actual  conductivity  contrast. 
In  this  study,  an  iterated  SMM  with  Tikhonov  regularization 
was  used  to  solve  the  nonlinear  imaging  problem  and  also  to 
reduce  the  artifacts,  especially  near  the  boundaries  close  to 
the  electrodes.  The  solution  to  Equation  [7]  is  found  by  solv¬ 
ing  the  least  squares  problem: 

min  ||  SAo  -  Ab  ||  [8] 

where  ||  .  ||  is  the  L2  norm.  For  most  cases,  the  problem 
above  is  ill-conditioned  and  a  regularization  method  that 
computes  an  approximate  solution  through  a  regularization 
parameter  is  required.  In  this  study,  the  linear  equation  was 
solved  using  the  conjugate  gradient  method  with  Tikhonov 
regularization,  where  Gcalc  was  found  by  introducing  an  ad¬ 
ditional  term  to  the  cost  function: 

min  {||  SAg  -  Ab||2  +  X  ||Ag||2}  [9] 

where  X  is  the  regularization  parameter.  This  approach  at¬ 
tempts  to  minimize  the  residual  norm  while  penalizing  large 
perturbations  in  the  solution  (Ag).  The  regularization  param¬ 
eter  X  was  chosen  such  that  the  calculated  Bz  generated  by 
the  reconstructed  conductivity  distribution  was  closest  to  the 
magnetic  flux  density  measured  by  MREIT: 

mm{|| Bfinai  -  [B(a)]0  =  ocafc(>J}  [10] 

Once  ocalc  is  calculated,  the  change  in  conductivity  Ag  was 
checked  against  a  predefined  threshold,  then  this  ocalc  was  as¬ 


signed  as  the  new,  updated  oi9  and  the  process  was  iterated. 

The  SMM  reconstructs  relative  conductivities.  Thus,  the  true 
conductivity  distribution  is  related  to  the  reconstructed  con¬ 
ductivity  distribution  by  a  constant  scaling  factor.  Determin¬ 
ing  the  absolute  conductivity  distribution  requires  an  addition¬ 
al  constraint,  such  as  a  surface  voltage  measurement  (17). 

In  Vivo  Experiments 

Ten  rats  bearing  malignant  tumors  were  imaged  in  this  study. 
Data  sets  from  two  animals  were  discarded  due  to  severe 
motion  artifacts.  The  tumors  were  either  R3230AC  tumor 
grafts  or  induced  by  the  carcinogen  ENU  (N-ethyl-N-nitro- 
sourea)  (18).  Animals  were  anesthetized  by  IV  injection 
of  ketamine  and  xylazine  prior  to  imaging.  All  procedures 
were  approved  by  the  IACUC. 

An  animal  holder  was  prepared  from  acrylic  sheets  for  in 
vivo  MREIT  experiments  (Fig.  2).  This  holder  helped  keep 
the  animal  and  the  electrodes  stationary  during  imaging. 
This  structure  also  enabled  placement  of  the  electrodes  in 
consistent  positions  across  different  animals  as  well  as  on  the 
same  animals  for  potential  longitudinal  studies.  5mm  x  5mm 
self-adhesive  copper  foil  was  taped  over  hollow  acrylic  tubes 
as  electrodes.  Those  tubes  were  filled  with  CuS04  solution 
as  an  MRI  marker  to  detect  the  electrode  positions  precisely 
in  the  images.  Precise  localization  of  electrode  positions  is 
essential  to  enter  correct  boundary  conditions  for  accurate 
reconstruction  of  impedance  images.  Otherwise,  artifacts 
were  seen  especially  along  the  periphery,  near  the  electrodes. 
Thin  copper  wires  that  carried  the  electrical  current  ran  along 
these  tubes,  and  the  tubes  were  aligned  in  z-direction  to  mini¬ 
mize  interference  from  the  magnetic  flux  density  generated 
by  current  in  the  wires.  The  electrodes  were  covered  with  a 
thin  layer  of  conductive  gel  to  provide  good  electrical  con¬ 
tact.  The  skin  areas  where  the  electrodes  would  be  placed 
were  shaved  for  better  conductance. 


Figure  2:  Animal  holder  with  a  rat  is  shown.  The  electrodes  are  placed  on 
the  tubes  and  current  carrying  wires  run  along  those  tubes. 


Structural  images  were  collected  using  a  T2  weighted  SE 
sequence.  Scan  parameters  were:  data  matrix  =  128  x  128; 
FOV  =10  cm;  slice  thickness  =  4  mm,  with  2  mm  gap.  TR 
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=  3  s,  TE  =  50  ms,  and  NEX  =  2  were  used.  MREIT  images 
were  collected  using  the  outlined  pulse  sequence  with  TR  = 
1000  ms,  TE  =  30  ms,  and  NEX  =  2;  64  x  64  data  matrix, 
FOV  =  10  cm,  slice  thickness  =  4  mm  with  2  mm  gap.  Two 
cycles  of  100Hz  current  with  1mA  peak  was  applied  for  all 
studies.  One  set  of  MREIT  data  was  collected  by  applying 
the  current  between  one  pair  electrodes.  Since  there  were 
four  electrodes,  subsequent  MREIT  data  were  collected  by 
applying  current  from  different  pairs,  generating  six  dif¬ 
ferent  current  profiles  (Fig.  3).  As  mentioned  previously, 
data  were  collected  with  both  ±  polarities  of  the  currents  to 
eliminate  phase  accumulation  from  other  sources.  In  two 
of  the  rats,  we  have  also  collected  MRIET  data  with  current 
amplitudes  of  0.5,  1,  and  2mA. 


observed  in  T2  weighted  MRI  images.  We  have  also  calcu¬ 
lated  the  ratio  of  standard  deviation  to  mean  conductivity  in 
each  tumor  region,  which  may  be  an  indication  of  conductiv¬ 
ity  heterogeneity  inside  the  tumor  volume,  rather  than  SNR 
(Figure  6).  As  seen  from  these  figures,  consistent  results  were 
obtained  from  these  eight  animals.  Average  conductivity  in¬ 
creased  by  roughly  2.2  times  in  the  tumor  compared  to  the 
rest  of  the  body.  The  conductivity  varied  typically  between 
10%  and  20%  within  the  tumor.  In  two  of  the  rats,  from 
which  we  have  collected  MRIET  data  with  current  ampli¬ 
tudes  of  0.5,  1,  and  2mA,  highly  consistent  results  were  seen 
especially  between  1  and  2mA  cases.  For  example,  when 
MREIT  images  of  1mA  and  2mA  cases  were  subtracted,  the 
mean  of  the  residual  was  only  1.3%  of  the  mean  conductivity 
in  the  whole  slice.  In  the  case  of  0.5mA  vs.  1mA,  the  mean 
of  the  residual  was  12%.  Since  the  SMM  gives  relative  con¬ 
ductivity  distributions,  the  reconstructed  conductivity  values 
have  an  arbitrary  scaling;  therefore,  the  mean  of  each  MREIT 
image  is  scaled  to  the  same  value  before  subtraction. 


Figure  3:  Six  different  current  injection  profiles  with  four  elec¬ 
trodes  used  in  the  experiments. 

Magnetic  flux  density  calculations  and  conduc¬ 
tivity  reconstruction  algorithms  utilizing  the  it¬ 
erative  SMM  with  Tikhonov  regularization  were 
developed  using  MATLAB.  For  reconstruction, 
the  FEM  mesh  boundary  was  determined  using  a 
standard  MR  magnitude  image.  Between  n  =  151 
to  1037  conductivity  elements  and  m  =  280  to  473 
measurement  points  were  utilized  in  reconstruc¬ 
tion,  with  the  exact  values  being  specific  to  the  ge¬ 
ometry  of  the  animal  being  analyzed. 

Results 


Figure  4:  Structural  and  MREIT  images  of  six  animals  are  illustrated.  Anatomic  (T2 
weighted)  scans  are  displayed  in  gray  levels  and  corresponding  impedance  (MREIT)  im¬ 
ages  are  depicted  in  color  right  below  the  T2  weighted  images.  Each  image  pair  shows 
axial  images  from  different  animals.  Tumor  areas  are  circled  with  red  lines.  Bright  ob¬ 
jects  outside  the  animals’  body  were  markers  to  identify  exact  location  of  electrodes. 

Discussion 


Figure  4  shows  T2  weighted  MRI  and  MREIT  images  of  six 
rats.  Tumor  areas  show  increased  conductivity  depicted  with 
yellow-red  colors  (see  colorbar).  On  the  MREIT  images  of 
the  animals,  separate  ROIs  (region  of  interest)  were  drawn 
over  the  tumor  region  and  the  rest  of  the  body  and  the  mean 
conductivity  values  in  these  ROIs  were  calculated.  Since  the 
SMM  yields  relative  conductivity  values,  the  ratio  of  mean 
conductivities  otumor/abody  was  calculated  for  each  animal  and 
the  graph  is  shown  in  Figure  5.  It  was  found  that  the  aver¬ 
age  of  these  conductivity  ratios  pooled  over  eight  animals 
was  2.17.  ROIs  were  drawn  manually  based  on  the  tumor 


In  this  study,  it  has  been  demonstrated  that  MRI  based  im¬ 
pedance  imaging  has  the  potential  to  investigate  malignant 
tumors  in  vivo.  This  method  can  be  used  to  identify  and  char¬ 
acterize  neoplasticities,  which  are  known  to  possess  higher 
electrical  conductivity  with  respect  to  healthy  tissues  and  be¬ 
nign  formations.  In  using  the  SMM,  only  relative  conductiv¬ 
ities  can  be  found.  However,  only  contrasts  in  conductivity 
were  required  to  differentiate  tumors  from  the  surrounding 
background  tissue.  Higher  conductivity  inside  the  tumor  was 
observed  consistently  in  all  the  animals  imaged.  Moreover, 
it  was  observed  that  the  method  was  sensitive  enough  to  re¬ 
veal  conductivity  heterogeneity  in  tumors.  The  same  heter- 
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Figure  5:  The  ratio  of  mean  conductivity  in  the  tumor  versus  the  rest  of 
the  body  over  eight  in  vivo  studies.  Mean  and  standard  deviation  of  these 
measurements  are  indicated  by  the  dash  and  the  error  bar,  respectively. 


Figure  6:  The  ratio  of  standard  deviation  of  conductivity  versus  its  mean 
in  the  tumor  in  eight  in  vivo  studies.  Mean  and  standard  deviation  of  these 
measurement  are  indicated  by  the  dash  and  the  error  bar,  respectively. 


ogenous  structure  was  observed  at  different  current  levels  in 
the  same  animal;  therefore,  the  variations  were  not  random 
noise.  There  is  still  possibility  that  some  systematic  errors 
in  data  acquisition  or  reconstruction  might  have  caused  such 
consistent  variations.  Even  though  such  patterns  were  not 
observed  in  phantom  studies,  the  much  more  complex  struc¬ 
ture  of  the  animal  might  have  amplified  such  systematic  er¬ 
rors.  This  needs  to  be  investigated  further. 


cation  was  implemented  to  reduce  computation  time.  How¬ 
ever,  the  quality  and  accuracy  of  conductivity  images  can  be 
further  improved  by  utilizing  3D  FEM  based  reconstruction. 
Some  of  the  residual  noise  or  artifacts  in  the  conductivity 
images  could  be  due  to  the  superposition  of  magnetic  fields 
generated  by  currents  flowing  outside  the  imaging  slice.  This 
can  only  be  modeled  by  3D  FEM  and  reconstructed  accu¬ 
rately,  which  is  currently  being  developed. 


We  have  previously  carried  out  extensive  phantom  studies  to 
optimize  the  method  and  assess  the  spatial  resolution,  con¬ 
trast,  and  linearity  (13).  More  robust  reconstruction  with 
regularization  was  developed  to  minimize  artifacts.  Software 
correction  algorithms  were  developed  to  minimize  electrode 
localization  errors,  which  helped  reduce  boundary  artifacts. 

Although  the  animal  model  of  breast  cancer  used  in  this  study 
is  a  suitable  one  to  perform  preliminary  studies  of  in  vivo 
MREIT,  it  also  introduced  motion  artifact  problems.  This 
could  be  easily  avoided  in  actual  human  studies.  In  breast 
imaging,  the  breasts  are  stationary  and  chest  motion  can  be 
kept  outside  the  breast  images  by  selecting  phase  encoding  in 
the  left-right  direction.  However,  in  future  tests  with  animal 
models,  we  plan  to  utilize  motion  correction  schemes  such 
as  navigator  echoes  to  minimize  the  errors  caused  by  mo¬ 
tion.  Moreover,  the  abdominal  region  of  a  rat  is  a  more  com¬ 
plicated  structure  compared  to  a  woman’s  breast,  when  all 
the  internal  organs  and  their  inherent  motion  are  considered. 
Therefore,  we  expect  the  MREIT  images  to  be  more  accurate 
in  potential  applications  of  MREIT  on  human  breast. 

This  study  utilized  a  2D  FEM  based  reconstruction,  where  the 
electrodes  were  placed  within  the  imaging  plane  and  the  cur¬ 
rent  flow  assumed  restricted  within  this  plane.  While  formu¬ 
lation  of  the  SMM  remains  the  same  for  3D,  the  increase  in 
the  size  of  the  matrices  and  added  complexity  of  the  3D  FEM 
presents  additional  challenges.  For  this  study,  a  2D  simplih- 


The  positioning  of  the  surface  electrodes  also  plays  a  crucial 
role,  since  this  direct  affects  the  current  density  distribution. 
Conceivably,  placing  the  electrodes  in  a  certain  arrangement 
may  improve  the  reconstruction  results  in  specific  regions 
( e.g .,  regions  with  higher  current  density)  while  reducing  the 
sensitivity  in  others.  In  this  study,  the  electrodes  were  ar¬ 
ranged  in  a  generic  manner  by  placing  them  within  the  imag¬ 
ing  plane  and  approximately  equidistant  around  the  animal. 
The  reconstruction  may  be  improved  by  using  a  priori  in¬ 
formation  to  determine  specific  regions  of  interest,  and  then 
placing  the  electrodes  accordingly.  Such  optimation  of  elec¬ 
trode  positioning  needs  to  be  investigated  further. 

At  this  stage  the  technique  may  not  be  used  as  a  standalone 
diagnostic  tool;  but  it  could  provide  useful  information  in 
characterizing  the  tumor,  once  a  suspicious  lesion  is  detected 
by  other  methods.  Physiological  and  structural  changes  in 
tumors  that  lead  to  such  changes  in  conductivity  will  be  in¬ 
vestigated  in  future  studies.  Note  that  good  quality  MREIT 
images  were  collected  with  biologically  safe  electrical  cur¬ 
rent  levels.  To  the  best  of  our  knowledge,  our  research  group 
is  the  only  one  to  demonstrate  potential  use  of  MR  based 
impedance  imaging  in  diagnosing  or  characterizing  tumors. 
The  method  is  still  in  its  early  development  stages.  Further 
refinements  of  this  technique  could  improve  accuracy  more 
and  may  reveal  important  information  about  the  tumor  char¬ 
acteristic  that  cannot  be  observed  with  other  techniques. 
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Abstract 

In  magnetic  resonance  electrical  impedance  tomography  (MREIT),  currents  are 
applied  to  an  object,  the  resulting  magnetic  flux  density  measured  using  MRI 
and  the  conductivity  distribution  reconstructed  using  these  MRI  data.  In  this 
study,  we  assess  the  ability  of  MREIT  to  monitor  changes  in  the  conductivity 
distribution  of  an  agarose  gel  phantom,  using  injected  current  pulses  of  900 
/xA.  The  phantom  initially  contained  a  distinct  region  of  high  sodium  chloride 
concentration  which  diffused  into  the  background  over  time.  MREIT  data 
were  collected  over  a  12  h  span,  and  conductivity  images  were  reconstructed 
using  the  iterative  sensitivity  matrix  method  with  Tikhonov  regularization. 
The  results  indicate  that  MREIT  was  able  to  monitor  the  changing  conductivity 
and  concentration  distributions  resulting  from  the  diffusion  of  ions  within  the 
agarose  gel  phantom. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

In  vitro  studies  have  shown  that  the  electrical  conductivity  of  malignancies  may  be  20  to 
40  times  higher  than  healthy  tissues  and  benign  formations  (Surowiec  et  al  1988,  Jossinet 
1998).  Therefore,  in  vivo  conductivity  imaging  could  aid  in  the  diagnosis  and  characterization 
of  suspicious  lesions.  Electrical  impedance  tomography  (EIT)  is  a  conductivity  imaging 
modality  that  consists  of  three  general  steps.  A  current  distribution  is  established  within  the 
object  of  interest,  the  effects  of  these  currents  are  measured,  and  then  these  measurement  data 
are  used  to  reconstruct  the  conductivity  distribution  within  the  object.  In  conventional  EIT, 
surface  electrodes  are  used  to  inject  current  into  the  object  and  to  acquire  voltage  measurements. 
This  technique  is  constrained  by  the  physically  limited  number  of  measurement  electrodes 
and  restriction  to  surface  measurements.  As  a  result,  it  suffers  from  poor,  non-uniform  spatial 
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Figure  1.  MREIT  pulse  sequence  diagram. 


resolution  and  low  sensitivity  of  the  voltages  to  conductivity  changes  in  the  interior  regions 
(Seagar  et  al  1987). 

In  magnetic  resonance  electrical  impedance  tomography  (MREIT),  MRI  is  used  to 
measure  the  effects  of  applied  currents  throughout  the  entire  object.  These  applied  currents 
generate  a  magnetic  flux  density,  of  which  the  component  parallel  to  the  main  static  field 
can  be  measured  using  an  appropriate  pulse  sequence  (Scott  et  al  1991,  Ider  and  Muftuler 
1997,  Mikac  et  al  2001).  Given  this  single  component  of  current-generated  magnetic  flux 
density,  several  algorithms  have  been  proposed  to  reconstruct  the  conductivity  distribution. 
These  include  the  sensitivity  matrix  method  (Ider  and  Birgul  1998,  Birgul  et  al  2003b),  V2Z?z 
algorithms  (Seo  et  al  2003,  Ider  and  Onart  2004)  and  WBZ  algorithms  (Park  et  al  2004a, 
2004b).  Several  MREIT  studies  using  biological  tissue  phantoms  have  been  reported  (Khang 
et  al  2002,  Lee  et  al  2003,  Oh  et  al  2005).  However,  these  studies  focused  on  reconstructing 
static  conductivity  distributions  using  injected  current  levels  of  at  least  12  mA.  The  ability  to 
monitor  changes  in  conductivity  over  time  could  provide  additional  diagnostic  information, 
such  as  in  monitoring  tumour  growth  and  treatment.  In  this  study,  we  assessed  the  ability  of 
MREIT  to  monitor  changes  in  conductivity  over  time  using  sub-milliamp  injected  currents 
more  appropriate  for  human  use.  Validating  this  ability  is  a  necessary  step  towards  the  imaging 
and  monitoring  of  human  subjects  using  MREIT. 

2.  Methods 

2.1.  Theory 

In  MREIT,  applied  currents  establish  a  conductivity-dependent  current  density  distribution 
within  an  object,  which  in  turn  generates  a  magnetic  flux  density  according  to  Ampere’s 
Law.  Thus,  magnetic  flux  density  measurements  contain  information  about  the  conductivity 
distribution.  In  this  study,  current  pulses  of  alternating  polarity  were  applied  to  a  test  phantom 
using  surface  electrodes.  For  alternating  currents,  the  component  of  magnetic  flux  density 
parallel  to  the  main  static  field  (z-component)  can  be  measured  by  MRI  using  a  modified 
spin-echo  sequence  (Mikac  et  al  2001).  The  pulse  sequence  used  in  this  study  is  shown  in 
figure  1 .  The  current-generated  magnetic  flux  density  introduces  a  phase  shift  (p  in  the  MR 
image,  given  as: 

<P(  r)  =  2]/Vrpuise£(r) 


(1) 
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where  y  is  the  gyromagnetic  ratio,  N  is  the  number  of  cycles  of  injected  current  such  that 
each  cycle  consists  of  one  positive  and  one  negative  pulse,  rpuise  is  the  duration  of  each  pulse, 
and  B(r)  is  the  amplitude  of  z-component  current-generated  magnetic  flux  density  at  point  r. 
Measurement  of  this  phase  shift  allows  for  calculation  of  B( r). 

For  a  given  set  of  electrodes  used  for  current  injection,  data  were  collected  twice,  each  with 
opposite  polarities  in  the  applied  current  pulses.  The  resulting  phase  maps  were  subtracted 
then  divided  by  two,  so  as  to  cancel  out  any  additional  phase  contributions,  such  as  those 
arising  from  small  imperfections  in  the  hardware  timing.  A  given  set  of  injecting  electrodes 
forms  an  ‘injection  profile’.  Data  from  at  least  two  different  injection  profiles  are  required  to 
reconstruct  a  unique  solution  (Kwon  et  al  2002).  Previous  MREIT  studies  have  shown  that 
using  just  two  profiles  is  sufficient  for  reconstructing  accurate  images.  This  is  a  significant 
advantage  over  conventional  EIT,  which  typically  requires  the  use  of  several  (16  or  more) 
electrodes. 

For  reconstructing  the  conductivity  distribution  from  the  z -component  magnetic  flux 
density  data,  this  study  utilized  the  sensitivity  matrix  method  (SMM).  The  object  domain  was 
discretized  into  a  mesh  of  first-order  triangular  elements,  chosen  for  compatibility  with  the 
finite  element  method  (Reddy  1993).  A  linear  relationship  between  conductivity  perturbations 
and  z-component  magnetic  flux  density  perturbations  was  then  assumed,  such  that: 

A  Bz  =  S  Act  (2) 

where  S  is  known  as  the  ‘sensitivity’  matrix.  Given  m  measurement  points  and  n  conductivity 
elements,  ABZ  is  an  m  x  1  vector,  Act  is  an  n  x  1  vector,  S  is  an  m  x  n  matrix,  and  the  matrix 
element  Sfj  is  the  change  in  the  zth  magnetic  flux  density  measurement  due  to  a  small  variation 
in  the  conductivity  of  the  7th  element.  An  initial  conductivity  distribution  crinitiai  was  assumed 
(e.g.  uniform  conductivity),  and  the  problem  linearized  around  this  initial  condition: 

final  ^initial  =  S(<7fmai  ^initial)  (3) 


where  initial  is  the  z-component  magnetic  flux  density  distribution  given  ainitiai  and  calculated 
using  the  finite  element  method  and  Biot-Savart  law,  B& nai  is  the  MRI  measured  z-component 
magnetic  flux  density,  crfinal  is  the  actual  (unknown)  conductivity  distribution,  and  S  is  the 
sensitivity  matrix  such  that: 


Sij  = 


dBi 

d(Tj 


(4) 


In  general,  S  is  an  ill-conditioned  matrix,  so  that  a  simple  least-squares  fit  cannot  be  used 
to  solve  for  afinal.  In  this  study,  the  linear  equation  was  solved  using  the  conjugate  gradient 
method  with  Tikhonov  regularization,  where  afinai  was  found  by: 


min{||SA<j  —  ABZ\\2  +  A||A<r||2}  (5) 

°final 


where  A  is  a  regularization  parameter  (Golub  et  al  1999).  The  regularization  parameter 
A  was  chosen  such  that  the  calculated  z-component  magnetic  flux  density  generated  by  the 
reconstructed  conductivity  distribution  was  closest  to  the  MRI-measured  magnetic  flux  density. 


mm{||5finai  -  [5(or)]a=afilul(A.)||}.  (6) 

After  finding  afinai,  if  the  change  in  conductivity  A  o  was  greater  than  some  predefined 
threshold,  then  this  <Tfinal  was  assigned  as  the  new,  updated  o'  initial,  and  the  process  was  iterated. 

The  SMM  reconstructs  relative  conductivities.  Thus,  the  true  conductivity  distribution 
a  tme  is  related  to  the  reconstructed  conductivity  distribution  crfinai  by  a  constant  scaling 
factor  K : 


^true  —  A  CTfina]  • 


(7) 
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Figure  2.  (a)  Side  view  and  (b)  top  view  of  the  phantom  used  to  monitor  the  diffusion  of  NaCl. 


Determining  the  absolute  conductivity  distribution  requires  an  additional  constraint,  such  as  a 
surface  voltage  measurement  (Birgul  et  al  2003a). 

2.2.  Phantom 

The  previously  outlined  methods  were  used  to  monitor  the  conductivity  distribution  within  an 
agarose  gel  phantom.  The  conductivities  of  several  4  mM  Q1SO4  and  2%  agarose  mixtures 
with  a  range  of  NaCl  concentrations  from  0%  to  1%  were  measured  using  the  4-electrode 
method  (Baumann  et  al  1997).  Within  this  range,  a  linear  relationship  was  found: 

a  =  AC +  E  (8) 

where  the  conductivity  a  is  in  units  of  S  m-1,  the  NaCl  concentration  C  is  in  units  of  g/ 100  ml 
water,  and  the  constants  A  and  E  are  1 .7  and  0.074  respectively.  For  larger  NaCl  concentrations, 
the  relationship  becomes  nonlinear  (Fuoss  and  Onsanger  1957). 

For  the  phantom,  a  hollow  acrylic  disc  with  an  inner  diameter  of  7  cm  was  filled  to  a 
thickness  of  1  cm  with  4  mM  Q1SO4  and  2%  (g/100  ml  water)  agarose.  Within  this  disc,  a 
smaller  cylindrical  region  of  12  mm  diameter  and  1  cm  thickness  was  filled  with  4  mM  Q1SO4, 
2%  agarose  and  0.83%  NaCl.  This  NaCl  concentration  was  selected  using  (8)  to  provide  a 
conductivity  contrast  of  20  relative  to  the  background.  An  actual  conductivity  contrast  of 
19.7  was  found  by  measuring  the  prepared  gels  using  the  4-electrode  method.  The  initial 
mixtures  within  each  of  these  two  regions  were  assumed  to  be  homogeneous.  The  phantom 
was  constructed  such  that  within  the  phantom,  the  overall  concentration  distribution  did  not 
vary  along  the  z-direction.  A  schematic  of  the  phantom  is  shown  in  figure  2. 

To  measure  the  conductivity  distribution  within  the  phantom  using  MREIT,  four  copper 
electrodes  each  6.35  mm  wide  were  placed  equidistant  along  the  inner  wall  of  the  acrylic  disc 
and  used  to  inject  currents  into  the  interior  region  (figure  2).  Using  (7)  and  (8),  measurement 
of  the  conductivity  distribution  was  used  to  map  the  NaCl  concentration  distribution.  The 
scaling  factor  K  was  determined  by  applying  the  constraint  that  the  total  amount  of  NaCl  in 
the  phantom  remained  constant: 

/  Cdv  =  constant.  (9) 

J  phantom 

Measurements  of  the  conductivity  distribution  over  time  were  used  to  monitor  the  diffusion 
process  occurring  within  the  phantom. 
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In  the  phantom’s  initial  state,  there  existed  a  distinct  region  containing  NaCl  and  a 
background  containing  no  NaCl.  Over  time,  the  NaCl  diffused  from  the  region  of  higher 
concentration  to  the  region  of  lower  concentration  as  governed  by  the  diffusion  equation: 


dc 

—  =  V  •  DVC 

dt 


(10) 


where  C  is  the  concentration  of  the  diffusing  substance  and  D  is  the  diffusion  constant.  Due 
to  the  geometry  of  the  phantom,  there  existed  no  concentration  gradient  along  the  z -direction 
within  the  phantom.  Thus  this  diffusion  process  was  independent  of  the  z -direction  and  can 
be  analysed  in  2D.  For  a  disc  of  radius  a  on  an  infinite  place  surface,  the  concentration  C  at 
radius  r  and  time  t  is  given  as: 


r  (  \ 

C(r)  =  - exp 

V  2Dt  F 


(11) 


where  C7  is  the  initial  uniform  concentration  of  the  disc  and  Iq  is  the  modified  Bessel  function 
of  the  first  kind  of  order  zero  (Crank  1975).  For  a  given  time  t  >  0,  the  peak  concentration  C 
is  located  at  the  centre  of  the  disc,  and  is  given  as 


C  =  Ci  \  1  —  exp 


£))■ 


(12) 


23.  Experiment 

The  agarose  gel  phantom  was  placed  within  a  4T  MRI  scanner  with  the  plane  of  the  disc 
aligned  perpendicular  to  the  main  static  field.  Currents  within  the  phantom  were  constrained 
within  this  plane  and  the  current  density  distribution  contained  a  negligible  z-component. 
As  a  result,  the  current-generated  magnetic  flux  density  within  the  phantom  was  primarily 
along  the  main  static  field  axis.  This  was  ideal  since  the  pulse  sequence  only  detects  the 
z-component  of  any  magnetic  flux  density.  Currents  were  delivered  to  the  electrodes  using 
thin  copper  wires.  These  wires  were  aligned  parallel  to  the  main  static  field  for  up  to  10  cm 
away  from  the  phantom,  such  that  the  magnetic  flux  density  they  generated  contained  a 
negligible  z-component  in  the  phantom,  and  thus  did  not  affect  the  MRI  measurements.  A 
current- source  circuit  was  constructed  and  used  to  generate  the  injected  currents  (Baumann 
et  al  1997).  Synchronization  of  the  circuit  output  with  the  pulse  sequence  was  controlled  by 
the  scanner  computer  through  a  digital  interface. 

An  MREIT  measurement  was  acquired  as  soon  as  possible  after  preparation  of  the  agarose 
gel  phantom  (30  min).  Additional  measurements  were  taken  at  1  h  intervals  after  phantom 
preparation  for  up  to  8  h,  followed  by  scans  at  10  and  12  h.  For  each  time  point,  data 
were  collected  for  two  injection  profiles,  using  pairs  of  electrodes  directly  opposite  of  each 
other  (figure  4).  Both  data  sets  were  used  simultaneously  in  the  conductivity  reconstruction 
algorithms.  For  each  injection  profile,  data  were  collected  using  the  previously  outlined 
pulse  sequence,  injecting  5  cycles  of  current  with  an  amplitude  of  900  /zA,  a  pulse  length  of 
3567  /zs  and  a  cycle  period  of  10  ms  (figure  1).  Other  scan  parameters  were:  TR  =  500  ms, 
TE  =  60  ms,  slice  thickness  =  5  mm,  FOV  =100  mm,  matrix  =  64  x  64,  NEX  =  8.  Raw  MRI 
data  were  exported  to  a  separate  computer  for  processing.  Magnetic  flux  density  calculations 
and  reconstruction  algorithms  utilizing  the  iterative  SMM  with  Tikhonov  regularization  were 
developed  using  MATLAB.  For  reconstruction,  a  circular  finite  element  mesh  with  1089  nodes 
and  2048  triangular  elements  was  utilized,  as  shown  in  figure  3(a).  For  each  timepoint,  the 
reconstruction  was  iterated  until  successive  iterations  resulted  in  not  more  than  a  1  %  change 
in  the  peak  conductivity  value. 
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(a)  (b) 


Figure  3.  (a)  Finite  element  mesh  used  for  reconstruction,  (b)  Spin  echo  image  of  the  phantom. 
The  smaller  circular  markers  around  the  periphery  were  used  to  locate  the  electrodes  during 
reconstruction. 
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Figure  4.  Reconstructed  z- component  magnetic  flux  density  distributions  at  the  first  time  point  for 
the  two  injection  profiles.  The  white  arrows  indicate  the  location  of  the  current  injecting  electrodes. 


3.  Results 

Figure  3(b)  shows  the  MR  magnitude  image  of  the  phantom  taken  using  a  standard  spin-echo 
sequence  just  after  at  the  first  time  point.  The  scan  parameters  were:  TR  =  1  s,  TE  =  20  ms, 
slice  thickness  =  5  mm,  FOV  =  100  mm,  matrix  =  256  x  256,  NEX  =  2.  Despite  two  distinct 
regions  of  different  NaCl  concentrations,  there  is  little  contrast  in  this  image. 

Figure  4  shows  the  calculated  z -component  magnetic  flux  density  at  the  first  time  point  for 
the  two  injection  profiles.  Figure  5  shows  the  reconstructed  conductivity  images,  obtained  by 
processing  the  magnetic  flux  density  data  with  the  iterative  SMM  with  Tikhonov  regularization. 
Unlike  the  MR  magnitude  image,  these  conductivity  images  show  a  clear  contrast  between 
regions  of  different  NaCl  concentrations.  For  each  time  point,  the  spatial  profile  centred 
on  the  initial  high  NaCl  concentration  region  was  calculated  by  averaging  the  radial  profiles 
taken  along  four  different  directions  as  shown  by  the  red  dotted  lines  in  figure  2(b).  The 
resulting  curves  are  plotted  in  figure  6.  Over  time,  the  higher  conductivity  region  broadened 
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Figure  5.  Reconstructed  conductivity  distribution  (a)  30  min,  (b)  1  h,  (c)  2  h,  (d)  3  h,  (e)  4  h, 
(f)  5  h,  (g)  6  h,  (h)  7  h,  (i)  8  h,  (j)  10  h  and  (k)  12  h  after  phantom  preparation. 
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Figure  6.  Comparison  of  the  spatial  profiles  of  the  inner  disc.  The  initial  profile  at  /  =  0  is  denoted 
by  the  dotted  line. 
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Figure  7.  Decay  of  the  peak  conductivity.  Reconstructed  values  denoted  by  the  blue  squares  were 
fitted  to  (12)  resulting  in  the  black  curve. 


and  decreased  in  magnitude,  and  the  border  between  initial  regions  became  less  distinct.  The 
change  in  the  conductivity  distribution  is  consistent  with  the  diffusion  of  NaCl  from  the  initial 
disc  region  into  the  background. 

To  quantify  this  diffusion,  the  peak  conductivity  value  for  each  time  point  was  scaled  using 
(9)  then  fitted  to  (12)  using  the  Curve  Fitting  Toolbox  of  MATLAB.  The  results  are  shown 
in  figure  7.  The  peak  values  appear  to  coincide  with  the  theoretical  curve,  indicating  that  the 
changing  reconstructed  conductivity  distributions  over  time  can  be  attributed  to  the  diffusion 
process.  An  experimental  diffusion  constant  of  D  =  7.68  x  10-6  cm2  s-1  was  obtained. 


4.  Discussion 

The  previous  results  indicate  that  MREIT  was  able  to  monitor  the  changing  conductivity 
distribution  resulting  from  the  diffusion  of  ions  within  the  agarose  gel  phantom.  The 
experimentally  determined  diffusion  constant  is  45%  smaller  than  the  previously  reported 
measurement  of  1.4  x  10-5  cm2  s-1  by  Schantz  and  Lauffer  (1962).  However,  their  study 
investigated  NaCl  in  an  estimated  1.5%  agarose  gel,  which  can  in  part  explain  their  faster 
diffusion.  The  smaller  diffusion  constant  for  the  MREIT  phantom  can  also  be  attributed  to 
its  finite  shape  and  confined  volume.  Equation  (11)  assumes  that  diffusion  occurs  on  an 
infinite  plane.  In  a  restricted  area,  NaCl  can  accumulate,  which  decreases  the  apparent  rate  of 
diffusion. 

In  using  the  SMM,  only  relative  changes  in  conductivity  can  be  found.  However,  in 
monitoring  the  phantom,  only  contrasts  in  conductivity  were  required  to  observe  the  ion 
diffusion.  In  applications  such  as  tumour  characterization,  differences  in  conductivity  between 
suspicious  lesions  and  surrounding  background  tissue  may  be  more  significant  than  absolute 
conductivity  values. 

While  data  were  collected  every  hour,  shorter  time  resolutions  are  possible  and  limited 
by  the  rate  at  which  data  can  be  acquired.  For  the  scan  parameters  used  in  this  study,  data 
for  one  injection  profile  were  collected  in  under  5  min.  In  addition  to  shorter  time  periods, 
MREIT  should  also  be  able  to  monitor  changes  in  conductivity  over  periods  much  larger  than 
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that  used  in  this  study.  Such  extended  durations  would  be  useful  when  monitoring  tumour 
growth  and  assessing  response  to  treatment.  For  example,  MREIT  could  be  used  to  evaluate  a 
patient’s  response  to  chemotherapy  by  measuring  the  (changing)  conductivity  of  a  suspicious 
lesion  over  several  weeks. 

Application  of  MREIT  to  human  studies  has  seemed  unpractical  due  to  the  large  amounts 
of  injected  currents  utilized  in  previous  studies.  Lower  current  levels  result  in  a  decrease 
in  the  signal-to-noise  ratio  (SNR)  of  the  acquired  magnetic  flux  density  maps  used  in  the 
reconstruction  process.  Oh  et  al  (2005)  recently  acquired  conductivity  images  of  a  biological 
tissue  phantom  in  a  3T  MRI  system  using  the  harmonic  Bz  algorithm.  Their  results  show 
significant  degradation  in  the  reconstruction  at  as  high  as  12  mA  of  injected  current.  This  is 
largely  due  to  the  use  of  the  derivative  operator  on  the  magnetic  flux  density  measurements, 
which  amplifies  errors  from  noise.  Other  reported  biological  tissue  phantom  studies  utilizing 
gradient-based  reconstruction  methods  used  currents  in  excess  of  24  mA.  The  results  of  this 
study  demonstrate  that  MREIT  is  possible  for  injected  currents  as  low  as  900  /jlA.  The  SMM 
appears  significantly  more  robust  when  reconstructing  lower  SNR  data.  Future  work  will 
search  for  the  minimum  levels  of  current  needed  for  reconstructing  accurate  conductivity 
images. 

Due  to  the  geometry  of  the  diffusion  phantom,  currents  were  restricted  to  flow  within 
a  slab  of  1  cm  thickness.  In  more  intricate  objects  such  as  the  human  body,  currents  may 
flow  throughout  the  entire  object.  Such  3D  current  flow  will  result  in  a  decrease  in  the 
measurable  magnetic  flux  density.  The  corresponding  reduction  in  the  SNR  presents  an 
additional  challenge.  To  improve  the  SNR  without  increasing  the  injected  current  level, 
future  work  will  optimize  hardware  and  pulse  sequence  parameters.  Use  of  a  reconstruction 
algorithm  capable  of  processing  lower  SNR  data,  such  as  the  SMM,  is  also  critical.  The 
SMM  has  already  been  applied  by  Muftuler  et  al  (2004)  in  reconstructing  the  conductivity 
distribution  of  an  in  vivo  rat  using  1  mA  of  injected  current. 

The  maximum  level  of  current  that  can  safely  be  applied  to  humans  depends  on  the 
frequency.  In  general,  the  safe  current  limit  increases  as  the  frequency  increases.  A  majority  of 
previously  reported  MREIT  studies  have  utilized  single  cycle,  long  current  pulses,  essentially 
classifying  them  as  dc  with  regards  to  safety  standards.  This  study  utilized  alternating  current 
in  the  form  of  multiple  cycles  of  shorter  alternating  polarity  pulses,  with  a  cycle  period  of 
10  ms  (primary  frequency  component  of  100  Hz).  While  safety  levels  at  100  Hz  are  only 
slightly  higher  than  dc,  the  concept  of  using  ac  over  dc  has  the  advantage  that  future  studies 
using  higher  frequencies  will  allow  for  safe  application  of  larger  current  levels.  Use  of 
higher  frequencies  will  also  enhance  the  capacitive  contributions  to  overall  impedance,  thus 
providing  an  additional  source  of  information  that  may  be  of  diagnostic  relevance.  Future 
work  will  test  this  pulse  sequence  and  develop  new  pulse  sequences  for  use  in  higher  frequency 
MREIT. 


5.  Conclusions 

By  utilizing  NaCl  diffusion  in  an  agarose  gel  phantom,  we  were  able  to  demonstrate  the 
ability  of  MREIT  to  monitor  changes  in  an  object’s  conductivity  distribution  over  time.  An 
injected  current  level  of  900  /xA  was  found  acceptable  for  reconstructing  accurate  images 
when  utilizing  the  iterative  sensitivity  matrix  method  with  Tikhonov  regularization.  This  is 
approaching  the  safe  current  levels  for  application  in  humans.  The  results  of  this  study  indicate 
that  human  MREIT  should  be  feasible  in  the  near  future,  and  that  this  modality  can  be  used 
not  only  for  acquiring  static  images,  but  also  used  for  performing  dynamic  studies. 
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Resolution  And  Contrast  In  Magnetic  Resonance 
Electrical  Impedance  Tomography  (MREIT) 
and  its  Application  To  Cancer  Imaging 

www.tcrt.org 

It  has  been  reported  that  the  electrical  impedance  of  malignancies  could  be  20-40  times 
lower  than  healthy  tissues  and  benign  formations.  Therefore,  in  vivo  impedance  imaging 
of  suspicious  lesions  may  prove  to  be  helpful  in  improving  the  sensitivity  and  specificity 
of  detecting  malignant  tumors.  Several  systems  have  been  developed  to  map  the  con¬ 
ductivity  distribution  inside  a  volume  of  tissue,  however  they  suffer  from  poor  spatial  res¬ 
olution  because  the  measurements  are  taken  only  from  surface  electrodes.  MRI  based 
impedance  imaging  (MREIT)  is  a  novel  method,  in  which  weak  electrical  currents  are 
injected  into  the  tissue  and  the  resulting  perturbations  in  the  magnetic  field  are  measured 
using  MRI.  This  method  has  been  shown  to  provide  better  resolution  compared  to  pre¬ 
vious  techniques  of  impedance  imaging  because  the  measurements  are  taken  from 
inside  the  object  on  a  uniform  grid.  Thus,  it  has  the  potential  to  be  a  useful  modality  that 
may  detect  malignancies  earlier.  Several  phantom  imaging  experiments  were  performed 
to  investigate  the  spatial  resolution  and  dynamic  range  of  contrast  of  this  technique.  The 
method  was  also  applied  to  a  live  rat  bearing  a  R3230  AC  tumor.  Tumor  location  was 
identified  by  contrast  enhanced  imaging. 

Key  words:  Electrical  Impedance  Tomography,  MRI,  Cancer  imaging. 

Introduction 

Breast  cancer  is  the  leading  cause  of  cancer  deaths  among  women  in  many 
parts  of  the  world  and  currently  no  single  imaging  modality  has  both  high  sen¬ 
sitivity  and  specificity  for  the  diagnosis  of  breast  cancer.  At  present,  well- 
established  breast  screening  methods  have  high  sensitivity  but  suffer  from 
poor  specificity.  X-ray  mammography,  which  is  accepted  as  the  gold  standard 
for  breast  cancer  screening,  provides  high  sensitivity  but  has  a  high  rate  of 
false-positives  (2-4).  Similarly,  the  sensitivity  of  breast  MRI  can  be  very  high, 
however  the  specificity  for  the  detection  of  abnormalities  is  variable  (5,6).  In 
addition  to  the  high  rate  of  false  positives,  some  of  these  techniques  also  fail 
to  detect  breast  cancer  in  some  cases,  such  as  patients  with  dense  breast  tis¬ 
sue,  which  is  common  among  younger  patients,  and  patients  undergoing  hor¬ 
mone  replacement  therapy  during  post-menopausal  period  (7).  New  imaging 
modalities  with  better  specificity  may  help  reduce  the  rate  of  false  positives, 
hence  eliminate  unnecessary  biopsies. 


Abbreviations:  MRI,  Magnetic  Resonance  Imaging;  EIT,  Electrical  Impedance  Tomography; 
MREIT,  Magnetic  Resonance  Electrical  Impedance  Tomography;  FEM,  Finite  Element  Method;  EIS, 
Electrical  Impedance  Scanning;  OPAMP,  Operational  Amplifier;  SVD,  Singular  Value 
Decomposition;  NEX,  Number  of  Excitations;  CE-MRI,  Contrast  Enhanced  MRI. 
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It  has  been  shown  that  electrical  properties  of  malignant  tis¬ 
sues  are  significantly  different  from  those  of  normal  and 
benign  tissues.  Surowiec  et  al.  (1)  have  reported  that  the 
electrical  impedance  of  malignant  tumors  decreases  by  a  fac¬ 
tor  of  20  to  40  with  respect  to  normal  or  benign  tissues.  This 
information  may  be  used  in  tumor  detection  and  characteri¬ 
zation.  Therefore,  the  spatial  distribution  of  conductivity 
measurements  (which  is  inversely  related  to  impedance)  can 
be  used  in  conjunction  with  other  imaging  modalities  to 
identify  tumors  and  achieve  higher  specificity  rates  com¬ 
pared  to  the  currently  used  techniques. 

Electrical  Impedance  Scanning  (EIS)  is  a  new  technique 
recently  introduced  to  aid  in  the  diagnosis  of  malignant  breast 
tumors  (12-15,  26).  The  device  that  uses  this  technique  gen¬ 
erates  a  map  of  the  conductivity  distribution  inside  the  breast 
tissue.  This  device  received  FDA  approval  recently  and  has 
been  introduced  to  the  market  (15).  A  probe  that  comprises 
of  an  array  of  electrodes  is  placed  over  the  breast  and  small 
amounts  of  current  are  injected  into  the  tissue  through  anoth¬ 
er  electrode  placed  in  the  palm  of  the  patient.  The  magnitude 
and  phase  of  the  current  is  measured  at  one  electrode  at  a  time 
and  all  electrodes  are  scanned  sequentially.  Using  the  distri¬ 
bution  of  electric  currents  in  this  electrode  array,  a  conduc¬ 
tivity  map  is  reconstructed;  but  it  does  not  provide  tomo¬ 
graphic  images.  Therefore,  the  conductivity  maps  recon¬ 
structed  using  this  technique  have  poor  spatial  resolution  and 
the  spatial  resolution  is  not  uniform  throughout  the  imaging 
region.  The  generated  maps  of  conductivity  are  insensitive  to 
tumors  deep  inside  the  breast,  which  are  located  3-3 .5cm 
away  from  the  surface.  Moreover,  tumors  just  under  the  nip¬ 
ple  are  confounded  by  the  high  conductivity  of  the  nipple. 

Electrical  Impedance  Tomography  (EIT)  was  developed  in 
the  1980s  to  reconstruct  the  conductivity  distribution  inside 
a  conducting  volume.  In  this  method,  a  current  distribution 
is  generated  inside  the  object  by  injection  or  induction  and 
peripheral  voltage  measurements  are  acquired  using  limited 
number  of  surface  electrodes  to  find  the  internal  conductivi¬ 
ty  distribution  (8).  This  method  also  suffers  from  poor  spa¬ 
tial  resolution  due  to  the  limited  number  of  voltage  measure¬ 
ments  taken  only  from  the  surface.  This  imaging  technique 
is  also  less  sensitive  to  conductivity  changes  deep  inside  the 
object.  Cherepenin  et  al.  (9)  have  recently  investigated  the 
potential  of  this  technique  in  breast  cancer  detection.  They 
have  developed  an  improved  EIT  system  with  256  electrodes 
placed  on  a  planar  surface  for  improved  spatial  resolution. 
They  collected  in  vivo ,  multi-slice  tomographic  impedance 
images  of  healthy  and  cancerous  breasts.  Two  remote  elec¬ 
trodes  were  placed  on  patient’s  extremities  for  current  injec¬ 
tion  and  voltage  reference.  This  allowed  using  the  assump¬ 
tion  of  spherical  equipotential  surfaces  near  the  electrode 
array.  Filtered  backprojection  along  these  equipotential  sur¬ 
faces  were  used  for  fast  3D  image  reconstruction.  This 


device  measured  only  the  magnitude  of  voltage  measure¬ 
ments,  so  the  phase  information  was  ignored.  This  setup 
helped  improve  image  quality  and  they  were  able  to  detect 
high  conductivity  tumor  formations  of  3-5  cm  in  size. 
However,  the  system  still  lacks  high  resolution  to  detect 
tumors  less  than  a  few  centimeters.  Moreover,  the  electrodes 
with  bad  contacts  cause  significant  image  artifacts.  They 
developed  a  thresholding  technique  to  detect  such  electrodes 
and  discard  information  from  them. 

Estrela  et  al.  (10),  and  Kerner  et  al.  (11)  have  studied  elec¬ 
trical  properties  of  breast  tumors  using  an  EIT  system  in 
which  both  the  magnitude  and  phase  of  input  current  and 
surface  voltage  measurements  were  recorded.  Sinusoidal 
input  currents  were  applied  where  the  frequency  was 
changed  between  lOKHz  up  to  about  1MHz  in  discrete  steps. 
This  allowed  imaging  both  conductivity  and  permittivity  of 
breast  tissues.  Typically,  16  surface  electrodes  were  used 
that  were  placed  on  a  circular  ring,  which  resulted  in  lower 
spatial  resolution  compared  to  (9).  However,  it  was  noted 
that  permittivity  information  was  more  useful  compared  to 
their  conductivity  counterparts  in  detecting  malignancies. 

MREIT  is  a  recently  developed  method  for  conductivity 
imaging  (16-21).  It  uses  the  magnetic  flux  density  measure¬ 
ments  acquired  from  MR  phase  images  to  reconstruct  con¬ 
ductivity.  Magnetic  flux  density  generated  by  applied  cur¬ 
rents  can  be  measured  with  high  spatial  resolution  using 
MRI.  Various  techniques  have  been  proposed  for  DC  (22), 
AC  (23,  24),  and  RF  (25)  currents.  Unlike  EIT,  the  spatial 
resolution  in  the  MREIT  is  position  independent.  However, 
it  should  be  noted  that  only  the  component  of  the  magnetic 
flux  density  in  the  direction  of  the  main  field  of  the  MRI  sys¬ 
tem  can  be  measured.  Therefore,  one  must  develop  a  tech¬ 
nique  to  solve  the  inverse  problem  of  finding  the  conductivi¬ 
ty  or  current  density  from  only  one  component  of  magnetic 
flux  density.  With  this  technique,  only  the  relative  conduc¬ 
tivity  values  can  be  reconstructed  from  using  only  the  mag¬ 
netic  flux  density  measurements.  In  order  to  find  the  absolute 
conductivity  values,  at  least  one  voltage  measurement  from 
the  boundary  is  required.  The  reconstruction  algorithms  can 
also  be  divided  into  two  groups  depending  on  the  data  type 
required.  The  first  group  uses  magnetic  flux  density  directly, 
whereas  in  the  second  one,  the  current  density  distribution  is 
required  in  image  reconstruction.  Algorithms  that  use  current 
density  distribution  requires  the  measurement  of  all  three 
components  of  the  magnetic  flux  density,  thus  rotation  of  the 
object  inside  the  magnet  is  required.  All  three  orthogonal 
components  of  the  magnetic  field  have  to  be  measured  in 
order  to  calculate  the  current  density  or  precise  conductivity 
distribution.  However,  this  is  not  a  practical  situation,  espe¬ 
cially  for  human  subjects.  Rotation  would  also  introduce 
problems  with  registration  of  measurements  taken  at  different 
orientations  of  any  object  under  investigation. 
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Thus,  a  method  that  uses  only  the  z-component  of  the  mag¬ 
netic  field  to  calculate  the  conductivity  distribution  (16,  20) 
was  adopted  in  this  study.  This  method  assumes  a  linear 
relationship  between  the  measured  magnetic  fields  and  rela¬ 
tive  conductivity  perturbations. 

In  the  presented  study,  the  MREIT  technique  was  tested  on 
various  conductivity  phantoms  to  investigate  its  spatial  reso¬ 
lution,  signal-to-noise  ratio  (SNR)  and  contrast- to-noise 
ratio  (CNR).  Various  phantoms  were  built  with  agarose  gels 
mixed  with  NaCl  to  create  compartments  with  known  con¬ 
ductivity  distributions.  Preliminary  data  was  also  collected 
in  vivo  from  Sprague-Dawley  rats  innoculated  with  the 
R3230  AC  tumor.  Contrast  Enhanced  MRI  images  were  also 
collected  to  localize  the  tumor  and  investigate  the  spatial 
correlation  between  the  outcomes  of  the  two  techniques. 

Materials  and  Methods 

Data  were  collected  by  a  4T  whole  body  MRI  system  that  has 
a  Magnex  magnet  (Magnex  Scientific  Inc.,  UK)  equipped 
with  a  whole  body  gradient  coil  set  (Tesla  Eng.  UK),  which 
provides  up  to  3G/cm  gradient  fields.  The  clear  bore  size  of 
the  magnet  with  gradient  assembly  is  650mm.  A  13  channel 
room-temperature  high-order  shim  system  with  MXA-13-4 
shim  power  supply  (Resonance  Research,  Billerica,  MA)  is 
also  available  to  minimize  field  inhomogeneities.  This  sys¬ 
tem  is  interfaced  with  a  MRRS  console  (Magnetic  Resonance 
Research  Systems,  Guildford,  UK)  that  has  broadband  RF 
transmit  and  receive  channels.  A  16  leg,  quadrature,  high- 
pass  birdcage  coil  with  10  cm  diameter  and  18  cm  length  was 
designed  and  built  in-house  for  the  MREIT  experiments. 

Pulse  Sequence 

The  pulse  sequence  used  for  the  MREIT  experiments  is  sim¬ 
ilar  to  a  fast  spin  echo  sequence,  where  a  train  of  180°  RF 
pulses  is  applied  following  a  90°  RF  pulse  (Fig.  1).  However, 
no  phase  encode  or  read-out  gradients  were  applied  between 
the  180°  RF  pulses  and  the  data  was  collected  with  a  single 
read-out  gradient  only  after  the  last  180°  RF  pulse.  Burst  sine 
wave  currents  were  injected  into  the  object  during  the  RF 
pulse  train  where  each  180°  RF  pulse  was  applied  at  the  zero 
crossing  of  the  sine  waves.  The  duration  of  each  180°  RF 
pulse  was  1 .3ms.  A  similar  sequence  was  proposed  by  Mikac 
et  al.  to  obtain  image  of  electric  currents  inside  an  object  (24). 

The  z-component  of  magnetic  field  (parallel  to  the  main  stat¬ 
ic  field)  generated  by  the  injected  currents  introduces  a  phase 
shift  that  accumulates  over  each  half-cycle  of  the  sine  wave. 
The  simplified  MRI  signal  equation  is  given  as: 

s(u,v)  =  Jj  M(x,v)  eiQ(x^  e>^xu  + y  v)  eiV$o(b(x,y)  cos  (®0)  dt  fajy 

[1] 


Figure  1:  Timing  diagram  of  the  pulse  sequence  used  in  the  MREIT 
experiments . 

M(x,y)  is  proportional  to  the  density  of  protons  within  a 
voxel  at  coordinates  (x,y).  The  phase  angle  0(xy)  accounts 
for  all  the  constant  phase  terms  including  static  field  inho¬ 
mogeneity  and  other  hardware  related  phase  delays.  The 
gyromagnetic  ratio  is  shown  by  y,  and  u  and  v  in  this  equa¬ 
tion  are  called  the  “spatial  frequencies”.  They  are  defined  as 
u=y*Gx*TGx  and  v=y*Gy*TGy,  where  TGx  and  TGx  repre¬ 
sent  the  duration  of  the  gradient  pulses  Gx  and  Gy,  respec¬ 
tively.  Therefore,  in  the  final  image  the  accumulated  total 
phase  at  a  pixel  at  location  (x,y)  due  to  the  magnetic  field 
generated  by  injected  currents  is: 

cp(x,y)  =  4  •  y  •  N  •  b(x,y)  /  oo  [2] 

where  N  is  the  number  of  cycles  of  injected  current,  b(x,y)  the 
amplitude  of  the  z-component  current-generated  magnetic 
field  at  point  (x,y),  and  oo  the  angular  frequency  of  the  inject¬ 
ed  current.  Hence,  measurement  of  this  phase  shift  allows  for 
calculation  of  the  (z-component)  magnetic  field  distribution. 

During  the  experiments,  the  data  was  collected  twice  with 
opposite  polarities  of  electric  current,  hence  changing  the 
polarity  of  current-generated  magnetic  field  distribution. 
When  resulting  phase  terms  were  subtracted,  those  coming 
from  the  static  field  inhomogeneities  and  delays  in  MRI 
hardware  were  eliminated,  leaving  only  the  terms  caused  by 
b{x,y).  To  obtain  the  magnetic  field  distribution  b{x,y)  from 
the  acquired  data,  the  MRI  images  were  first  reconstructed 
using  the  FFT.  The  resulting  phase  terms  from  the  two  sets 
of  images  were  subtracted.  Then,  the  magnetic  field  term 
was  calculated  using  the  expression  given  in  Equation  [2] . 

Hardware:  An  HP  ESG-4400B  signal  generator  produced 
the  sine  waves.  These  signals  were  synchronized  to  the 
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pulse  sequence  by  a  TTL  pulse  generated  by  the  scanner 
computer.  A  transconductance  amplifier  was  designed  and 
built  using  three  LM741  OPAMP  circuits  to  convert  the  volt¬ 
age  from  the  signal  generator  into  a  current  output.  This  cir¬ 
cuit  is  described  in  detail  in  (28),  which  was  replicated  from 
(29).  The  output  current  was  calibrated  to  lma/2V  and  an 
independent  measurement  of  current  was  not  used  in  the 
present  study.  The  whole  experimental  setup  is  illustrated  in 
Figure  2.  MRRS  console  controls  the  whole  experiment, 
which  generates  the  pulse  sequence,  acquires  the  incoming 
data,  and  synchronizes  the  external  units.  After  data  collec¬ 
tion,  raw  MRI  data  (k- space)  was  exported  to  another 
Pentium  PC  for  off-line  processing. 


Figure  2:  Schematic  of  the  experimental  setup.  AH  is  the  animal  holder,  which 
is  placed  inside  the  MRI  RF  coil. 


Conductivity  Image  Reconstruction 

Sensitivity  Matrix  Method:  Since  the  relation  between  the 
conductivity  distribution  and  the  magnetic  fields  generated  by 
the  internal  current  distribution  is  nonlinear,  one  has  to  either 
employ  iterative  techniques  or  try  to  linearize  around  an  initial 
value.  In  the  presented  study,  the  second  approach  was  adopt¬ 
ed,  in  which  a  linear  relationship  between  the  conductivity 
perturbations,  A  a,  and  magnetic  field  perturbations,  AB,  was 
assumed  (16).  This  relationship  is  given  in  Equation  [3]  as: 


In  the  first  step,  the  electric  potential  distribution  cj)(v,y,z)  in 
the  imaging  slice  is  determined  by  solving  the  Poisson’s 
equation  with  Neumann  boundary  value  conditions.  This 
problem  is  nonlinear  in  terms  of  conductivity  since  the  elec¬ 
tric  potential  itself  is  also  conductivity  dependent.  The 
Finite  Element  Method  (FEM)  was  used  to  calculate  the  dis¬ 
tribution  of  the  potential.  Once  is  calculated,  the  electric 
field  and  current  density  distribution  inside  the  imaging 
region  may  be  found  using  the  equations: 


E  =  -V(|) 
J=  a  E 


[4] 


The  magnetic  field  generated  by  this  ohmic  current  can  be 
calculated  by  using  Equation  [4]  and  the  Biot-Savart  law. 
The  differential  magnetic  field  caji  be  written  in  terms  of 
the  differential  current  element  Idl  as, 


cffxR  \ 


[5] 


where  \i0  is  the  permeability  constant  and  R  is  the  vector 
from  the  source  point  at  (x’ ,  y’ ,  z’)  to  the  field  point  (x,  y, 
z).  The  relationship  between  the  current  and  magnetic 
field  is  linear  and  the  total  magnetic  field  can  be  obtained 
by  integrating  Equation  [5]  over  all  sources.  Using 
Equations  [4-5],  the  magnetic  field  can  be  calculated  for 
an  assumed  conductivity  distribution. 


The  sensitivity  matrix  can  be  calculated  using  either  a 
numerical  or  a  semi-analytical  approach.  In  the  numerical 
approach,  each  column  of  the  sensitivity  matrix  is  calculat¬ 
ed  separately  by  solving  the  forward  problem  by  changing 
the  conductivity  of  a  single  element.  This  approach  requires 
the  repetitive  solution  of  the  forward  problem  and  increases 
the  execution  time.  For  faster  reconstruction,  the  semi-ana¬ 
lytical  method  was  used  for  the  calculation  of  the  sensitivi¬ 
ty  matrix  as  explained  in  (20). 


AB  =  S  •  Aa  [3] 

If  there  are  n  conductivity  elements  and  m  measurement 
points,  then  Aa  is  an  n  x  1  vector,  AB  is  an  m  x  1  vector  and 
the  “sensitivity  matrix”  S  is  an  m  x  n  matrix.  The  entry  of 
the  sensitivity  matrix  at  ith  row  and  jth  column  denotes  the 
change  in  the  ith  field  measurement  due  to  a  small  change  in 
conductivity  of  the  jth  element. 

In  order  to  calculate  the  “sensitivity  matrix”  S  one  must  solve 
the  forward  problem  for  a  given  conductivity  distribution  and 
current  injection  scheme.  First,  the  initial  magnetic  field  Bj  is 
calculated  from  an  initial  estimate  of  conductivity  (Oj).  Then, 
S  is  calculated  by  perturbing  each  element  of  Qj  one  by  one 
and  finding  the  corresponding  magnetic  field  distribution. 


Once  the  sensitivity  matrix  S  is  calculated  for  certain 
geometry,  the  conductivity  distribution  can  be  approxi¬ 
mated  by, 

Aa  =  S"1  •  AB  [6] 

Where  AB  =  Bmeas  -  B-  Bmeas  is  the  magnetic  field  meas¬ 
urements  obtained  from  MRI  and  Bj  is  the  initial  value  of 
the  magnetic  field  corresponding  to  the  initial  conductiv¬ 
ity  distribution  aj.  The  generalized  inverse,  S"1,  is  found 
by  using  singular  value  decomposition  (SVD).  This 
method  has  been  tested  both  by  simulation  and  experi¬ 
mental  phantom  studies  and  shown  that  a  good  approxi¬ 
mation  to  the  relative  conductivity  distribution  map  can 
be  obtained  (16,  20). 
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Phantom  Experiments 

A  series  of  phantom  studies  were  conducted  to  test  the  spa¬ 
tial  resolution  and  dynamic  range  of  contrast  for  MREIT. 
Several  agarose  gel  phantoms  with  different  conductivity 
distributions  were  prepared  for  the  tests.  The  gels  were 
placed  inside  an  acrylic  cylinder  with  an  inner  diameter  of 
7cm  and  height  of  1cm.  The  conductive  gels  consisted  of 
2%  (g/lOOmL)  agarose  and  varying  concentrations  of  NaCl. 
During  imaging,  the  axis  of  the  cylinder  was  placed  parallel 
to  the  z-axis  (direction  of  the  MRI  magnetic  field).  Four 
electrodes  made  of  copper  foil,  each  3mm  wide,  were  placed 
at  0°,  90°,  180°  and  270°  along  the  inner  wall  and  used  to 
inject  currents  into  the  interior  region.  For  each  phantom, 
data  were  collected  twice  for  two  current  injection  schemes 
and  used  simultaneously  in  conductivity  reconstruction.  The 
first  time,  current  was  injected  between  the  electrodes  at  0° 
and  180°,  and  the  second  time  between  the  ones  at  90°,  and 
270°.  In  the  preliminary  studies  conducted,  it  has  been 
empirically  verified  that  this  scheme  improves  the  sensitivi¬ 
ty  and  spatial  resolution  compared  to  using  a  single  pair  of 
electrodes.  In  all  phantom  studies,  current  carrying  wires 
were  placed  precisely  in  the  z-direction  to  eliminate  any  con¬ 
tribution  from  the  currents  flowing  in  those  wires.  Wires 
were  mounted  on  acrylic  support  beams  for  rigidity  and  sev¬ 
eral  RF  chokes  and  low  pass  filters  were  placed  along  these 
wires  to  suppress  RF  coupling. 

For  contrast  studies,  a  16mm  diameter  inner  disk  of  gel  sur¬ 
rounded  by  a  background  of  1%  NaCl  was  used.  Two  exper¬ 
iments  were  carried  out  with  two  similar  contrast  phantoms, 
where  the  inner  disk  contained  either  4%  or  10%  NaCl.  In 
the  preliminary  tests,  it  has  also  been  experimentally  verified 
that  conductance  scales  sufficiently  linearly  with  NaCl  con¬ 
centration.  Therefore,  the  conductivity  contrast  values  were 
approximately  1:4  and  1:10,  respectively. 

To  study  spatial  resolution,  two  different  phantoms  were 
imaged.  The  first  phantom  was  prepared  by  placing  two  hol¬ 
low  nylon  disks  (each  16mm  in  diameter)  inside  the  larger 
disk.  They  were  separated  by  17mm  center  to  center.  The 
nylon  shells  acted  as  insulators,  and  each  small  disk  was 
filled  with  the  same  gel  as  the  surrounding  background.  The 
second  phantom  was  prepared  with  the  same  geometry  as  the 
first  one  but  this  time  the  inner  disks  were  agar  gels  with  a 
conductivity  ten  times  higher  than  the  background. 

For  each  phantom  configuration,  two  separate  data  sets  were 
collected  by  injecting  6  cycles  of  2mA  or  4mA  (rms)  200Hz 
current  into  the  phantom.  All  data  were  collected  using  the 
previously  outlined  pulse  sequence  with  the  following 
parameters:  TR=500ms,  TE=40ms,  slice  thickness=lcm, 
FOV=10cm,  data  matrix  64  x  64  and  NEX=8.  A  single  slice 
was  collected  for  each  experiment.  The  z-component  of  cur¬ 


rent-generated  magnetic  field  distributions  were  calculated 
from  the  resulting  data  and  the  conductivity  distributions 
computed  using  the  method  outlined  above. 

Effect  of  Initial  Conductivity  Estimate 

As  described  in  the  methods  section,  conductivity  reconstruc¬ 
tion  requires  an  initial  estimate  of  conductivity  distribution. 
We  performed  simulation  studies  to  investigate  the  effect  of 
this  initial  estimate  on  reconstructed  conductivity  and  con¬ 
trast.  Four  cylindrical  phantoms  were  generated,  each  with  a 
small  disc  inside,  which  was  placed  off-center.  All  four  phan¬ 
toms  were  generated  with  a  background  conductivity  of 
0.002S/cm  and  the  conductivity  of  the  small  disk  was  assigned 
either  0.5,  2,  4  or  10  times  the  background  conductivity. 

Each  of  the  four  phantom  data  were  reconstructed  three 
times  with  uniformly  distributed  initial  conductivity  esti¬ 
mates  of  0.001  S/cm,  0.002S/cm  and  0.004S/cm,  respective¬ 
ly.  Singular  values  below  1%  were  truncated  and  1089 
nodes  and  2048  first  order  triangular  elements  were  used  in 
all  finite  element  models  used  in  reconstructions. 

Animal  Experiments 

An  MREIT  experiment  was  also  performed  for  in  vivo 
impedance  imaging  of  a  tumor-bearing  rat.  For  this  pur¬ 
pose,  a  special  animal  holder  was  prepared  from  acrylic 
sheets.  This  holder  helped  keep  the  animal  stationary  dur¬ 
ing  imaging.  It  also  allowed  for  the  placement  of  the  elec¬ 
trodes  in  consistent  positions  for  longitudinal  studies. 
Electrodes  were  placed  on  acrylic  hollow  tubes  filled  with 
CuS04  solution  to  detect  the  electrode  positions  precisely 
in  the  images.  Initial  experience  has  shown  that  the  precise 
localization  of  electrode  positions  is  essential  for  accurate 
reconstruction  of  impedance  images.  Current  carrying 
wires  ran  along  these  tubes,  which  were  in  z-direction. 
This  is  essential  to  minimize  interference  from  the  magnet¬ 
ic  fields  generated  by  current  in  the  wires.  Prior  to  imag¬ 
ing,  the  animal  was  anesthetized  by  IV  injection  of  keta¬ 
mine  and  xylazine  and  then  placed  inside  the  holder.  The 
electrodes  were  covered  with  a  thin  layer  of  conductive  gel 
to  provide  good  electrical  contact.  The  skin  areas  of  con¬ 
tact  were  shaved  for  better  conductance. 

An  anatomical  image  was  collected  using  FSE  sequence 
prior  to  the  MREIT  images.  A  single  slice  with  6mm  thick¬ 
ness  was  collected  from  the  same  anatomical  location  as  the 
MREIT  image.  The  data  matrix  was  256  x  256,  FOV=10cm, 
TR=4s,  TE=20ms/  100ms  NEX=4.  MR-EIT  images  were 
collected  using  the  previously  outlined  pulse  sequence  with 
TR=500ms,  TE=30ms,  NEX=8  (signal  averages),  64  x  64 
data  matrix,  FOV=10cm,  slice  thickness=6mm,  with  an  AC 
current  of  1mA  peak,  200Hz,  4  cycles. 
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Figure  3a:  Spin-Echo  image  of  1:4  contrast  phantom.  The  higher  con¬ 
ductivity  region  is  marked  with  red  dashed  lines  and  it  is  barely  distin¬ 
guishable  in  the  MR  image. 


Figure  3c:  Plots  of  reconstructed  conductivity  ratio  versus  actual  conduc¬ 
tivity  ratio  obtained  from  MREIT  images.  Plots  are  given  for  both  of  the 
injected  current  amplitudes.  Error  bars  are  also  included  in  these  plots. 


Figure  3b:  Conductivity  images  of  1:10  contrast  (top  row)  and  1:4  contrast 
(bottom  row)  with  2mA  and  4mA  current  injection  cases. 


Figure  3d:  Plots  of  reconstructed  conductivity  ratio  versus  actual  conduc¬ 
tivity  ratio  in  the  phantom  data  used  in  MREIT  simulations.  Three  plots  are 
given  for  three  different  initial  conductivity  estimates  of  0.001,  0.002  and 
0.004S/cm.  Actual  background  conductivity  was  given  as  0.002S/cm. 


Table  I 

Mean  and  standard  deviation  of  relative  conductivity  (oir,  o0r)  values  in  contrast  phantoms 
calculated  by  the  MREIT  method.  Results  are  reported  for  both  injected  current  amplitudes. 

a0:Gi  =  1:4  phantom  (m±std)  g0:o,  =  1:10  phantom  (m±std) 
I=2mA  I=4mA  I=2mA  I=4mA 

Inside  small  disk  (a,)  235.4±10.5  229.6±18.4  233.4±15.4  226.6±21.2 

Background  (Gor) 81.2±18.3  63.8+14.7  50.1±17.9  36.1±12.3 
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Figure  4a:  Conductivity  images  of  resolution  phantom  with  insulating  discs.  Both, 
2mA  and  4mA  current  injection  cases  were  shown.  Small,  dark  circular  areas  show 
high  resistance  compartments  (insulating  disks). 
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Figure  4b:  A  profile  taken  diagonally  from  the  conductivity  image 
of  the  first  resolution  phantom  going  through  the  centers  of  the 
small  disks.  This  profile  is  given  for  the  4mA  current  injection  case 
and  the  location  is  shown  with  blue  dashed  line  in  Figure  4a. 


Figure  4c:  Conductivity  images  of  resolution  phantom  with  high  conduc¬ 
tivity  discs.  Only  the  4mA  current  injection  case  is  shown.  Yellow-white 
circular  areas  show  low  resistance  compartments. 
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Figure  4d:  A  profile  taken  diagonally  from  the  conductivity  image  of  the 
second  resolution  phantom  shown  in  Figure  4c. 


Figure  5:  Axial  slices  from  an  SD  rat  inoculated  with  a  R3230  AC  tumor 
showing:  (a)  Conductivity;  (b)  CE-MRI;  and  (c)  spatial  correlation  of  CE- 
MRI  and  MREIT.  The  tumor  had  spread  to  multiple  foci  surrounding  the 
body.  In  (a),  red  shows  high  conductivity  and  blue  shows  low  conductivity 
regions;  whereas  red  regions  in  (b)  shows  enhancement  by  Gd-DTPA. 
Areas  that  have  both  high  conductivity  and  also  enhanced  by  contrast  agent 
are  shown  in  hot  colors  in  (c). 


Figure  5d:  FEM  mesh  constructed  from  the  anatomical  MR  Image. 
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Contrast  enhanced  MR  images  were  also  collected  (CE- 
MRI)  using  Gd-DTPA,  which  is  a  well-established  method 
to  detect  malignant  tissues.  Since  vascular  growth  is  great¬ 
ly  enhanced  in  tumor  structures,  contrast  agent  uptake  of 
tumor  sites  increases  with  respect  to  normal  tissues.  In 
addition  to  that,  interstitial  compartment  in  tumors  are  large 
compared  to  normal  structures,  so  the  wash-out  of  contrast 
agent  is  also  slower  in  tumors.  Therefore,  if  two  images  are 
collected,  one  pre-contrast  and  the  other  post-contrast,  the 
difference  image  yields  enhanced  pixel  intensity  in  the  areas 
of  tumor  growth.  In  this  study,  the  goal  was  to  verify  the 
potential  of  MR-EIT  to  detect  tumors.  Contrast  enhanced 
images  were  collected  using  a  GE  sequence  with  a  64  x  64 
data  matrix,  FOV=10cm,  6mm  slice  thickness,  TR=150ms, 
TE=5ms,  and  45°  flip  angle.  One  pre-contrast  image  was 
acquired  before  Gd-DTPA  injection  and  a  post-contrast 
image  was  collected  3  minutes  after  injection. 

Results 

Conductivity  images  were  reconstructed  using  the  method 
outlined  in  the  Methods  section.  In  the  calculation  of  the 
inverse  of  S  using  SVD,  the  highest  200  singular  values  were 
used  and  the  rest  were  truncated  to  minimize  noise  while 
preserving  essential  information.  This  value  was  also  deter¬ 
mined  experimentally  to  be  an  optimum  level  for  truncation. 
In  the  FEM  mesh,  512  triangular  elements  were  used. 

Impedance  Contrast  Phantoms 

Figure  3a  illustrates  the  standard  spin-echo  MRI  image 
acquired  from  the  contrast  phantom  (with  1 :4  conductivity 
contrast).  As  seen  in  the  Figure,  it  is  difficult  to  distinguish 
the  two  compartments  in  the  MRI  image.  The  top  row  in 
Figure  3b  shows  conductivity  images  obtained  from  the 
1:10  contrast  phantom  with  two  separate  current  ampli¬ 
tudes.  Similarly,  the  bottom  row  of  Figure  3b  shows  con¬ 
ductivity  image  of  the  1:4  contrast  phantom.  Higher  con¬ 
ductivity  regions  are  greatly  enhanced  in  the  resulting 
images.  The  means  of  the  reconstructed  relative  conductiv¬ 
ity  values  in  the  two  compartments  inside  the  contrast  phan¬ 
toms  were  calculated  and  given  in  Table  I. 

Figure  3c  shows  the  relationship  between  relative  object 
conductivity  values  and  relative  MREIT  image  conductivity 
values  as  measured  from  the  images  given  in  Figure  3b  (and 
Table  I).  Since  the  background  conductivity  is  kept  the  same 
for  1:4  and  1:10  contrast  phantoms,  all  conductivity  values 
reported  in  Table  I  were  normalized  by  their  corresponding 
background  values  (aor)  for  this  graph.  Thus,  the  conduc¬ 
tivity  image  contrast  is  obtained  by  dividing  the  mean  con¬ 
ductivity  in  small  disk  (air)  by  the  background  mean  value 
(aor).  Using  these  measurements,  the  relative  changes  in 
object  conductivity  and  corresponding  conductivity  contrast 


in  the  MREIT  images  were  plotted.  The  plots  are  given  for 
both  2mA  and  4mA  injected  current  cases.  A  non-linear  rela¬ 
tionship  is  observed  in  these  plots.  Moreover,  the  difference 
between  the  two  curves  was  statistically  significant. 

Effect  of  Initial  Conductivity  Estimate  on  Contrast 

The  reconstructed  images  in  simulation  studies  were  analyzed 
to  assess  the  effect  of  initial  estimate  in  conductivity  image 
contrast.  Mean  values  in  the  background  and  inner  disc  were 
calculated  for  each  of  the  generated  1:0.5,  1:2,  1:4  and  1:10 
contrast  phantoms.  Those  values  are  calculated  for  each  of 
the  three  different  initial  estimates  of  0.001,  0.002  and 
0.004S/cm.  Reconstructed  conductivity  ratio  is  defined  as  the 
mean  value  of  the  inner  disk  divided  by  the  mean  value  of  the 
background  in  the  reconstructed  images.  For  each  of  the  three 
initial  estimates,  a  graph  of  actual  conductivity  ratio  versus 
reconstructed  conductivity  ratio  was  plotted  in  Figure  3d. 

Spatial  Resolution  Phantoms 

Figure  4a  shows  the  conductivity  image  of  the  “spatial  reso¬ 
lution  phantom”  with  two  identical  nylon  disks.  The  results 
of  two  experiments  with  2mA  and  4mA  current  injections  are 
given  in  this  Figure.  The  gradually  decreasing  distance 
between  the  two  small  disks  was  used  to  assess  the  spatial 
resolution  capability  of  conductivity  images.  A  cross  sec¬ 
tional  profile  taken  diagonally  from  Figure  4a  is  shown  in 
Figure  4b.  This  profile  was  taken  from  a  cut  that  passed 
through  the  centers  of  the  small  disks.  Even  though  the  con¬ 
trast  was  somewhat  reduced  due  to  low-pass  filtering  of  the 
point- spread  function,  the  1mm  separation  could  still  be 
observed.  On  the  other  hand,  the  conductivity  phantom  with 
ten  times  higher  conductivity  does  not  provide  the  same  pre¬ 
cision  of  detection  capability  (Figure  4c).  The  1mm  separa¬ 
tion  is  hardly  distinguishable  in  this  figure. 

In  vivo  Imaging 

The  first  step  in  animal  imaging  is  the  construction  of  the  finite 
element  mesh  using  an  anatomical  high  resolution  MRI  image. 
The  FEM  mesh  shown  in  Figure  5d  consists  of  702  nodes  and 
1268  first  order  triangular  elements.  Exact  electrode  locations 
are  found  using  the  markers  in  the  animal  setup  and  used  for  the 
boundary  condition  determination.  The  relative  conductivity 
distribution  was  computed  as  described  above  and  overlaid  on 
the  anatomical  image  (Figure  5a).  The  resulting  images  show 
the  higher  conductivity  regions  with  high  contrast.  Similarly, 
the  contrast  enhancement  by  Gd-DTPA  is  illustrated  in  Figure 
5b.  To  investigate  how  the  two  images  spatially  correlate, 
another  image  was  generated  by  masking  the  CE-MRI  image 
with  the  MREIT  image  regions  that  have  conductivity  values  in 
the  upper  one-thirds  of  the  full  range.  Figure  5c  shows  regions 
that  have  high  conductivity  and  also  enhanced  by  Gd-DTPA. 
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Discussion 

In  this  preliminary  study,  the  potential  of  the  MREIT  method 
for  obtaining  in  vivo  conductivity  information  with  good 
spatial  resolution  and  sensitivity  was  demonstrated,  which 
may  be  helpful  in  identifying  malignant  tumors. 

In  the  results  presented  here,  a  fixed  frequency  of  200Hz 
was  used  for  the  applied  electrical  currents.  This  choice 
was  based  on  previous  studies  that  investigated  electrical 
properties  of  tumors  (7).  Those  studies  have  reported  that  a 
higher  conductivity  contrast  was  obtained  with  the  applica¬ 
tion  of  200Hz  currents. 

In  contrast  phantom  tests,  conductivity  ratios  of  1 :4  and  1 : 10 
were  used.  The  method  presented  here  clearly  detected  these 
conductivity  perturbations  in  the  objects.  Based  on  the  range 
of  conductivity  changes  of  malignancies  compared  to  those 
of  benign  or  healthy  tissues  that  were  reported  in  the  litera¬ 
ture  (1),  MREIT  should  offer  sufficient  sensitivity  for  the 
detection  of  malignancies.  Although  the  minimum 
detectable  conductivity  perturbation  by  MREIT  has  not  been 
tested,  from  the  results  presented,  it  may  be  anticipated  that 
much  smaller  changes  could  also  be  detected. 

Since  the  reconstruction  method  used  here  is  based  on  the 
linearity  assumption  between  conductivity  perturbations  and 
measured  magnetic  fields,  it  only  gives  precise  relative  con¬ 
ductivity  measurements  for  small  conductivity  perturbations 
within  the  object  (or  body).  In  the  results  presented  here,  the 
ratio  of  conductivity  values  measured  in  various  compart¬ 
ments  in  the  reconstructed  images  did  not  reflect  the  actual 
conductivity  contrast.  Since  relatively  large  conductivity 
changes  were  used  in  the  phantoms,  it  was  observed  that  the 
conductivity  contrast  in  the  images  was  lower  than  the  actu¬ 
al  ratios  due  to  the  nonlinear  nature  of  the  actual  relationship. 
The  non-linear  relationship  between  conductivity  perturba¬ 
tions  in  the  object  versus  perturbations  in  reconstructed  EIT 
images  was  studied  in  detail  by  Ider  et  al.  (27).  These  prob¬ 
lems,  arising  from  the  nonlinearity  of  reconstruction  method, 
can  be  minimized  by  using  iterative  reconstruction  tech¬ 
niques  with  regularization.  One  such  technique  called  the 
weighted  regularized  least  square  method  (WRLSM)  was 
used  by  Oh  et  al.  (30)  to  reconstruct  conductivity  images 
using  MRI.  Although  they  applied  this  method  to  harmonic 
Bz  algorithm ,  it  can  be  easily  adopted  to  our  reconstruction 
technique.  Harmonic  Bz  algorithm  requires  the  second  par¬ 
tial  derivatives  of  the  magnetic  field  Bz  (V2#z),  therefore,  it 
is  vulnerable  to  the  noise  in  the  Bz  measurements.  Spatial 
dependence  of  accuracy  was  observed  in  Oh’s  study,  where 
the  two  profiles  from  different  sections  of  conductivity 
images  revealed  different  accuracy  with  respect  to  the  true 
conductivity  profiles.  Moreover,  they  have  used  significant¬ 
ly  higher  current  amplitude  of  26mA,  compared  to  the  val¬ 


ues  of  2  and  4mA  used  in  the  presented  study.  One  advan¬ 
tage  of  our  reconstruction  technique  is  that  it  is  less  vulnera¬ 
ble  to  noise  compared  to  the  harmonic  Bz  algorithm,  because 
of  the  fact  that  it  does  not  require  derivatives  of  Bz. 
Therefore,  an  iterative  WRLSM  technique  is  expected  to 
improve  accuracy  of  relative  conductivity  images  obtained 
in  our  studies.  Our  group  is  currently  working  on  the  imple¬ 
mentation  of  such  a  technique.  In  addition  to  that,  it  can  also 
be  observed  that  the  contrast  in  the  conductivity  images 
decreased  when  the  amplitude  of  the  current  was  decreased 
from  4  to  2mA.  Since  the  non-linear  relationship  between 
magnetic  field  distribution  and  conductivity  values  is 
assumed  to  be  linear  in  the  sensitivity  matrix  method,  the 
degree  of  deviation  from  linearity  will  vary  with  changing 
current  amplitude.  As  seen  in  the  Figure,  the  curve  with 
2mA  current  amplitude  is  slightly  more  linear  compared  to 
the  4mA  current  injection  case  but  with  reduced  slope,  i.e., 
reduced  contrast.  This  aspect  could  be  an  important  consid¬ 
eration  in  longitudinal  studies  where  the  change  in  conduc¬ 
tivity  in  a  tumor  is  investigated  over  time. 

Simulation  studies  have  shown  that  the  relative  conductivity 
contrast  in  reconstructed  images  is  also  dependent  on  the  ini¬ 
tial  conductivity  estimate.  Slight  deviations  in  curves  shown 
in  Figure  3d  were  observed  within  the  range  of  initial  values 
tested.  Therefore,  the  initial  estimate  of  conductivity  in  sen¬ 
sitivity  matrix  reconstruction  will  have  an  impact  on  the 
resultant  image  contrast. 

Two  simple  phantoms  were  prepared,  where  the  separation 
between  two  closely  separated  non-conductive  or  highly 
conductive  disks  were  used  to  assess  resolution.  It  was 
found  that  the  closest  distance  of  1mm  was  still  detectable 
in  the  images  where  non-conducting  disks  were  used.  On 
the  other  hand,  the  1mm  separation  was  hardly  distinguish¬ 
able  in  the  phantom  with  highly  conductive  disks. 
Therefore,  the  relative  conductivity  of  inner  structures 
plays  an  important  role  in  the  resolution,  as  expected. 
Current  density  is  an  important  factor  that  determines  sen¬ 
sitivity  in  this  imaging  technique.  High  current  density  in  a 
particular  area  will  increase  sensitivity  in  that  area.  Thus, 
in  the  non-conducting  object  case,  there  will  be  high  current 
density  between  the  two  insulators,  whereas,  in  the  high- 
conductivity  object  case,  most  of  the  current  density  will  be 
concentrated  inside  the  relatively  large  objects  and  suppress 
the  perturbation  of  small  separation.  Moreover,  a  certain 
amount  of  diffusion  between  the  high-conductivity  objects 
and  the  low-conductivity  separation  can  be  expected,  which 
will  adversely  affect  the  resolution.  To  minimize  that 
effect,  images  are  collected  immediately  after  the  phantom 
preparation.  For  better  assessment  of  resolution,  diffusion 
has  to  be  taken  into  account.  Similarly,  it  may  be  expected 
that  the  detectability  of  the  object  separation  will  also  be 
affected  by  the  size  of  the  inner  disks.  Construction  of 
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sophisticated  high-resolution  phantoms  with  agarose  gel  is 
a  mechanically  challenging  task  due  to  the  lack  of  rigidity 
in  the  gel  structure.  Small  structures  can  easily  melt  and 
blend  into  each  other  making  boundaries  somewhat 
ambiguous.  Although  more  complicated  phantoms  could 
be  built  with  different  materials  to  investigate  the  spatial 
resolution  of  the  method,  these  two  simple  phantoms  illus¬ 
trate  how  well  the  separation  between  objects  of  different 
conductivity  could  be  resolved.  Generally,  a  high  contrast 
object  is  required  to  assess  the  spatial  resolution  of  an  imag¬ 
ing  technique,  since  it  is  affected  less  by  noise;  therefore, 
the  test  with  insulators  appears  to  be  a  better  assessment  of 
the  resolution  in  the  presented  study. 

Similarly,  conductivity  images  collected  from  animals  in 
vivo  showed  significantly  increased  conductivity  in  tumor 
areas.  Tumor  location  was  identified  by  contrast  enhanced 
MRI.  It  is  seen  from  these  results  that  there  is  high  corre¬ 
lation  between  the  conductivity  images  and  contrast 
enhanced  MR  images.  Although  they  do  not  completely 
overlap,  this  is  expected  because  the  two  methods  empha¬ 
size  different  properties  of  tumors.  For  example,  edema 
will  most  likely  show  high  conductivity  but  will  not  show 
signal  enhancement  with  the  CE-MRI.  To  exclude  edema¬ 
tous  regions  from  the  conductivity  images,  one  can  use  T2 
weighted  sequences  that  highlights  edema.  In  this  prelimi¬ 
nary  experiment,  an  animal  with  a  large  tumor  size  was 
imaged,  which  was  available  at  the  time.  As  seen  in  the  CE- 
MRI  images,  the  tumor  had  spread  into  several  foci  that 
encircled  the  animal’s  trunk.  In  large  tumors,  various  com¬ 
partments  like  edema,  necrosis  and  viable  tumor  cells  exist 
and  their  conductivity  and  contrast  agent  enhancement  will 
be  different.  Currently,  a  longitudinal  study  is  planned  to 
observe  changes  in  conductivity  as  well  as  Gd-DTPA  based 
contrast  enhancement  in  tumor  structures  as  the  tumor 
grows.  The  CE-MRI  will  be  used  to  verify  the  MREIT 
results  in  vivo.  At  the  end  of  the  study,  tumors  will  be 
excised  and  undergo  histologic  analysis.  Contrast  enhanced 
and  T2  weighted  images  together  with  the  MREIT  maps 
will  allow  for  the  assessment  of  how  conductivity  correlates 
with  different  compartments  in  tumors.  The  positions  and 
numbers  of  electrodes  could  be  another  factor  that  may  con¬ 
found  the  detection  of  some  low  conductivity  structures 
close  to  surface.  The  effects  of  electrode  placement  will 
also  be  investigated  in  future  studies. 

In  all  phantom  studies,  as  well  as  in  vivo  experiments,  the  cur¬ 
rent  carrying  wires  were  placed  precisely  in  the  z-direction 
with  acrylic  support  beams.  Those  wires  were  attached  to  the 
electrodes  without  any  bends.  Therefore,  the  magnetic  fields 
generated  by  the  currents  in  these  lead  wires  would  have 
mainly  Bx  and  By  components  and  Bz  will  be  negligible. 
Therefore,  any  significant  contribution  from  the  currents  in 
the  lead  wires  to  the  Bz  measurements  is  not  expected. 


In  general,  currents  flowing  in  an  object  generate  a  3D  mag¬ 
netic  field.  The  currents  flowing  inside  an  imaging  slice 
generate  a  magnetic  field  that  has  only  a  component,  which 
is  perpendicular  to  the  plane  of  the  slice.  In  this  preliminary 
study,  the  magnetic  fields  generated  by  currents  outside  the 
imaging  slice  were  ignored,  reducing  the  reconstruction  to 
2D.  Since  most  of  the  current  flow  will  be  concentrated 
between  the  electrodes,  this  approximation  was  deemed  to 
be  adequate  for  the  pilot  studies.  Currently,  3D  reconstruc¬ 
tion  methods  are  being  be  developed  in  our  laboratory  that 
will  account  for  such  out-of- slice  effects. 
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Purpose 

The  difference  between  the  conductivity  of  malignant  tumors  and  the  surrounding  normal  tissue  provides  a  possible  way  to  detect  and  specify  different  tumor  types. 
The  conductivity  distribution  can  be  found  using  magnetic  resonance-electrical  impedance  tomography  (MREIT)  with  a  resolution  higher  than  other  conductivity 
imaging  techniques  [1].  In  MREIT  a  current  distribution  is  created  within  the  object,  resulting  magnetic  flux  density  is  measured  using  MRI,  and  an  inverse  problem  is 
solved  to  find  the  conductivity  distribution  using  magnetic  flux  density  measurements.  MREIT  reconstruction  algorithms  can  be  grouped  into  three:  (1)  current  density 
based  algorithms  (that  calculates  current  density  from  magnetic  field  during  reconstruction),  (2)  gradient  Bz  algorithms  (that  uses  first  or  second  order  derivative  of 
magnetic  flux  density)  and  (3)  sensitivity  matrix  based  algorithms  (that  uses  magnetic  field  directly).  First  group  of  algorithms  requires  rotation  of  the  object  inside  the 
magnet,  which  is  unpractical.  Second  group  of  algorithms  are  very  sensitive  to  noise  in  the  phase  due  to  the  derivative  operator  used.  Sensitivity  matrix  approach  is 
shown  to  perform  better  in  noise;  however,  it  requires  storage  and  processing  of  large  matrices.  Several  results  have  been  presented  for  phantom  and  animal 
experiments  using  2D  sensitivity  approach  [1,2].  In  2D  approximation,  it  is  assumed  that  injected  current  is  confined  to  the  slice  of  interest  and  magnetic  flux  density 
measured  from  the  same  slice  is  used  for  reconstruction.  In  most  experiments,  however,  injected  currents  will  be  distributed  in  all  directions.  Depending  on  the 
conductivity  values  of  the  adjacent  slices,  off-slice  currents  will  contribute  to  the  magnetic  flux  measurement  and  must  be  included  in  modeling.  In  this  study,  we 
conducted  experiments  with  3D  phantoms  and  implemented  a  3D  sensitivity  matrix  reconstruction  algorithm. 


Methods 

When  low  amplitude  sinusoidal  current  is  injected  into  an  object,  the  resulting  magnetic  field  accumulates 
additional  phase  in  the  MR  images.  A  modified  fast  spin-echo  sequence  [3],  where  a  90°  RF  pulse  is  followed  by  a 
number  of  180°  (7t)  RF  pulses,  is  used  to  measure  the  magnetic  field.  Phase  induced  in  positive  and  negative  current 
cycles  would  cancel  out  each  other  unless  multiple  7t  are  used.  By  synchronizing  successive  7t  pulses  to  half  cycles  of 
the  current,  this  phase  shift  accumulates  and  is  given  in  the  final  image  as  (p(r)  =  4-y-N-b(r)  /  co,  where  y  is  the 
gyromagnetic  ratio,  N  the  number  of  cycles  of  injected  current,  b(r)  the  amplitude  of  z-component  current- generated 
magnetic  field  at  point  r,  and  co  the  angular  frequency  of  the  injected  current. 

Once  the  magnetic  field  is  measured,  the  inverse  problem  of  calculating  the  internal  conductivity  distribution  from 
these  measurements  is  solved.  Sensitivity  based  reconstruction  algorithm  is  implemented  in  3D.  Uniform  conductivity 
distribution  is  assumed  and  sensitivity  matrix  (S)  is  calculated  using  a  3D  finite  element  method  (FEM).  Sensitivity 
matrix  gives  the  relation  between  Ab  (the  difference  between  measured  magnetic  flux  density  and  the  magnetic  flux 
density  corresponding  to  initial  distribution)  and  Aa  (the  change  with  respect  to  initial)  as  Ab  =  SAo.  Including 
Tikhonov  regularization  parameter,  X,  the  matrix  equation  becomes  (SrS  +^I)  Aa  =  STAb  where  I  is  the  identity 
matrix.  The  matrix  equation  is  solved  for  different  values  of  X  using  conjugate  gradient  method  and  the  optimum 
regularization  value  is  selected  as  the  one  minimizing  the  difference  between  measured  and  calculated  magnetic  flux 
densities. 
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Fig  1.  Pulse  Sequence 


Results 

A  3D  phantom  was  prepared  using  agarose  powder  (lgr/lOOmL)  as  background  and  a  spherical  shell  filled  with  the  same  agarose  gel.  Due  to  insulator  shell,  the 
whole  sphere  acted  as  a  nonconducting  object.  Four  electrodes  of  size  0.5cm  x  0.5cm  were  placed  at  the  center  of  the  cylinder  every  90°  as  shown  in  Fig  2.  Opposite 
pairs  were  used,  giving  two  current  injection  profiles.  A  FEM  mesh  with  2399  nodes  and  11781  tetrahedral  elements  was  constructed.  The  magnetic  flux  density  is 
calculated  at  a  rectangular  grid  of  24x24x5,  where  2300  of  grid  points  falling  into  the  cylindrical  region.  A  total  of  4600  measurements  were  used  in  image 
reconstruction.  MR  structural  images  for  the  5  slices  used  are  given  in  Fig.  3  (top  row).  Corresponding  reconstructed  conductivity  images  are  presented  in  bottom  row, 
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Fig  3.  MR  anatomical  images  (top  row)  and  corresponding  reconstructed  conductivity  images  (bottom  row) 


Discussion 

In  this  study,  we  demonstrated  the  feasibility  of  using  3D  sensitivity  matrix  reconstruction  in  multi-slice  magnetic  resonance  electrical  impedance  imaging. 
Although  sensitivity  matrix  based  reconstruction  methods  are  computationally  expensive,  they  perform  well  in  the  presence  of  noise.  Previously,  we  used  sensitivity 
matrix  based  reconstruction  using  2D  FEM.  That  approach  ignores  the  contributions  of  currents  in  off-center  slices  leading  to  various  artifacts  in  the  conductivity 
images;  conductivity  values  in  the  reconstructed  slice  could  be  under  or  overestimated,  spatial  resolution  and  point  spread  function  could  be  compromised.  In  some 
cases  the  ‘ghosts’  of  objects  in  other  slices  may  even  appear  in  the  reconstructed  slice  depending  on  how  the  conductivity  and  resulting  3D  current  density  is  distributed 
in  the  object. 
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Purpose 

Low  level  injected  currents  generate  a  magnetic  field  within  a  body  that  can  be  measured  using  MRI.  Magnetic  resonance-electrical  impedance  tomography 
(MREIT)  utilizes  these  measurements  due  to  injected  currents  to  construct  the  conductivity  distribution  within  an  object.  Electrical  conductivities  of  tissues  are  different 
among  normal,  malignant,  and  benign  tissues.  In-vivo  impedance  imaging  of  suspicious  lesions  could  aid  in  the  diagnosis  and  specification  of  malignant  tumors. 
Several  studies  using  phantoms  and  small  animal  bearing  tumors  demonstrated  the  efficacy  of  the  algorithm  [1,  2].  The  magnetic  flux  density  is  extracted  from  MRI 
phase  images  and  one  of  the  challenges  in  MREIT  is  acquiring  phase  information  with  sufficient  SNR  while  maintaining  the  injected  currents  at  safe  levels.  Sadleir  et. 
al.  carried  out  noise  analysis  in  magnetic  flux  density  measurements  at  3T  and  11T  and  they  suggest  that  the  noise  level  in  magnetic  flux  density  is  reduced  as  field 
strength  is  increased  by  a  factor  approximately  proportional  to  the  increase  in  the  field  strength  [3].  In  this  study,  we  carried  out  experiments  at  4T  and  7T  field 
strengths  to  compare  the  effect  of  increased  magnetic  flux  density  SNR  in  the  final  reconstructed  conductivity  images. 
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Methods 

For  the  test  phantom,  a  hollow  acrylic  tube  and  agarose  gels  with  different 
NaCl  concentrations  were  used.  2gr/100mL  agarose  powder  and  4mM  CuS04  were 
kept  the  same  for  all  regions.  Amount  of  NaCl  was  set  to  lgr/lOOmL  and 
4gr/100mL  for  background  and  objects  (01,02),  respectively.  The  conductivity 
values  were  measured  in  an  independent  system  as  S.bBSm'1  and  l.blSm'1  (3.3:1 
contrast).  The  diameter  of  the  both  objects  were  6mm  and  their  centers  were  located 
lcm  apart  giving  4mm  gap  between  their  boundaries.  The  thickness  of  the  phantom 
was  lcm  and  in-slice  dimensions  of  the  phantom  are  defined  in  Fig  1.  Four  copper 
electrodes  each  3  mm  wide  were  placed  equidistantly  along  the  inner  acrylic  wall 
and  used  to  probe  the  interior  region.  Bipolar  current  injection  scheme  was  used  the 
resulting  magnetic  flux  density  distribution  measured  using  a  modified  spin-echo 
pulse  sequence  as  shown  in  Fig  2.  [4].  Opposite  electrode  pairs  were  used  to  get 
two  current  injection  profiles.  Iterated  sensitivity  matrix  method  (SMM)  with 
Tikhonov  regularization  was  used  for  conductivity  reconstruction  [1]. 

Results 

MREIT  experiments  were  carried  out  using  two  phantoms  with  the  same 
distribution  in  4T  and  7T  systems.  Since  ion  diffusion  alters  the  initial  conductivity 
distribution  [5] ,  we  prepared  the  phantom  immediately  before  each  experiment  to 
minimize  these  effects  in  the  comparison.  The  pulse  sequence  parameters  for  both 
experiments  were  as:  TR  =  500ms,  TE  =  60ms  FOV  =  7cm,  slice  thickness  =  2mm, 
field  image  matrix  =  64x64,  NEX  =  4.  The  sampling  bandwidth  was  set  to  33.3KHz 
and  1.25MHz  for  4T  and  7T  systems,  respectively.  The  amplitude  of  the  injected 
current  was  1mA  and  the  duration  of  the  positive  and  negative  cycles  were  27.3ms 
and  27.7ms.  Reconstructed  conductivity  images  and  profiles  along  p-line  (as 
defined  in  Fig  1)  at  both  field  strengths  are  given  in  Fig  3.  Although  the  phantom  is 
slightly  rotated  between  two  cases,  this  does  not  affect  the  relative  positions  of  the 
objects  with  respect  to  each  other.  We  calculated  full-width-at-half-maximum 
(FWHM)  for  both  objects  through  x  and  y  cross-sections  and  also  we  looked  at 
peak  contrast  in  the  object.  Note  that  the  SMM  reconstructs  relative  conductivity 
values  if  magnetic  flux  density  measurements  are  used  alone  and  the  values 
presented  here  are  normalized  to  the  background  conductivity.  The  peak 
reconstructed  contrast  and  FWHM  values  are  summarized  in  Table  1  for  both  field 
strengths  and  objects.  The  objects  were  not  fully  resolved  in  the  p-line  profile  for 
4T  case  but  more  separable  in  the  7T  case.  For  all  cases,  the  constructed  contrasts 
are  lower  than  the  expected  ones,  especially  for  the  object  located  at  the  center. 

Table  1.  Comparison  of  reconstructed  peak  contrast  and  FWHM 


peak  contrast 

FWHM-x  (cm) 

FWHM-y  (cm) 

Ol 

02 

Ol 

02 

Ol 

02 

True 

3.3 

3.3 

0.6 

0.6 

0.6 

0.6 

4T 

1.6 

1.8 

1.3 

1.0 

1.4 

0.9 

7T 

1.8 

2.2 

1.0 

0.8 

1.1 

0.8 

RF-O— 

Gz-^tr 


n 


AC- 


Fig  1.  Phantom  definitions 


Fig  2.  MREIT  Pulse  Sequence 


(a) 


(b)  Position  (cm) 
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(d)  Position  (cm) 


Fig  3. (a)  image  at  4T  (b)  4T  profile  (along  p-line) 

Discussion  (C)  imag£  a‘  ?T  (d)  7T  Pr°flle  (al°ng  P'Hne) 

In  this  study,  we  carried  out  MREIT  experiments  at  two  different  field  strengths  to  understand  how  the  improvement  in  magnetic  flux  density  SNR  affects  the 
reconstructed  conductivity  values  and  target  detection.  When  the  reconstructed  images  are  compared,  it  is  seen  that  given  the  same  amount  of  experimentation  time,  the 
accuracy  and  the  spatial  resolution  of  the  conductivity  images  at  7T  is  better  compared  to  4T.  Although  the  reconstructed  contrast  values  are  lower  than  the  expected 
values  due  to  ion  diffusion  in  all  cases,  the  comparison  between  different  field  strengths  shows  13%-22%  improvement  at  7T  compared  to  the  4T  case.  Major 
improvement  is  noticable  in  the  resolution  of  the  objects,  especially  for  the  one  at  the  center,  when  higher  field  strength  is  used. 
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Purpose 

In  Magnetic  Resonance  Electrical  Impedance  Tomography  (MREIT),  electrical  currents  are  injected  into  an  object  and  the  resulting  magnetic  flux 
density  distribution  measured  using  MRI.  Various  spin  echo  based  pulse  sequences  have  been  proposed  to  perform  this  measurement,  where  a  single 
line  of  k- space  data  is  collected  per  excitation.  In  this  study,  we  investigate  the  use  of  a  single  shot,  spin  echo,  echo  planar  imaging  (SS-SEPI)  pulse 
sequence  for  acquiring  MREIT  data. 

Methods 

For  the  test  phantom,  a  hollow  acrylic  disk  with  an  inner  diameter  of  7cm  and  thickness  of  1cm  was  filled  with 
2%  agarose  and  4mM  CuS04.  Within  this  disk,  a  smaller  cylindrical  region  of  13mm  diameter  was  filled  with  1% 

NaCl,  2%  agarose,  and  4mA  CuS04  to  generate  a  high  conductivity  region  (Figure  1).  The  plane  of  the  disk  was 
placed  perpendicular  to  the  main  static  MRI  field.  Two  copper  electrodes  each  6mm  wide  were  placed  opposite  of 
each  other  along  the  inner  acrylic  wall  and  used  to  inject  currents  into  the  interior  region. 

A  4mA  bipolar  current  pulse  was  injected  into  the  phantom  and  the  resulting  magnetic  flux  density  distribution 

measured  using  a  modified  SS-SEPI  pulse  sequence 
(Figure  2).  The  component  of  current- generated 
magnetic  flux  density  parallel  to  the  main  static 
MRI  field  (z-component)  introduced  a  phase  shift  in 
the  MR  image  given  as:  cp(r)  =  y-(TA+TB)-Bz(r). 

The  scan  parameters  were:  TA  =  35ms,  TB  =  10ms,  TR  =  3ms,  TE  =  60ms,  matrix  = 
64X64,  FOV  =  20cm,  and  single  slice  thickness  =  5mm.  Data  was  collected  twice, 
each  with  opposite  polarities  in  the  applied  current  waveform.  The  resulting  phase 
maps  were  subtracted  then  divided  by  two,  so  as  to  cancel  out  any  additional  phase 
contributions,  such  as  those  arising  from  small  imperfections  in  the  hardware  timing. 

To  reconstruct  the  conductivity  distribution  using  the  MRI  measurements,  the 
Sensitivity  Matrix  Method  was  utilized  [Birgul  et  al,  Phys  Med  Biol  51:  5035-5049 
(2006)]  in  which  the  relationship  between  conductivity  and  z-component  magnetic 
flux  density  is  linearized  around  an  initial  conductivity  (i.e.  uniform  distribution) 
and  formulated  as  a  matrix  equation.  This  equation  is  then  solved  for  the  true 
conductivity  distribution  using  Tikhonov  regularization.  The  resulting  conductivity 
can  then  be  substituted  back  into  the  linearized  equation  as  the  new,  updated  initial 
condition,  and  the  process  iterated  to  improve  the  reconstruction. 
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Fig.  2.  SS-SEPI  pulse  sequence  for  MREIT 


Fig.  1.  Phantom  Schematic 


Results 

Data  was  collected  using  a  4T  MRI  system.  For  comparison,  data  was  also  acquired  using  no  injected  current,  and  also  using  the  original  SE- 
based  measurement  sequence  [Scott  et  al,  IEEE  TMI 10:  362-374  (1991)].  Relative  conductivities  were  reconstructed  for  the  various  data  using  5 
iterations  of  the  Sensitivity  Matrix  Method  (Figure  3).  As  part  of  the  reconstruction  process,  a  Finite  Element  Method  mesh  was  aligned  to  the  MR 
image,  and  the  resulting  conductivity  map  linearly  transformed  back  into  a  circular  geometry. 


Discussion 

Comparison  of  the  0mA  and  4mA  cases  shows  that  the  SS-SEPI  pulse  sequence  was  able  to  detect  the  magnetic  flux  density  generated  by  the 
injected  current.  Comparison  with  the  SE  results  shows  that  the  SS-SEPI  pulse  sequence  was  also  able  to  map  out  the  general  variation  in  the 

magnetic  flux  density  due  to  the  high  conductivity 
perturbation.  However,  the  SS-SEPI  images  suffer  from 
geometric  distortions  resulting  from  the  low  sampling 
bandwidth  in  the  phase  encode  direction  inherent  in  EPI- 
based  pulse  sequences.  As  a  result,  the  reconstructed 
conductivity  distribution  also  suffers  from  some  geometric 
distortion.  Never  the  less,  the  high  conductivity  perturbation 
could  still  be  reconstructed  from  the  SS-SEPI  data.  We  could 
further  improve  reconstruction  by  applying  some  sort  of 
geometric  correction  to  the  MRI  data. 

An  advantage  of  the  SS-SEPI  pulse  sequence  is  the 
significant  reduction  in  data  acquisition  time.  This  could 
allow  for  increased  signal  averaging  in  an  allotted  study  time, 
which  may  be  required  for  improving  the  SNR  when  using 
low  amplitude  injected  currents.  Single  shot  pulse  sequences 
are  also  less  susceptible  to  motion  artifacts,  which  may  be  of 
significance  when  applying  MREIT  to  in  vivo  studies. 
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Fig.  3.  Magnetic  flux  density  and  reconstructed  conductivity  distributions 
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Purpose 

In  Magnetic  Resonance  Electrical  Impedance  Tomography  (MREIT),  electrical  currents  are  injected  into  an  object  and  the  resulting  magnetic  flux 
density  distribution  measured  using  MRI.  These  MRI  measurements  are  then  used  to  reconstruct  the  conductivity  distribution  within  the  object.  We 
have  recently  reported  on  the  measurement  of  ion  diffusion  with  MREIT  using  injected  current  pulses  of  900jlA  [Hamamura  et  al,  Phys  Med  Biol 
51:  2753-2762  (2006)].  To  our  knowledge,  this  is  the  lowest  level  of  current  successfully  used  in  MREIT  to  date.  However,  this  level  still  remains 
above  the  required  safety  limits  for  application  in  human  studies.  The  IEC  601  standard  limits  “patient  auxiliary  currents”  to  lOOpA  at  low 
frequencies.  In  this  study,  we  assess  the  efficacy  of  MREIT  using  such  biologically  safe  injected  current  levels. 


Fig.  1.  Phantom  Schematic 


Methods 

Lower  injected  current  levels  reduce  the  SNR  of  the  acquired  magnetic  flux  density  maps  used  in  the  reconstruction 
process.  To  improve  the  SNR,  we  utilized  a  7T  high  field  system  with  increased  signal  averaging  to  acquire  our  MRI 
measurements.  For  the  test  phantom,  a  hollow  acrylic  disk  with  an  inner  diameter  of  4.445cm  and  thickness  of  1cm 
was  filled  with  2%  agarose  and  4mM  CuS04.  Within  this  disk,  a  smaller  plastic  shell  of  6mm  diameter  was  placed 
slightly  off-center  to  simulate  a  low  conductivity  region  (Fig.  1).  The  plane  of  the  disk  was  placed  perpendicular  to  the 
main  static  MRI  field.  Four  copper  electrodes  each  3mm  wide  were  placed  equidistant  along  the  inner  acrylic  wall  and 
used  to  inject  currents  into  the  interior  region. 

A  bipolar  current  pulse  was  injected  into  the  phantom  and  the  resulting  magnetic  flux  density  distribution  measured 
using  a  modified  spin-echo  pulse  sequence  (Fig.  2)  [Scott  et  al,  IEEE  TMI 10:  362-374  (1991)].  The  scan  parameters 
were:  TR  =  500ms,  TE  =  60ms,  matrix  =  128X128,  FOV  =  14cm,  and  single  slice  thickness  =  1mm.  Data  was 

collected  for  two  different  current  injection  schemes  (in  pairs  of  electrodes  directly 
opposite  of  each  other)  and  used  simultaneously  in  conductivity  reconstruction.  To 
reconstruct  the  conductivity  distribution  using  the  MRI  measurements,  the  Sensitivity 
Matrix  Method  was  utilized  [Birgul  et  al,  Phys  Med  Biol  51:  5035-5049  (2006)]  in 
which  the  relationship  between  conductivity  and  magnetic  flux  density  is  linearized 
around  an  initial  conductivity  (i.e.  uniform  distribution)  and  formulated  as  a  matrix 
equation.  This  equation  is  then  solved  for  the  true  conductivity  distribution  using 
Tikhonov  regularization.  The  resulting  conductivity  can  then  be  substituted  back  into 
the  linearized  equation  as  the  new, 
updated  initial  condition,  and  the 
process  iterated  to  improve  the 
reconstruction. 
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Results 

Data  was  collected  with  various 
amounts  of  signal  averaging  using  a 
current  amplitude  of  lOOpA.  For 
comparison,  high  SNR  data  was 
collected  by  using  5mA  of  injected 
current.  Relative  conductivities 
were  reconstructed  for  the  various 
data  using  5  iterations  of  the 
Sensitivity  Matrix  Method  (Fig.  3). 


-0>- 


ACQ 


Fig.  2.  MREIT  Pulse  Sequence 


Discussion 

The  results  of  this  study  demonstrate  that  MREIT  is  capable  of  reconstructing 
conductivities  using  biologically  safe  injected  current  levels.  Even  with  a  single 
average,  we  can  observe  the  lower  conductivity  perturbation  within  our  test  phantom. 
However,  lower  SNR  data  generated  considerable  variation  in  the  (uniform) 
background.  Improving  the  SNR  through  increased  signal  averaging  reduced  these 
artifacts. 

When  applying  MREIT  to  human  studies,  the  larger  object  size  and  3D  current  flow 
will  result  in  a  decrease  in  the  measurable  magnetic  flux  density.  This  study  illustrates 
that  increased  signal  averaging  can  compensate  for  this  corresponding  decrease  in 
SNR.  In  practice,  the  presence  of  correlated  noise  decreases  the  efficiency  of 
additional  averaging,  and  places  an  upper  bound  on  SNR  gain.  This  limitation  must  be 
investigated  further,  as  well  as  other  techniques  to  improve  SNR  when  using  low 
amplitude  injected  currents. 
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Fig.  3.  Magnetic  flux  density  and  reconstructed  conductivity 
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Introduction: 

As  an  adjunct  to  existing  cancer  detection  techniques  such  as  x-ray  and  MRI,  Magnetic  Resonance-Electrical  Impedance  Tomography  (MREIT)  has 
the  potential  to  improve  the  specificity  of  current  techniques  because  it  is  an  imaging  modality  that  measures  tissue  conductivity  and  it  has  been  shown 
that  the  conductivity  differs  significantly  among  malignant,  benign  and  normal  tissues.  A  sensitivity  based  reconstruction  algorithm  with  Tikhonov 
regularization  has  been  successfully  employed  to  resolve  a  complex  and  nested  conductivity  distributions  [1].  In  this  work  a  mutual  information  based 
MREIT  reconstruction  algorithm  was  investigated  to  reconstruct  the  irregular  conductivity  distributions  with  improved  accuracy. 

Methods: 

The  injected  currents  inside  an  object  will  generate  magnetic  flux  density;  the  component  of  the  magnetic  flux  density  in  the  main  static  field 
direction  (bz)  can  be  calculated  from  the  MRI  phase  images  using  a  modified  spin  echo  pulse  sequence  [2] .  To  solve  the  conductivity  distribution  from 
the  magnetic  field  measurement,  the  relationship  between  the  perturbation  of  conductivity  distribution  (Ao)  around  an  initial  estimate  (Gj)  and  the 
resultant  perturbation  of  the  magnetic  flux  density  (Abz)  can  be  linearized  (Abz=SAo)  for  a  given  current  injection  scheme.  Here  S  is  the  sensitivity 
matrix  that  can  be  calculated  analytically  [3].  Due  to  the  ill-posed  nature  of  the  sensitivity  matrix  S  and  noise  in  the  data,  regularization  is  required  to 
solve  the  inverse  problem  of  calculating  Ao  from  Abzand  S.  Including  the  Tikhonov  regularization  parameter  X,  the  matrix  equations  becomes  (STS+  /J) 
Ao  =  STAbz  where  I  is  the  identity  matrix,  and  the  equation  can  then  be  solved  using  the  conjugate  gradient  method.  Previously  X  was  chosen  to  be  the 
value  that  minimizing  the  difference  between  the  measured  magnetic  flux  density  and  the  magnetic  flux  density  calculated  from  the  reconstructed 
conductivity.  In  this  study  we  investigated  mutual  information  based  reconstruction  methods  [4]  for  MREIT.  In  the  first  approach,  the  X  value  that 
maximized  the  mutual  information  between  the  calculated  conductivity  distribution  image  and  the  MR  image  was  found.  This  is  valid  because  mutual 
information  is  the  measurement  of  nonlinear  statistical  dependency  of  two  systems,  and  since  both  the  MR  image  and  conductivity  image  contain 
structural  information  about  the  tumor  and  the  normal  tissue,  the  mutual  information  between  the  two  images  should  be  maximized  when  the 
conductivity  image  is  reconstructed  most  accurately. 

Since  the  location  and  shape  of  the  suspected  lesions  can  be  detected  by  MRI,  one  can  generate  an  initial  conductivity  distribution  with  a 
two -compartment  model  that  separates  the  suspected  lesion  and  the  normal  tissue.  This  model  will  be  a  better  approximation  of  the  initial  conductivity 
distribution  (gO  instead  of  the  standard  approach  of  using  a  uniform  Gj.  For  the  second  approach,  we  first  solve  the  inverse  problem  to  find  the 
conductivity  values  in  these  two  compartments  only,  so  a  single  conductivity  value  is  calculated  for  each  compartment.  The  resulting  two-compartment 
conductivity  map  is  used  as  the  initial  conductivity  distribution  in  the  second  step  and  the  equation  above  is  solved  to  find  the  conductivity  value  for 
each  pixel.  Since  this  two-compartment  Gi  is  a  better  approximation,  starting  the  iterative  reconstruction  from  an  initial  condition  closer  to  the  actual 
values  should  yield  more  accurate  conductivity  maps  and  converge  to  the  solution  faster. 

Results  and  discussion: 

The  same  data  acquired  in  the  experiment  described  in  [1]  were  used  in  the  new  reconstruction  method.  The  reconstructed  conductivity  distribution 
using  the  original  method  in  [1],  the  mutual  information  based  reconstruction  method  with  optimum  X  (uniform  initial  conductivity),  and  the  mutual 
information  based  reconstruction  method  with  two -compartment  initial  conductivity  distribution  are  shown  in  Fig.l.  The  two  compartments  were 
defined  from  the  MR  images  that  provide  the  contrast  between  compartments  I  and  II  in  Fig .2.  The  reconstructed  conductivities  in  the  three  cases  were 
similar  in  general  except  in  the  regions  marked  by  the  red  circles  where  substantial  improvement  was  seen.  Moreover,  calculated  conductivities  from 
the  mutual  information  based  reconstructions  were  much  closer  to  the  actual  values  compared  to  the  one  calculated  from  the  original  method,  where 
conductivity  values  turned  out  higher.  Also  note  that  the  third  compartment  was  successfully  recovered  although  the  Gi  was  modeled  with  only  two 
compartments.  Therefore,  it  can  be  concluded  that  maximization  of  mutual  information  can  be  used  as  a  criterion  to  reduce  the  uncertainty  (such  as 
estimation  of  the  regularization  parameter)  in  the  MREIT  conductivity  reconstruction,  thus  the  accuracy  of  the  conductivity  measurement  can  be 
improved  by  using  the  MR  images  as  a  priori  information. 


Figl.  The  reconstructed  conductivity  distribution  using  the  original  method  a),  the  mutual  information  based 
reconstruction  method  with  optimum  X  and  uniform  initial  conductivity  distribution  b),  and  the  mutual  information 
based  reconstruction  method  with  two  compartment  averaged  initial  conductivity  distribution  c). 


Fig  2.  The  schematic  of  the  phantom 
with  three  nested  conductivity 
compartments  labeled  I,  II  and  III 
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Introduction :  The  goal  of  this  study  was  to  investigate  the  effects  of  limited  volume  coverage  on  the  accuracy  of  conductivity  reconstruction  in 
MR-Electrical  Impedance  Tomography  (MREIT).  Since  the  currents  injected  into  an  object  will  be  distributed  in  the  whole  volume,  limited  volume 
coverage  in  data  acquisition  or  reconstruction  will  impact  the  accuracy  of  resulting  conductivity  maps.  On  the  other  hand,  the  magnitude  of  current 
density  will  decrease  rapidly  as  the  distance  from  the  injecting  electrodes  increases.  Moreover,  if  one  is  interested  in  a  particular  volume  of  interest 
(VOI),  the  magnetic  fields  generated  by  currents  away  from  the  VOI  could  be  negligible.  Therefore,  we  investigated  the  decrease  in  magnitude  of 
current  density  as  the  distance  from  the  electrodes  increases.  We  also  investigated  the  contribution  of  these  weaker  currents  to  the  magnitude  of  the 
magnetic  fields  generated  within  the  VOI. 

Methods :  In  MREIT,  weak  electrical  currents  are  injected  into  an  object  that  generate  magnetic  fields,  the  z- 
component  of  which  induces  additional  phase  information  in  MR  images.  A  modified  spin-echo  sequence  was 
used  with  several  7t  pulses  applied  during  the  zero- crossings  of  an  alternating  current  and  a  phase  shift 
accumulates,  which  is  given  in  the  final  image  as  cp(r)  =  4  y-N  bz(r)/co,  (y:  gyromagnetic  ratio;  N:  #  cycles  of 
injected  current;  bz(r):  the  current- generated  magnetic  field  at  point  r;  co:  angular  frequency  of  the  injected 
current).  Once  bz(r)  is  calculated  from  the  phase  (p(r)  measurements,  one  can  devise  a  method  to  calculate  the  E  z 
conductivity  map  [Muftuler  LT  et  al  TCRT  v  3,  599-610,  2004].  Since  only  the  z-component  of  the  magnetic  field  is 
used  in  MREIT,  the  transverse  components  of  the  current  density,  Jx(r)  and  Jy(r)  are  of  interest.  We  performed  several 
simulation  studies  to  calculate  the  current  density  and  magnetic  fields  inside  a  cylindrical  volume  (d=4.5cm  h=7cm) 
with  uniform  conductivity  distribution.  Two  electrodes  were  placed  at  135°  and  315°  around  the  central  transaxial  slice 

(slice  0).  First,  the  electric  potential  distribution  was  determined  by  solving  the  Poisson’s  equation  with  Neumann  boundary  value  problem  (Eq.l) 
using  Finite  Element  Method  (FEM).  Here,  <7 is  the  3D  conductivity  distribution  and  (f>  is  the  electric  potential.  The  current  density  and  the  magnetic 
flux  density  are  calculated  using  equations  2  and  3,  respectively.  In  3D  FEM  model  13773  nodes  and  73452  of  tetrahedral  elements  were  used. 

Results:  Fig. la  shows  the  magnitude  of  the  transverse  current 
density  (Zx(r)2+  Jy{ r)2)1/2  across  15  slices  for  regions  labeled  13-19 
in  Fig. lb.  bz(r)  was  calculated  for  five  cases  on  these  15  transverse 
planes  that  were  5mm  apart.  For  each  case,  the  current  density  only 
within  a  transverse  slab  of  various  thicknesses  was  used  for  bz(r) 
calculations:  easel :  4.2cm,  case2:  3cm,  case3 :  2.4cm,  case4. 

1.8cm  and  case5 :  1.2cm.  Fig.2  shows  the  bz(r)  contour  maps  in 
slice  0.  Fig.  2a  is  the  resulting  magnetic  field  when  the  transverse 
current  density  in  4.2cm  thick  slab  was  used  {easel).  Similarly, 

Fig. 2b  is  the  bz(r)  for  case5.  Fig. 2c  is  the  difference  field  when  2b 
is  subtracted  from  2a.  When  easel  is  taken  as  the  reference,  the 
field  difference  maps  (Abz(r))  between  easel  and  others  give  us  a 
measure  of  the  errors  made  in  the  calculation  of  bz(r)  when  the 
current  density  outside  the  selected  slab  is  ignored.  Therefore,  we 
defined  a  region  of  interest  (ROI)  that  encompassed  the  white 
contour  band  in  the  upper  right  hand  corner  of 
Fig. 2c.  The  mean  bz(r)  inside  the  ROI  for 
easel  is  taken  as  the  reference  and  the  mean  of 
Abz(r)  for  each  case  is  calculated  and  divided 
by  this  reference  mean.  Results  are  summarized 
in  table  1 . 

Discussion.  The  result  in  Fig  la  shows  that 
the  current  density  in  the  plane  that  is  21mm 
away  from  the  electrodes  has  significant 
magnitude,  (e.g.  in  region  19,  the  magnitude  of 
the  current  is  roughly  40%  of  that  of  slice  0). 

Moreover,  the  magnetic  field  in  the  ROI  in 

slice  0  is  underestimated  by  almost  30%  if  only  currents  within  a  12mm  slab 
are  taken  into  account.  This  will  lead  to  underestimation  of  conductivity 
maps  in  this  slice.  Therefore,  one  has  to  acquire  MREIT  data  that  covers 
sufficiently  large  volume  to  obtain  accurate  conductivity  maps.  On  the  other 
hand,  reconstruction  of  conductivity  maps  from  large  data  sets  may  require 
18-20  hours  of  computing  time.  Therefore,  optimum  volume  coverage  has  to 
be  found  to  obtain  a  balance  between  acceptable  accuracy  and  computational  efficiency  by  carrying  out  simulations  and  experiments.  The 
preliminary  results  presented  here  were  obtained  from  a  uniform  conductivity  phantom.  The  results  from  objects  of  nonuniform  conductivity  will  be 
different  and  models  should  be  developed  accordingly  to  find  the  optimum  volume  coverage  for  studies  such  as  in  vivo  experiments. 
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Fig.l.  (a)  magnitude  of  current  density  across  slices  in  selected  regions  shown  in  the 
field  map  (b). 


Fig.2.  bz(r)  contour  maps  for  case  1  (a),  case  5  (b)  are  illustrated.  The  difference  of  the  two  bz(r)  maps  is 
shown  in  (c).  Current  was  injected  between  electrodes  placed  at  135°  and  315°. 
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Table  1 .  Mean  of  Abz(r)  for  each  case  divided  by  the  reference  mean 
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Purpose 

It  is  possible  to  detect  locations  of  lesions  in  breast  cancer  using  techniques  such  as  x-ray  mammography  or  MRI  accurately,  however,  the  specificity  of  current 
techniques  is  low  [1].  Since  the  conductivity  values  of  malignant,  benign,  and  normal  tissues  are  significantly  different,  this  information  can  be  used  in  classification  to 
improve  specificity.  Magnetic  resonance-electrical  impedance  tomography  (MREIT)  is  an  imaging  modality  that  reconstructs  conductivity  images  from  magnetic  field 
measurements  generated  due  to  a  current  distribution  in  a  volume  conductor.  In  MREIT  low  amplitude  sinusoidal  current  is  injected  into  an  object  and  the  resulting 
magnetic  field  accumulated  additional  phase  in  MR  images.  A  modified  fast  spin-echo  sequence  is  used  to  measure  this  magnetic  field.  These  measurements  are  used  to 
solve  the  inverse  problem  of  finding  the  conductivity  distribution  inside  the  object  using  an  iterated  sensitivity  reconstruction  algorithm.  In  most  cases,  the  conductivity 
distribution  is  expected  to  vary  within  the  tumor.  There  are  several  phantom  studies  in  literature  that  assess  the  performance  of  MREIT  using  simple  cases  but  none 
investigates  whether  it  is  capable  of  detecting  complex  and  nested  conductivity  distributions,  which  models  the  real  life  cases  more  accurately.  In  this  study,  we  used  a 
new  criterion  for  selection  of  optimum  regularization  parameter  in  image  reconstruction  and  showed  that  it  is  possible  to  resolve  a  6mm  inhomogeneity  within  an 
irregular  region  using  1mA  peak  current. 

Methods 

Reconstruction  of  conductivity  involves  two  basic  steps.  The  first  step  is  the  measurement  of  magnetic  flux  density  using  magnetic  resonance  imaging.  This  step 
involves  MRI  data  acquisition  using  a  modified  spin  echo  pulse  sequence  [2]  and  generation  of  magnetic  flux  density  images  from  MRI  phase  images  using  scaling. 
The  component  of  the  magnetic  flux  density  in  the  direction  of  the  main  static  field  of  MRI  is  sufficient  in  image  reconstruction.  In  the  second  step,  these  images  are 
used  as  input  data  in  the  inverse  problem  of  finding  conductivity  from  magnetic  field  information.  Sensitivity  based  reconstruction  algorithm  is  implemented  for  the 
solution  of  the  inverse  problem.  Uniform  conductivity  distribution  is  assumed  and  sensitivity  matrix  is  calculated  analytically  [3].  Resulting  matrix  equation  is  given  as 
Ab  =  SAg  where  Ab  is  the  difference  between  measured  magnetic  flux  density  and  the  magnetic  flux  density  corresponding  to  initial  distribution,  Ag  is  the  change  with 

respect  to  initial  and  S  is  the  sensitivity  matrix  that  gives  the  relation  between  changes  in  magnetic  field  and  conductivity.  Including  Tikhonov  regularization  parameter, 

T  T 

X,  the  matrix  equation  becomes  (S  S  +AI)  Ag  =  S  Ab  where  I  is  the  identity  matrix.  The  matrix  equation  is  solved  for  different  values  of  X  using  conjugate  gradient 
method  and  the  optimum  regularization  value  is  selected  as  the  one  minimizing  the  difference 

minEK««  ~BcalcMi 

where  m  is  the  total  number  of  measurement  points,  Bmeas  is  the  measured  magnetic  flux  density,  and  Bcaic  is  the  flux  density  calculated  using  reconstructed 
conductivity.  Calculated  conductivity  distribution  is  assigned  as  the  initial  value  and  the  steps  starting  with  sensitivity  matrix  calculation  are  repeated  until  the  change  in 
conductivity  two  consecutive  iterations  are  below  a  defined  threshold. 


Results 

A  three  region  conductivity  phantom  is  generated  using  agarose  and  different  concentrations  of  NaCl  (Fig  la).  For  the  background,  1%  (g/lOOmL)  agarose  and  1% 
NaCl  is  used.  Agarose  amount  is  kept  constant  in  all  regions  but  NaCl  is  increased  to  2%  and  4%  for  regions  II  and  III  respectively.  Corresponding  conductivity  values 
are  measured  using  a  4  electrode  conductivity  cell  and  the  true  values  are  found  to  be  1.54,  2.94,  and  5.92mS/cm.  Four  electrodes  are  placed  around  the  object  giving  6 
different  current  injection  profiles.  First,  high-resolution  anatomical  images  are  acquired  to  determine  the  exact  locations  of  the  electrodes  to  be  used  in  the  image 
reconstruction  (Fig  lb).  Since  the  water  content  is  same  in  all  regions,  the  irregular  region  and  the  6mm  inclusion  cannot  be  seen  in  the  MR  image.  Then  4  cycles  of 
1mA  (peak)  100Hz  current  were  injected  with  TR  =  1000ms,  TE  =  50ms,  and  NEX  =  4  and  field  is  measured  for  6  cases  (Fig  lc).  Measurements  around  the  outer  ring 
are  masked  out  to  eliminate  spurious  boundary  effects  in  images  reconstruction.  Reconstructed  conductivity  image  is  given  in  Fig  Id.  where  the  6mm  region  inside  the 
irregular  image  is  resolved.  Peak  reconstructed  values  for  each  region  are  1.23,  2.55,  and  4.93mS/cm  for  each  region. 


Figure  la 


Figure  lb 


Figure  Id 


Figure  1  (a)  Schematic  of  the  multi- compartment  conductivity  phantom  (b)  MR  magnitude  images  (c)  Scaled  and  masked  MR  phase  images  for  6  different  current 
injection  case  (d)  reconstructed  conductivity  image  (red  indicating  higher  conductivity  and  blue  indicating  lower  conductivity 


Discussion 

In  this  study,  we  have  shown  that  it  is  possible  to  reconstruct  complex  conductivity  distributions  within  an  object  using  MREIT  technique  and  iterated  sensitivity 
reconstruction  algorithm.  The  irregular  object  and  6mm  inclusion,  which  cannot  be  seen  in  the  anatomical  MR  images,  are  clearly  resolved  in  the  conductivity  image 
with  16.7%  peak  error  at  current  levels  which  are  acceptable  for  human  imaging  applications. 

References 

[1]  Malich  A,  et.  al,  Eur.  Radiol,  10:  1555-1561  (2000),  [2]  Muftuler  L  T,  et.  al., TCRT December  2004  [3]  Birgul  O,  et.  al.  PMB,  48:  3485-2504  (2003) 
Acknowledgments 

This  research  is  supported  in  part  by  Department  of  Defense  Award  W81XWH-04- 1-0446  and  NIH/NCI  Award  R01  CA1 14210. 


Dynamic  MREIT  Using  Sub-Milliamp  Currents 


M.  J.  Hamamura1,  L.  T.  Muftuler1,  O.  Birgul1,  O.  Nalcioglu1 

lrTu  &  Yuen  Center  for  Functional  Onco- Imaging,  University  of  California,  Irvine,  CA,  United  States 

Purpose 

In  Magnetic  Resonance  Electrical  Impedance  Tomography  (MREIT),  electrical  currents  are  injected  into  an  object  and  the  resulting  magnetic  flux  density 
distribution  measured  using  MRI.  These  MRI  measurements  are  then  used  to  reconstruct  the  conductivity  distribution  within  the  object.  Previous  MREIT  studies  have 
focused  primarily  on  reconstructing  static  conductivity  distributions.  However,  the  ability  to  detect  changes  in  conductivity  over  time  could  provide  additional 
diagnostic  information.  We  previously  reported  on  the  qualitative  monitoring  of  ion  diffusion  in  agarose  gel  over  four  time  points  using  MREIT  with  10mA  injected 
currents.  In  this  study,  we  perform  a  more  thorough  monitoring  using  sub-milliamp  injected  currents  more  appropriate  for  human  use. 


Methods 

For  the  test  phantom,  a  hollow  acrylic  disk  with  an  inner  diameter  of  7cm  and  thickness  of  1cm  was  filled  with  2% 
agarose  and  4mM  CUSO4.  Within  this  disk,  a  smaller  circular  region  of  12mm  diameter  was  filled  with  2%  agarose,  1% 
NaCl,  and  4mM  CuS04  (Figure  2a).  Over  time,  the  NaCl  diffused  from  the  smaller  region  into  the  background,  and  the 
conductivity  distribution  changed.  A  linear  relationship  between  conductivity  and  NaCl  concentration  (of  1%  or  less) 
was  found  after  performing  a  range  of  conductivity  measurements  using  the  4-electrode  technique.  The  plane  of  the 
disk  was  placed  perpendicular  to  the  main  static  MRI  field.  Four  copper  electrodes  each  6mm  wide  were  placed 
equidistantly  along  the  inner  acrylic  wall  and  used  to  inject  currents  into  the  interior  region. 

A  finite  alternating  current  pulse  waveform  with  an  amplitude  of  900uA  was  injected  into  the  phantom  and  the 
resulting  magnetic  flux  density  distribution  measured  using  a  modified  spin-echo  pulse  sequence  (Figure  2)  [Mikac  et  al, 
MRI  19:  845-856  (2001)].  The  scan  parameters  were  TR=500ms,  TE=60ms,  NEX=4,  Matrix=64X64,  FOV=10cm,  and 
single  slice  thickness  =  5mm.  Data  was  collected  for  two  different  current  injection  schemes  (in  pairs  of  electrodes 
directly  opposite  of  each  other)  and  used  simultaneously  in  conductivity  reconstruction.  To  reconstruct  the  conductivity 
distribution  using  the  MRI  measurements,  the  Sensitivity  Matrix  Method  was  utilized  [Birgul  et  al,  Phys  Med  Bio  48: 
3485-3504  (2003)]  in  which  the  relationship  between  conductivity  and  magnetic  flux  density  is  linearized  around  an 


Figure  1.  Pulse  sequence  used  in  MREIT 


initial  conductivity  (i.e.  uniform  distribution)  and  formulated  as  a  matrix  equation.  This  equation  is  then  solved  for  the  true  conductivity  distribution  using  Tikhonov 
regularization.  The  resulting  conductivity  can  then  be  substituted  back  into  the  linearized  equation  as  the  new,  updated  initial  condition,  and  the  process  iterated  to 
improve  the  reconstruction. 


Figure  2a  Figure  2b  Figure  2c  Figure  2d  Figure  2e  Figure  2f 


Figure  2g  Figure  2h  Figure  2i  Figure  2j  Figure  2k  Figure  21 


Figure  2.  (a)  Schematic  of  the  phantom;  (b)  Conductivity  after  20  minutes;  (c)  1  hour;  (d)  2  hours;  (e)  3  hours; 
(f)  4  hours;(g)  5  hours;  (h)  6  hours;  (i)  7  hours;  (j)  8  hours;  (k)  10  hours;  ( l )  12  hours 


Results 

Data  was  collected  at  various  time  points,  and  the  conductivity  distributions  reconstructed  (Figures  2b-l).  The 
resulting  images  clearly  show  a  change  in  the  conductivity  distribution  consistent  with  the  diffusion  of  NaCl  from 
the  higher  concentration  region  to  the  lower  concentration  region.  For  a  disk  of  radius  a  on  an  infinite  plane 
surface,  the  theoretical  concentration  C  at  the  center  of  the  disk  is  given  as  C=C0(l-exp(-a2/4Dt )).  The  peak 
conductivity  values  were  fitted  to  this  equation,  and  an  experimental  diffusion  constant  of  D  =  5. 8x1  O'6  cm2  sec'1 
obtained  (Figure  3).  This  value  is  smaller  than  the  previously  reported  measurement  of  1.4xl0'5  cm2  sec'1 
[Schantz  et  al,  Biochem  1:  658-663  (1962)].  The  slower  apparent  diffusion  in  this  phantom  can  in  part  be 
attributed  to  its  finite  shape  and  confined  volume. 

Discussion 

The  results  of  this  study  demonstrate  that  MREIT  can  monitor  changes  in  conductivity  over  time  using  sub- 
milliamp  injected  currents.  Validating  this  ability  is  a  necessary  step  towards  the  imaging  and  monitoring  of  in 
vivo  subjects  using  MREIT. 
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Purpose 

In  Magnetic  Resonance  Electrical  Impedance  Tomography  (MREIT)  electrical  currents  are  injected  into  an  object  and  the  resulting  magnetic  flux  density 
distribution  measured  using  MRI.  These  MRI  measurements  are  then  used  to  reconstruct  the  conductivity  distribution  within  the  object.  Many  of  the  reported  MREIT 
reconstruction  algorithms  utilize  numerical  calculation  of  the  magnetic  flux  density  for  a  given  conductivity  distribution  using  the  boundary  conditions  applied  to  the 
real  object.  This  corresponds  to  matching  the  position  of  the  electrodes  in  the  numerical  computation  to  that  of  the  actual  position  of  the  electrodes  on  the  object.  Near 
an  injecting  electrode,  there  exists  a  large  variation  in  the  magnetic  flux  density.  As  a  result,  any  misalignment  in  the  position  of  an  electrode  can  result  in  significant 
errors  when  calculating  the  difference  between  the  computed  magnetic  flux  density  distribution  and  the  MRI-measured  magnetic  flux  density  distribution.  Such  errors 
may  generate  artifacts  in  the  final  reconstructed  conductivity  distribution.  In  this  study,  we  investigate  various  correction  algorithms  to  reduce  these  artifacts. 

Method 

For  the  test  phantom,  a  hollow  acrylic  disk  with  an  inner  diameter  of  7cm  and  thickness  of  1cm  was  filled  with  2%  agarose,  0.1% 

NaCl,  and  4mM  CuS04.  Within  this  disk,  a  smaller  circular  region  of  12mm  diameter  was  filled  with  2%  agarose,  1%  NaCl,  and 
4mM  CuS04  (Figure  1).  The  conductivities  of  the  different  regions  were  measured  using  the  4-electrode  method  and  found  to  yield  a 
contrast  ratio  of  1  to  7.4.  The  plane  of  the  disk  was  placed  perpendicular  to  the  main  static  MRI  field.  Four  copper  electrodes  each 
6mm  wide  were  placed  equidistant  along  the  inner  acrylic  wall  and  used  to  inject  currents  into  the  interior  region. 

A  finite  alternating  current  pulse  waveform  with  an  amplitude  of  900uA  was  injected  into  the  phantom  and  the  resulting  magnetic 
flux  density  distribution  measured  using  a  modified  spin-echo  pulse  sequence  (Figure  2)  [ Mikac  et  al,  MRI  19:  845-856  (2001)].  The 
scan  parameters  were  TR=500ms,  TE=60ms,  NEX=4,  Matrix=64X64,  FOV=10cm,  and  single  slice  thickness  =  5mm.  Data  was 

collected  for  two  different  current  injection  schemes  (in  pairs  of  electrodes 
directly  opposite  of  each  other)  and  used  simultaneously  in  conductivity 
reconstruction.  To  reconstruct  the  conductivity  distribution  using  the  MRI 
measurements,  the  Sensitivity  Matrix  Method  was  utilized  [Birgul  et  al,  Phys 
Med  Bio  48:  3485-3504  (2003)]  in  which  the  relationship  between 
conductivity  and  magnetic  flux  density  is  linearized  around  an  initial 

conductivity  (i.e.  uniform  distribution)  and  formulated  as  a  matrix  equation.  This  equation  is  then  solved  for  the 
true  conductivity  distribution  using  Tikhonov  regularization.  The  resulting  conductivity  can  then  be  substituted 
back  into  the  linearized  equation  as  the  new,  updated  initial  condition,  and  the  process  iterated  to  improve  the 
reconstruction. 

Three  different  electrode  misalignment  correction  algorithms  can  be  implemented  during  reconstruction.  For 
the  first  algorithm  (MASK),  magnetic  flux  density  measurements  within  1cm  of  the  electrodes  were  discarded 
and  not  used  during  reconstruction.  For  the  second  algorithm  (SHIFT),  the  position  of  each  electrode  assigned 
during  numerical  computation  was  perturbed  to  find  the  best  location.  The  difference  between  the  MRI- 
measured  magnetic  flux  density  and  the  calculated  magnetic  flux  density  given  the  initial  condition  was 
minimized  as  a  function  of  electrode  position.  For  the  third  algorithm  (REG),  conductivity  perturbations  within 
lcm  of  the  electrodes  were  suppressed  by  a  factor  of  2  when  compared  to  the  rest  of  the  phantom  through 
increased  weighting  in  the  Tikhonov  regularization. 


Figure  2.  Pulse  sequence  used  in  MREIT 


Figure  1.  Schematic 
of  the  phantom 


Contrast  =  4.3555  Contrast  =  7.1568  Contrast  =  7.2036  Contrast  =  6.7546  Contrast  =  7.5216 


Figure  3a  Figure  3b  Figure  3c  Figure  3d  Figure  3e 

Figure  3.  (a)  Reconstructed  conductivity  distribution  using  no  correction;  (b)  MASK;  (c)  SHIFT;  (d)  REG;  (e)  all  3  algorithms 


Results 

Data  was  collected  using  a  4T  MRI  system.  Conductivity  distributions  were  reconstructed  for  various  combinations  of  electrode  correction  algorithms  using  5 
iterations  of  the  Sensitivity  Matrix  Method  (Figures  3a-e).  For  each  reconstructed  image,  the  contrast  between  the  high  and  low  conductivity  regions  was  calculated  by 
finding  the  maximum  conductivity  value  and  dividing  it  by  the  mean  conductivity  of  the  background.  The  background  conductivity  was  calculated  by  finding  the 
average  conductivity  of  the  phantom  excluding  a  2cm  diameter  disk  centered  on  the  high  conductivity  region. 

Discussion 

The  results  of  this  study  indicate  that  electrode  misalignment  affects  the  reconstructed  conductivity  distribution  throughout  the  phantom.  Regions  next  to  the 
electrodes  contain  artifacts  as  circled  in  Figure  3a,  while  objects  in  the  interior  region  suffer  from  diminished  contrast.  Each  of  the  proposed  correction  algorithms 
improves  the  reconstruction,  with  the  best  result  occurring  when  all  three  methods  are  applied. 
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Purpose 

In  Magnetic  Resonance  Electrical  Impedance  Tomography  (MREIT)  electrical  currents  are  injected  into  an  object 
and  the  resulting  magnetic  flux  density  distribution  measured  using  MRI.  These  MRI  measurements  are  then  used  to 
reconstruct  the  conductivity  distribution  within  the  object.  In  order  to  determine  the  conductivity  distribution  uniquely, 
data  from  at  least  two  different  current  distributions  satisfying  lli(x,y)>^2(x,y)  I  ^0  must  be  acquired  [Kwon  O.  et  al, 
IEEE  Trans  on  BME  49:  160-167  (2002)].  Typically,  two  electrodes  are  used  to  provide  one  current  distribution,  and 
two  additional  electrodes  are  used  to  provide  a  second  current  distribution,  for  a  total  of  four  electrodes  and  two  current 
injection  schemes.  However,  with  four  electrodes,  one  can  apply  up  to  six  different  current  injection  schemes  using 
different  pairs  of  electrodes  (Figure  1).  In  this  study,  we  assess  the  effects  of  utilizing  these  additional  current  injection 
schemes. 


Method 


For  the  test  phantom,  a  hollow  acrylic  disk  with  an  inner  diameter  of  7cm  and  thickness  of  1cm 
was  filled  with  2%  agarose,  0.1%  NaCl,  and  4mM  CUSO4.  Within  this  disk,  a  three  smaller  circular 
regions  each  1cm  in  diameter  were  filled  with  2%  agarose,  1%  NaCl,  and  4mM  CuS04  (Figure  3). 

The  conductivities  of  the  different  regions  were  measured  using  the  4-electrode  method  and  found  to 
yield  a  contrast  ratio  of  1  to  7.4.  The  plane  of  the  disk  was  placed  perpendicular  to  the  main  static 
MRI  field.  Four  copper  electrodes  each  6mm  wide  were  placed  equidistant  along  the  inner  acrylic 
wall  and  used  to  inject  currents  into  the  interior  region. 

A  finite  alternating  current  pulse  waveform  with  an  amplitude  of  900uA  was  injected  into  the 
phantom  and  the  resulting  magnetic  flux  density  distribution  measured  using  a  modified  spin-echo 
pulse  sequence  (Figure  2)  [Mikac  et  al,  MRI  19:  845-856  (2001)].  The  scan  parameters  were 
TR=500ms,  TE=60ms,  NEX=4,  Matrix=64X64,  FOV=10cm,  and  single  slice  thickness  =  5mm.  To 
reconstruct  the  conductivity  distribution  using  the  MRI  measurements,  the  Sensitivity  Matrix  Method 
was  utilized  [Birgul  et  al,  Phys  Med  Bio  48:  3485-3504  (2003)]  in  which  the  relationship  between 
conductivity  and  magnetic  flux  density  is  linearized  around  an  initial  conductivity  (i.e.  uniform 

distribution)  and  formulated  as  a  matrix  equation.  A  separate  equation  is  generated  for  each  current  injection  scheme,  and  various  combinations  of 
these  equations  can  be  combined  to  solve  for  the  conductivity  distribution.  The  solution  is  obtained  using  Tikhonov  regularization.  The  resulting 
conductivity  can  then  be  substituted  back  into  the  linearized  equation(s)  as  the  new,  updated  initial  condition,  and  the  process  iterated  to  improve  the 
reconstruction. 


Figure  1.  Multiple  Current 
Injection  Schemes 


Figure  2.  Pulse  sequence  used  in  MREIT 
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5.3522 

Region  C 

3.6669 

3.6327 

3.6439 

Table  1.  Peak  conductivities 


Figure  2.  Schematic  of  phantom 


Results 

Data  was  collected  for  each  of  the  6  current  injection  schemes. 
Conductivity  distributions  were  reconstructed  for  different 
combinations  of  data  sets  using  5  iterations  of  the  Sensitivity  Matrix 
Method  (Figures  4a-c).  For  each  reconstructed  image,  the  peak 
conductivity  in  each  of  the  three  inner  regions  A,  B,  and  C  was  found 
(relative  to  the  background  conductivity  of  1),  and  the  results 
compiled  in  Table  1.  A  closer  view  of  region  C  was  also  extracted 
(Figures  5a-c). 


Fig  4.  Reconstructed  conductivity  using  injection  schemes:  (a)  1,2;  (b)  1,2,3;  (c)  1-6 


Figure  5a  Figure  5b  Figure  5c 

Figure  5.  Conductivity  of  Region  Cfor  injection  schemes:  (a)  1,2;  (b)  1,2,3;  (c)  1-6 


Discussion 

The  results  indicate  that  reconstructing  periphery  regions  away 
from  the  injecting  electrodes  presents  difficulties,  regardless  of  the 
current  injection  schemes  used.  From  inspection  of  Figure  5, 
including  data  from  the  periphery  injection  schemes  (3-6  in  Figure  1) 
improves  the  overall  shape  of  object  C.  In  particular,  adding  scheme 
3  provides  the  largest  current  density  to  the  region  of  interest  and 
appears  to  best  improve  the  spatial  distribution.  However,  overall 
contrast  is  best  when  using  only  the  standard  orthogonal  injection 
schemes  (1  and  2).  Using  additional  injection  schemes  further 
reduces  contrast  throughout  the  object. 

Selecting  which  current  injection  schemes  to  use  requires 
balancing  overall  contrast  with  improved  spatial  resolution  in  the 
periphery  regions,  and  will  thus  depend  on  the  object  to  be  imaged 
and  the  desired  information.  The  increased  scan  time  required  to 
collect  data  from  addition  current  injection  schemes  must  also  be 
considered. 
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Puryose 

Several  studies  have  shown  that  the  electrical  impedance  of  malignant  tissues  is  significantly  different  from  those  of  normal  and  benign  tissues 
[1,2].  Therefore,  in-vivo  impedance  imaging  of  suspicious  lesions  has  the  potential  to  improve  the  sensitivity  and  specificity  of  detecting  malignant 
tumors.  MR-Electrical  Impedance  Tomography  (MREIT)  has  been  recently  introduced,  in  which  weak  electrical  currents  are  injected  into  the  tissue 
and  the  resulting  perturbations  in  magnetic  field  were  measured  using  phase  information  in  MR  images.  We  have  reported  our  preliminary  studies 
with  phantoms  as  well  as  two  in  vivo  experiments  [3].  Here,  we  present  the  results  of  MREIT  done  on  eleven  animals.  Parameters  like  variance  and 
mean  in  the  tumor  versus  the  rest  of  the  body  were  investigated.  The  goal  is  to  verify  potential  of  MR-EIT  to  aid  in  the  diagnosis  of  tumors. 


Methods 

Weak  electrical  currents  that  are  injected  into  an  object  generate  magnetic  fields,  the  z-component  of  which  induces  additional  phase  information 
in  MR  images.  If  a  modified  spin-echo  sequence  was  used  with  several  n  pulses  applied  during  the  zero-crossings  of  the  alternating  current,  the  phase 
shift  accumulates  across  these  n  pulses  and  is  given  in  the  final  image  as  cp(r)  =  4  y-N  bz(r)  /  co,  (y:  gyromagnetic  ratio;  N:  the  number  of  cycles  of 
injected  current;  bz(r):  the  amplitude  of  current-generated  magnetic  field  at  point  r;  co:  angular  frequency  of  the  injected  current).  Here  bz(r)  is 
calculated  from  the  phase  cp(r)  measurements.  We  have  implemented  an  iterative  reconstruction  with  Tikhonov  regularization  to  reconstruct  the 
conductivity  images  from  bz(r):  STAbz(r)  =  (STS+M)  Ao(r’).  Here  the  sensitivity  matrix  S  is  calculated  using  Finite  Element  Method;  Abz(r)  is  the 


change  in  magnetic  field  at  point  r  for  a  given  current  injection  scheme  resulting  from  a 
regularization  parameter  and  I  is  the  identity  matrix.  The  details  were  given  in  [4]. 

Data  were  collected  in  a  whole  body  4T  MRI  system  with  a  MRRS  console.  Eleven 
rats  were  imaged,  ten  of  which  were  bearing  malignant  tumors  that  were  either  R3230AC 
tumor  grafts  or  induced  by  the  carcinogen  ENU.  Animals  were  anesthetized  prior  to 
imaging  and  all  procedures  were  approved  by  the  IACUC.  Two  data  sets  were  discarded 
due  to  severe  motion  artifacts.  Structural  images  were  collected  using  T2  weighted  SE 
sequence  prior  to  MREIT  images.  The  data  matrix  was  128X128,  FOV  =  10cm,  slice 
thickness  =  4mm,  with  2mm  gap.  TR  =  3s,  TE  =  50ms  and  NEX  =  2  were  used.  MREIT 
images  were  collected  using  the  outlined  pulse  sequence  with  TR= 1000ms,  TE=30ms,  and 
NEX=2,  64X64  data  matrix,  FOV  =  10cm,  slice  thickness  =  4mm  with  2mm  gap.  Two 
cycles  of  100Hz  current  with  1mA  rms  was  applied  sequentially  through  different  pairs  of 
four  electrodes,  generating  six  different  current  profiles.  Data  were  collected  with  both  + 
polarities  of  the  currents  to  eliminate  phase  accumulation  from  other  sources. 

Results  and  Discussion 

Fig.l  shows  T2  weighted  MRI  and  MREIT  images  of  two  rats.  Tumor  areas  show 
increased  conductivity  depicted  with  yellow-red  colors.  On  the  MREIT  images  of  the 
animals,  separate  ROIs  (region  of  interest)  were  drawn  over  the  tumor  region  and  the  rest 
of  the  body  and  the  mean  conductivity  values  in  these  ROIs  were  calculated.  Since  MREIT 
yields  relative  conductivity  values,  the  conductivity  ratio  atum0r/c>body  was  calculated  for 
each  animal  and  the  graph  is  shown  in  Fig.2.  It  was  found  that  the  average  of  these 
conductivity  ratios  was  2.17.  ROIs  were  drawn  manually  based  on  the  tumor  seen  in  MRI 
T2  images.  We  have  also  calculated  the  ratio  of  standard  deviation  to  mean  conductivity  in 
each  tumor  region,  which  may  be  an  indication  of  conductivity  heterogeneity  inside  the 
tumor  volume,  rather  than  SNR  (Fig. 3).  As  seen  from  these  figures,  consistent  results  were 
obtained  from  these  eight  animals.  Average  conductivity  increased  by  roughly  2.2  times  in 
the  tumor  compared  to  the  rest  of  the  body.  The  conductivity  varied  typically  between 
10%  and  20%  within  the  tumor.  In  two  of  the  rats,  we  have  also  collected  MRIET  data 
with  rms  currents  of  0.5,  1  and  2mA.  Highly  consistent  results  were  seen  especially 
between  1  and  2mA  cases.  For  example,  when  MREIT  images  of  1mA  and  2mA  cases 
were  subtracted,  the  mean  of  the  residual  was  only  1.3%  of  the  mean  conductivity  in  the 
whole  slice.  In  the  case  of  0.5mA  vs  1mA,  the  mean  of  the  residual  was  12%. 

In  this  study,  it  has  been  demonstrated  that  MRI  based  impedance  imaging  has  the 
potential  to  investigate  malignant  tumors  in  vivo.  At  this  stage  the  technique  may  not  be 
used  as  a  standalone  diagnostic  tool;  but  it  could  provide  useful  information  in 
characterizing  the  tumor,  once  a  suspicious  lesion  is  detected  by  other  methods  to  improve 
specificity.  Physiological  and  structural  changes  in  tumors  that  lead  to  such  changes  in 
conductivity  will  be  investigated  in  future  studies.  Note  that  good  quality  MREIT  images 
levels. 


change  Aa(r’)  in  the  conductivity  at  point  r\  X  is  the 


Fig.l.  Results  from  two  animals  are  illustrated.  T2 
weighted  scans  are  displayed  above  and  corresponding 
impedance  (MREIT)  images  are  depicted  in  color  right 
below.  Tumor  areas  are  marked  with  red  circles. 


Fig.2.  The  ratio  of  mean  conductivity  in  the  tumor 
versus  the  rest  of  the  body  over  eight  in  vivo  studies. 


Fig.3.  The  ratio  of  standard  deviation  of  conductivity 
versus  its  mean  in  the  tumor  in  eight  in  vivo  studies. 

were  collected  with  biologically  safe  electrical  current 
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Purpose 

Several  researchers  have  demonstrated  that  the  electrical  properties,  specifically  the  impedance  of  malignant  tissues  is  significantly  different  from 
those  of  normal  and  benign  tissues  [Malich  et  al,  Clin.Rad.,5 6:278-283,  2001;  Silva  et  al.  Med.  Biol.  Eng.  and  Comp.,  38:26-30,  2000].  Therefore,  in 
vivo  impedance  imaging  of  suspicious  lesions  could  aid  in  improving  the  sensitivity  and  specificity  of  detecting  malignant  tumors.  MRI  based 
impedance  imaging  is  a  novel  method,  in  which  weak  electrical  currents  are  injected  into  the  tissue  and  the  resulting  perturbations  in  magnetic  field 
were  measured  using  phase  information  in  MR  images.  We  have  reported  our  preliminary  studies  with  phantoms  as  well  as  some  in  vivo  experiments. 
However,  in  those  studies,  reconstruction  was  done  with  the  non-iterative  sensitivity  matrix  approach,  which  assumed  a  linear  relationship  between 
the  measured  field  perturbations  and  relative  conductivity  distribution.  For  large  conductivity  variations,  this  linearity  assumption  failed  and  the 
resulting  images  underestimated  the  conductivity  variations.  This  resulted  in  poorer  spatial  specificity  and  sensitivity.  We  refined  our  reconstruction 
technique  with  an  iterative  approach  to  account  for  the  nonlinearities.  The  goal  is  to  verify  potential  of  MR-EIT  to  aid  in  the  diagnosis  of  tumors. 


Methods 

When  a  sinusoidal  current  is  injected  into  an  object,  the  z-component  of  the  resulting  magnetic 
fields  induce  additional  phase  information  in  MR  images.  A  modified  spin-echo  sequence  was 
used,  where  several  n  pulses  were  applied  during  the  zero-crossings  of  the  sinusoidal  current 
(Fig.l)  [Mikac  et  al,  MRI  19:  845  856  (2001)].  This  phase  shift  accumulates  across  these  n  pulses 
and  is  given  in  the  final  image  as  cp(r)  =  4-y-Nbz(r)  /  co,  where  y  is  the  gyro  magnetic  ratio,  N  the 
number  of  cycles  of  injected  current,  bz(r)  the  amplitude  of  z-component  of  current- generated 
magnetic  field  at  point  r,  and  co  the  angular  frequency  of  the  injected  current.  Therefore,  once  the 
phase,  cp(r)  is  measured,  bz(r)  can  be  calculated.  In  the  original  approach,  a  linear  approximation 
AB(r)  =  S(r,r’)  Aa(r’)  was  assumed  to  reconstruct  conductivity  image.  Here  AB(r)  is  the  change  in 
magnetic  field  at  point  r  for  a  given  current  injection  scheme  resulting  from  a  change  Aa(r’)  in  the 
conductivity  at  point  r\  S  is  calculated  using  Finite  Element  Method  (FEM).  Recently,  we  have 
adopted  an  iterative  reconstruction  with  Tikhonov  regularization,  which  improved  the  accuracy  as  0  5  i°  Timers)  20  25  30 

well  as  spatial  specificity  of  the  MREIT  images.  Details  of  this  method  and  its  assessment  with  Fig.  1 .  Pulse  sequence  used  in  MR-EIT 

phantom  studies  were  presented  in  another  manuscript  that  is  being  submitted  concurrently.  In 

summary,  the  matrix  equation  was  modified  with  the  addition  of  a  regularization  parameter ,  X,  and  became  (STAB  =  (STS+AI)  Ag),  where  I  is  the 
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identity  matrix.  This  equation  was  solved  for  different  values  of  X  using  conjugate  gradient  method  and  the  optimum  value  of  X  was  found.  The 
conductivity  distribution  that  is  calculated  after  each  iteration  was  assigned  as  the  initial  value  for  the  next  iteration  and  the  iterations  are  repeated 
until  the  difference  between  the  results  of  two  consecutive  iterations  are  below  a  predefined  threshold.  Current  was  injected  from  four  electrodes 
placed  across  the  animal’s  body  and  six  different  current  distributions  were  generated  by  using  different  combinations  of  electrode  pairs.  This  scheme 
also  improved  the  uniformity  and  sensitivity  of  the  MREIT  images. 

Results 

Data  were  collected  in  a  whole  body  4T  MRI  system  with  a 
MRRS  console.  Two  rats  were  imaged,  which  were  bearing  malignant 
tumors  induced  by  the  carcinogen  ENU.  Animals  were  anesthetized 
prior  to  imaging.  All  procedures  were  approved  by  the  IACUC. 

Structural  images  were  collected  using  SE  sequence  prior  to  MREIT 
images.  The  data  matrix  was  128X128,  FOV  =  10cm,  slice  thickness 
=  4mm,  with  2mm  gap.  TR  =  3s,  TE  =  50ms  and  NEX  =  2  were  used. 

MREIT  images  were  collected  using  the  previously  outlined  pulse 
sequence  with  TR=500ms,  TE=30ms,  and  NEX=12,  64X64  data 
matrix,  FOV  =  10cm,  slice  thickness  =  4mm  with  2mm  gap.  2  cycles 
of  100Hz  current  with  2mA  rms  was  applied  through  different  pairs  of 
four  electrodes.  The  resulting  images  clearly  show  the  higher 
conductivity  regions.  Especially  in  animal- 1,  higher  conductivity  map 
of  the  tumor  precisely  overlaps  with  the  one  seen  in  the  structural 
image.  The  mean  conductivity  in  the  tumor  region  was  2.8  times 
higher  than  the  mean  conductivity  in  the  rest  of  the  body  for  animal- 1 . 

Discussion 


Fig. 2.  Structural  (top)  and  MREIT  (bottom)  images  of  two  animals.  Left  and  right 
columns  show  the  images  of  animal- 1  and  animal-2,  respectively.  Tumor  areas  are 
circled  with  red  lines.  Bright  objects  outside  the  animals  were  markers  to  identify 
exact  location  of  electrodes. 


In  this  study,  it  has  been  demonstrated  that  MRI  based  impedance 
imaging  has  the  potential  to  investigate  malignant  tumors  in  vivo. 

Improvements  in  current  injection  and  utilization  of  iterative  reconstruction  with  regularization  significantly  improved  the  accuracy  and  spatial 
specificity  of  MREIT.  In  our  experiments  it  was  also  observed  that  precise  localization  of  electrode  placement  was  critical  to  define  boundary 
conditions  accurately,  which  in  turn  affected  the  quality  of  MREIT  images.  In  the  presented  study,  2D  FEM  formulation  was  used.  Currently  3D 
FEM  is  being  developed  to  improve  the  accuracy  of  MREIT  images  further. 
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Purpose 

It  has  been  reported  that  the  electrical  impedance  of  malignancies  is  20-40  times  lower  than  healthy  tissues  and  benign  formations  [ Malich  A.  et  al,  Eur.  Radiol.  10: 
1555-1561  (2000)].  Therefore,  in  vivo  impedance  imaging  of  suspicious  lesions  could  aid  in  the  diagnosis  of  malignant  tumors.  In  MR-EIT,  electrical  currents  are 
injected  into  an  object  and  the  resulting  magnetic  field  perturbations  are  measured  using  MRI.  These  measurements  are  then  used  to  reconstruct  the  conductivity 
distribution  within  the  object.  Previous  studies  on  MR-EIT  have  focused  on  reconstructing  static  conductivity  distributions.  However,  the  ability  to  detect  changes  in 
conductivity  over  time  could  provide  additional  diagnostic  information,  such  as  in  monitoring  tumor  growth.  In  this  study,  we  assess  the  ability  of  MR-EIT  to  detect 
these  changes. 


Methods 

Sinusoidal  current  is  injected  into  an  object  and  the  resulting  magnetic  fields  are  measured  using 
a  modified  spin-echo  sequence  (Figure  1)  [Mikac  U.  et  al,  MRI  19:  845-856  (2001 J].  The  component 
of  current- generated  magnetic  field  parallel  to  the  main  static  field  (z-component)  introduces  a  phase 
shift.  By  synchronizing  successive  7t  pulses  to  half  cycles  of  the  current,  this  phase  shift  accumulates 
and  is  given  in  the  final  image  as  (p(r)  =  4-y-N-b(r)  /  co,  where  y  is  the  gyromagnetic  ratio,  N  the 
number  of  cycles  of  injected  current,  b(r)  the  amplitude  of  z-component  current- generated  magnetic 
field  at  point  r,  and  co  the  angular  frequency  of  the  injected  current.  Hence,  measurement  of  this 
phase  shift  allows  for  calculation  of  the  (z-component)  magnetic  field  distribution. 

To  reconstruct  the  conductivity  distribution  from  the  magnetic  field  distribution,  a  linear 
approximation  AB(r)  =  S(r,r’)  Aa(r’)  is  assumed,  where  AB(r)  is  the  change  in  magnetic  field  at 
point  r  for  a  given  current  injection  scheme  resulting  from  a  change  Aa(r’)  in  the  conductivity  at 
point  r\  The  matrix  component  Sy  is  the  change  in  magnetic  field  dBi  of  element  i  with  respect  to  a 
change  in  the  conductivity  Ogj  of  element  j.  An  initial  conductivity  distribution  G^aiis  assumed  (e.g. 
uniform  conductivity),  the  conductivity  of  a  given  element  j  perturbed  by  Agj,  the  resulting  AB 
calculated  using  the  Finite  Element  Method  (FEM),  and  the  matrix  components  approximated  as  Sy  = 

ABi  /Acjj.  The  linear  approximation  can  be  rewritten  as  (Bfinai  -  Binitiai )  =  S  (Ofjnai  -  (Uritiai),  where  o^tiai 
is  the  assumed  initial  (uniform)  conductivity  distribution,  Binitial  the  FEM  calculated  magnetic  field 
distribution  given  Ojnjtiai ,  Bfinai  the  MRI  measured  magnetic  field  distribution,  and  Gfinai  the  actual 
conductivity  distribution.  The  equation  is  solved  for  Gfinai  using  Tikhonov  regularization.  This  Gfinai 
then  substituted  back  into  the  linear  approximation  as  the  new,  updated  Ginitiai,  and  the  process  is  iterated  until  the  change  in  conductivity  between  successive  iterations  is 
below  some  predefined  threshold.  Details  of  this  reconstruction  method  can  be  found  in  another  abstract  submitted  concurrently. 

The  previously  outlined  method  was  used  to  measure  the  conductivity  distribution  of  an  agarose  gel  phantom.  For  the  phantom,  a  hollow  acrylic  disk  with  an  inner 
diameter  of  7cm  and  thickness  of  1cm  was  filled  with  1%  (g/lOOmL)  agarose  and  1%  NaCl.  Within  this  disk,  a  smaller  circular  region  of  12mm  diameter  was  filled 
with  1%  agarose  and  20%  NaCl  (Figure  2a).  Over  time,  NaCl  diffused  from  the  region  of  higher  concentration  to  the  region  of  lower  concentration.  Conductivity  is 
assumed  proportional  to  NaCl  concentration,  hence  the  conductivity  distribution  changed  as  a  result  of  this  diffusion.  The  plane  of  the  disk  was  placed  perpendicular  to 
the  z-axis.  Four  recessed  electrodes  each  3mm  wide  were  placed  equidistant  along  the  circular  acrylic  wall  and  used  to  inject  currents  into  the  interior  region.  Data  was 
collected  for  two  different  current  injection  schemes  (in  pairs  of  electrodes  directly  opposite  of  each  other)  and  used  simultaneously  in  conductivity  reconstruction. 


Time  (ms) 

Figure  1 .  Pulse  sequence  used  in  MR-EIT 


Results 

2  cycles  of  10mA  (rms)  100Hz  current  were  injected  into  the  phantom  using  the  previously 
outlined  pulse  sequence  with  the  parameters  TR=500ms,  TE=30ms,  and  NEX  =  4.  Scans  were  taken 
20  minutes,  1  hour,  6  hours,  and  24  hours  after  the  creation  of  the  phantom.  The  z-component 
current- generated  magnetic  field  distributions  were  calculated  from  the  resulting  data  and  the 
conductivity  distributions  computed  (Figures  2d-e).  The  resulting  images  clearly  show  a  change  in 
the  conductivity  distribution  consistent  with  the  diffusion  of  NaCl  from  the  higher  concentration 
region  to  the  lower  concentration  region.  Over  time,  the  higher  conductivity  region  broadens,  and  the 
border  between  the  initial  regions  becomes  less  distinct  (Figure  3). 

Discussion 

In  this  study,  we  have  shown  that  MR-EIT  can  be  used  to  monitor  changes  in  conductivity  over 
time.  While  we  assume  that  conductivity  is  directly  proportional  to  NaCl  concentration  for  this 
phantom,  future  studies  will  verify  this  by  correlating  this  method  with  Sodium  MRI.  We  also  plan  to 
monitor  malignancies  in  live  animals  by  performing  MR-EIT  measurements  over  the  span  of  several 
weeks. 
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Figure  3.  Profile  taken  across  the  dotted  red  line  of  Figure  2a 


